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13.  ABSTRACT 


1  I  2  ■  SPONSORING  MILITARY  ACTIVITY 


Air  Force  Flight  Dynamics  Laboratory 
Wright-Pntterson  AFB,  OH  45433 


Sonic  fatigue  testing  at  elevated  temperatures  has  become  increasingly 
important  as  aerospace  vehicle  structures  are  designed  for  the  combined  environments. 
Work  is  described  in  four  areas  related  to  thermal/acoustic  testing;  acoustic  simu¬ 
lation,  thermal  simulation,  specimen  mounting  affects,  and  instrumentation  and 
measurements. 

Acoustic  field  studies  considered  the  directional  properties  of  various 
fields,  as  well  as  coupling  of  the  acoustic  fields  with  structural  specimens. 

Thermal  environment  studies  dealt  with  predicting  temperatures  in  a  heated  struc- 
tural  specimen.  The  design  of  heating  systems  was  also  discussed.  In  mecimen 
mounting  effects,  equations  were  presented  to  show  how  thermal  loads  entt  •  into  the 
general  dynamic  equations.  Panels  with  free  edges  and  panels  with  fixed  adges  were 
treated  in  detail.  Measurement  methods  related  to  sonic  fatigue  testing  at  elevated 
temperatures  were  surveyed  for  availability  of  devices  which  operate  in  the  com¬ 
bined  environment,  / 

Results  presented  here  indicate  that  good  estimates  of  specimen  response 
can  be  made  for  basic  specimen  geometries,  mounting  conditions  and  orientations 
in  the  excitation  field.  Methods  are  suggested  for  extending  the  results  to  more 
complex  cases. 
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ABSTRACT 


Sonic  fatigue  testing  at  elevated  temperatures  has  become  increasingly  important  as 
aerospace  vehicle  structures  are  designed  for  the  combined  environments.  Work  is  described 
in  four  areas  related  to  thermal/acoustic  testing:  acoustic  simulation,  thermal  simulation, 
specimen  mounting  effects,  and  instrumentation  and  measurements. 

Acoustic  field  studies  considered  the  directional  properties  of  various  fields,  as  well  as 
coupling  f'fthe  acoustic  fields  with  structural  specimens.  Thermal  environment  studies  dep't 
with  predicting  temperatures  in  a  heated  structural  specimen.  The  design  of  heating  systems 
was  also  discussed.  In  specimen  mounting  effects,  equations  were  presented  to  show  how 
thermal  loads  enter  into  the  general  dynamic  equations.  Panels  with  free  edges  and  panels  with 
fixed  edges  were  treated  in  detail.  Measurement  methods  related  to  scnic  fatigue  testing  at 
elevated  temperatures  were  surveyed  for  availability  of  devices  which  operate  in  the  combined 
environment. 

Results  presented  here  indicate  that  good  estimates  of  specimen  response  can  be  made 
for  basic  specimen  geometries,  mounting  conditions,  and  orientations  in  the  excitation  field. 
Methods  are  suggested  for  extending  the  results  to  more  complex  cases. 
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SECTION  1 
INTRODUCTION 


Sonic  fatigue  testing  at  elevated  temperatures  has  become  a  matter  of  considerable 
importance  in  recent  years.  The  application  of  new  materials  and  advanced  fabrication 
methods  to  high  performance  aircraft  has  brought  a  need  for  evaluation  of  the  structural 
elements  in  the  combined  high  intensity  noise/high  temperature  environment. 

The  study  effort' reported  here  considered  all  aspects  of  high  temperature,  sonic 
fatigue  testing.  Characteristics  of  the  acoustic  environment  are  discussed,  along  with  coup¬ 
ling  of  the  acoustic  excitation  with  the  test  specimen.  Thermally  and  dynamically  induced 
stresses  are  considered  in  detail,  including  the  effect  of  the  specimen  mounting  on  the 
resultant  specimen  response.  Information  is  reported.on  methods  and  devices  for  measuring 
ali  the  prameters  that  enter  into  combined  thermal/acoustic  testing. 

Following  a  brief  summary  of  the  above  findings,  recommendations  are  made  for 
applications  to  combined  environment  tests,  and  areas  for  further  study  are  indicated. 


SECTION  2. 

ACOUSTIC  ENVIRONMENT  SIMULATION 

Sonic  fatigue  tests  are  usually  conducted  in  facilities  that  produce  either  progressive 
wave  or  reverberant  (diffuse)  acoustic  fields.  One  phase  of  this  program  was  to  investigate 
the  use  of  such  facilities  for  sonic  fatigue  tests  at  elevated  temperatures.  Since' the  purpose 
of  a  sonic  fatigue  test  is  to  determine  the  ability  of  a  structural  specimen  to  withstand  acoui 
tically  induced  stresses,  a  study  of  acoustic  environment  simulation  should  also  consider  the 
coupling  between  the  acoustic  field  and  the  specimen.  Accordingly,  this  section  first  con-: 
siders  mathematical  representations  for  progressive  wave  and  reverberant  fields.  Then,  dtliis?. 
ing  the  dynamic  properties  of  a  structural  specimen,  a  method  is  developed  for  estimating 
the  response  of  a  structure  to  an  acoustic  field.  The  effects  of  temperature  on  the  acoustic 
field  in  a  representative  test  application  are  then  discussed.' 

2.1  ACOUSTIC  FIELDS  •  : 

The  analysis  of  the  progressive  wave  and  reverberant  fields  follows  the  usual  assumpr  .  ~ 
tions  of  elementary  acoustics.  The  sound  pressure  and  particle  velocities  are  assumed  so  be 
small  compared  to  the  static  pressure  and  sound  veloci fy.s’especfiyeiy .  Attenuation  in  the  • 
medium  will  be  neglected.  ^ 

-  '•  '-V**-.  *  .  •  ‘ 

2.1.1  Progressive  Wave  Fields 

A  commonly  used  sonic  fatigue  test  enclosure  is  the  progressive  woyedube,  consisting  . 
of  a  constant  cross-section  duct  with  an  acoustic  source  at  one  end,  and  sfim^proYisian 
for  mounting  a  structural  specimen  in  the  wall  of  the  duct.  The  acoustic  pmssures  in  such- 
an  enclosure  are  considered  to  be  uniform  across  any  cross-section-so  thaL.t.he  progrer-sh 
wave  field  is  a  one-dimensional  acoustic  field.'  Therefore,  the  one"dimensi(Sit5i;^a.yt>;c-“.' 
equation  . 

d2 p(x,  t)  _  _L  9^p^x’  ^ 

3x^  c2  at2  ;*“[)- 


can  be  used  to  describe  the  progressive,  wave  field.  Tile  solution  of  equation  (2.1 )  depends  ft- 
on  the  boundary  conditions  at  the  end  of  the  enclosure.  For  the  arrangement  shown  in 
Figure  1  for  example,  the  solution  for  a  perfectly  anechoic  termination  (i.e.,  no  presf,tve 
wave  reflected  from  the  end)  is  of  the  form 

jtot  +  tJi)  .  . 

p(x,  t)  P0e  *  (2.2) 

This  solution  is  for  a  simple  harmonic  excitation  at  an  angular  frequency  Co  and, re  presents 
a  pressure  wave  traveling  in  the  -x  direction.  The  pressure  amplitude  envelope  is  equal  to 
P0  and  is  constant  along  the  length  of  the  tube,  while  the  phase  of  the  pressure  between  . 
two  points  varies  linearly  with  the  distance  between  these  points. 

Other  types  of  terminations  will  yield  different  pressure  functions  for  the  progressive 
wave  enclosure.  If  the  anechoic  termination  were  replaced  by  an  infinitely  rigid  closure  at 
the  x=0  plane,  then  a  pressure  wave  is  reflected  back  toward  the  source,  and  the  pressure 
in  the  enclosure  is  given  by2 


P(x,  t)  = 


sin  kC 


2ttx  .* 
cos  “  ejwt 


2 


Neglecting  the: time-varying  component  ofthe  pressure,  we  note  that. the  envelope  of 
the  magnitude  of  the  standing  pressure  wave  is'a  cosine  function,  The  maximum  value  of 
the  pressure  envelope  Is  a  function  of  the  length  of  the  enclosure,  if  and  the  frequency  at 
which  the  test  section  is  excited.  The  frequency  dependence  is  contained  in  the  wave 
number  k,  defined  as  (vk  Vv  v..  ' 


Equation  (2.3)  does  not  consider  any  losses;  and  therefore  suggests  that  the  pressure  be¬ 
comes  infinite  at  frequencies  and  enclosure  lengths  such  that 

k§  -  mt,  n  =  1 , 2,  3, , . ^‘h:  '  (2.5) 

In  practice,  some  damping  is  always  present  to  limit  thq  pressure  to  a  finite  value,  although 
the  standing  wave  rqtio  is  a  maximum  when  the  conditions  of  equation  (2.5)  are  met. 

'  -  A  third  form  of  the  pressure  function  results  if  the.  enclosure  is  allowed  to  open 
directly  into  free  space  at  the  x-Q  plane.  In  this  case,  there  will  he  a  pressure  antinode  at 
the  interface,  with  tire  pressure  distribution  in  the  duct  being 

o  2mi.  i,  ,*  -  . 

p(x,  t)  =  7~  .sin  ~ r  e}<Jt  (2.6) 

-  ’  *■.  COS  KJc  A 

Comments  relative  to  the  peak  value  of  pressure  for  the  rigid  termination  case  are  also 
applicable  to  the  enclosure  with  an  open  end,  except  that  pressure  maxima  occur  at 

•'  '..k2  =  “,h=l,3,5  .  (2.7) 

_ 

The  simple  trigonometric  functions  in  equations  (2.2),  (2.3),  and  (2.6)  are  convenient 
to  describe  the  pressure  distributions  in  progressive  wave  test  sections,  since  the  test  speci¬ 
mens  are  mounted  in  the  wall  of  the  test  section.  The  frequency  range  encompassed  by 
sonic  fatigue  tests,  combined  with  the  cross-sectional  areas  likely  to  be  used  for  progressive 
wave  sections,  may  yield  conditions  such  that  the  acoustic  waves  do  not  travel  along  the,  • . 
enclosure  with  plane  fronts.  This  factor  is  not  important  for  the  present  discussion,  since 
the  pressure  distribution  along  the  surface  of  the  enclosure  is  still  described  by  the  trig¬ 
onometric  functions.  Still  another  effect  may  be  encountered  at  various  combinations  of 
test  frequency  and  enclosure  size.  An  acoustic  field  may  be  established  somewhere  between 
that  which  would  '  e  expected  for  an  anechoic  termination  and  an  abrupt  termination.  In 
such  cases,  a  combination  of  trigonometric  functions  may  be  needed  to  represent  the 
acoustic  field.  Liter  in  this  section,  wc  will  utilize  trigonometric  representations  of  acoustic 
fields  for  estimating  structural  response  to  the  fields. 

2. 1 .2  Reverberant  Fields 

Another  common  type  of  acoustic  test  enclosure  is  a  large  room  that  is  characterized 
by  low  acoustic  absorption  at  the  internal  surfaces  of  the  room.  In  such  a  room,  a  large  . 
percentage  of  the  acoustic  energy  incident  upon  a  given  surface  is  reflected  from  that 
surface.  The  resulting  reverberant  (diffuse)  acoustic  field  has  uniform  energy  density 
throughout  the  volume  of  the  enclosure,  and  all  directions  of  propagation  are  equally 
probable. 
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An  analysis  of  the  acoustic  field  in  a  reverberant  room  begins  with  the  general  wave 
.  equation,  although  the  equation  must  now  account  for  propagation  in  three  directions.  In 
its 'general  form,  the  wave  equation  can  be  written  as: 


2  - 


VP 


1  a2p 


(18) 


where  the  V2  operator  indicates  the  second  partial  derivatives  with  respect  to  the  x,  y,  and 
z  directions.  Solutions  which  satisfy  equation  (2.8)  can  be  of  the  form: 

jwt 


p(x,  y,  z,  t)  =  P0(cos  kxx)(cos  kyy)(cos  kzz)e 
provided  the  constants  kx,  ky(  and  kz  satisfy 


-  k  =  7=  /k2+k2+k 2 


(2.9) 


(2.10) 


These  constants  are  also  related  to  “component"  wavelengths,  the  wavelengths  associated 
with  the  three  coordinate  axes,  by 


2ir 


=  nx* 


(2.11a) 


7  ny7r 


(2.11b) 


•  2n 

h 


(2.ilc) 


■  I 


The  characteristic  frequencies  in  the  room  are  restricted  to  values  such  that” 

/«  \  2 


■ 


7 


(2.12) 


This  analysis -may  appear-  to  contradict  the  previously  stated  reverberant  field  condi¬ 
tions  of  uniform  energy  density  and  absence  of  a  preferred  direction  of  propagation.  Equa¬ 
tion  (2.9)  does,  in  fact,  show  that  the  pressure  waves  in  a  reverberant  enclosure,  at  a  normal 
mode  frequency  of  the  enclosure,  are  highly  directional.  Tlie  acoustic  energy  density,  which 
is  proportional  to  the  pressure  squared,  will  also  vary  widely  in  a  room  excited  at  a  single 
normal  mode  frequency.  Therefore,  in  order  to  obtain  a  diffuse  field,  the  room  must  be 
excited  e‘  4  number  of  normal  mode  frequencies.  Stated  differently,  a  diffuse  field  can  be 
established  in  a  reverberant  room  only  if  the  room  is  excited  by  a  signal  that  has  energy 
distributed  in  a  band  of  frequencies;  e.g.,  by  a  random  signal.: 

The  exact  number  of  modos  that  are  excited  in  a  given  band  of  frequencies  can  bo 
determined  from  equation  (2.12).  As  an  example,  consider  a  room  with  dimensions  £X=3G, 
8y=>40,  and  feet.  The  (l,l,l)modc  (ux=ny=nznl  )is  seen  to  be  25.2  Hz»  the  (1,1,2) 
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mode  is  3 1 .5  Hz,  and  so  on.  Each  of  these  normal  mode  frequencies  can  be  considered  to  be 
a  vector  in  a  frequency  space  with  orthogonal  axes  fx>  fy,  and  fz.  The  total  number  of  modes 
below  a  given  frequency,  f  is  then  approximated  by 


N  = 


4tt 

I? 


VP 


(2.13) 


where  V  is  the  volume  of  the  room.  Differentiation  of  equation  '2.13)  yields  an  expression 
for  the  number  of  modes,  AN,  within  a  band  of  frequencies  of  Af  centered  at  a  frequency  f: 


4ttV  o 
AN  = - f2 


Af 


(2.14) 


A  minimum  of  10  modes  should  be  excited  within  a  given  band  in  order  to  achieve  uniform 
energy  density  throughout  the  volume  of  the  room,  and  thereby  establish  a  reasonably  dif¬ 
fuse  field  in  the  room 

The  present  study  is  concerned  with  sonic  fatigue,  so  we  need  some  means  to  relate  the 
above  properties  of  the  test  room  to  the  response  of  the  structure  undergoing  sonic  fatigue 
testing.  The  test  enclosure/structural  response  dependence  begins  with  the  principle  that 
most  of  the  energy  acceptance  of  a  structure  occurs  near  the  normal  mode  frequencies  of 
the  structure.  Just  how  near  the  excitation  frequency  must  be  to  the  normal  mode  fre¬ 
quency  is  determined  by  the  damping  of  the  structure.  Panel  structures  of  the  type 
commonly  used  in  aerospace  structures  have  low  damping,  typically  f  <0.02.  The  bandwidth 
of  a  normal  mode  frequency,  Af,  for  a  panel  with  damping  ratio  f,  is  approximated  from 


Af  =  2ff 

This  expressi  combined  with  equation  (2.14),  yielding 
8tt  Vf  f3 


N  =  ■ 


,3 


(2.15) 


(2.16) 


Equation  (2.16)  thus  shows  the  number  of  room  inodes  that  will  be  excited  within  a 
structural  resonant  bandwidth,  for  a  resonance  with  a  center  frequency  f  and  damping  ratio 
f.  If  tire  previously  stated  criterion  of  10  modes  is  applied,  equation  (2.16)  may  be  solved 
for  V,  the  room  volume.  Figure  2  presents  a  family  of  curves,  for  various  values  of  (\  for 
the  lowest  frequency  at  which  10  modes  are  excited  in  a  given  room  volume.  At  frequencies 
below  this  minimum  frequency,  the  acoustic  field  will  become  more  directional  and  will 
have  properties  sii ir  to  the  progressive  wave  fields.  Above  the  minimum  frequency,  the 
acoustic  field  is  diffuse,  with  a  corresponding  absence  of  directional  properties. 

Pressure  distributions  in  a  large  room  excited  by  a  random  source  cannot  be  described 
with  exact  expressions,  as  was  the  case  for  simple  harmonic  sources.  Hie  pressures  must  in¬ 
stead  be  discussed  in  terms  of  their  power  spectra  or  correlation  functions.  For  example,  the 
normalized  cross-power  spectrum  for  u  bandwidth-limited  random  accoustic  field  is  given 


by 
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:Pi’  Pj  (ri‘  rj’  X)  = 


sm 


J2 


2jt  (r  j  ~-rj) 

T 


(2.17) 
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where  >  is  the  wavelength  associated,  with  the  upper  frequency  of  the  band,  and  the  indices 
i  and  j  represent  locations  in  the  field.  This  expression  is  valid  for  a  diffuse  field  in  which  a 
large  number  of  room  modes  are  excited  within  the  bandwidth  of  the  structure.  The  value 
of  the  function  depends  on  the  separation  between  two  points  in  the  field  rather  than  the 
exact  location  at  which  the  function  is  measured.  If  the  size  of  the  enclosure  is  such  that 
only  a  very  few  normal  modes  of  the  enclosure  are  excited,  then  the  random  acoustic  field 
becomes  directional.  For  this  case,  the  normalized  cross-power  spectral  density  can  be 
approximated  by 


where  kx,  ky,  and  kz  represent  the  wave  numbers  of  the  characteristic  frequencies  of  the 
enclosure  in  the  x,  y,  and  z  directions,  as  defined  in  equation  (2. 1 1 ). 

Studies  described  in  the  remainder  of  this  section  will  consider  acoustic  fields  which  can  be 
described  by  the  simple  trigonometric  functions. 

2.2  ACOUSTIC  COUPLING 

A  description  of  the  acoustic  field  in  a  sonic  fatigue  test  is  one  ol  two  elements 
needed  to  estimate  the  response  of  a  structural  specimen  in  an  acoustic  environment.  The 
other  factor  involves  the  dynamic  properties  of  the  structure  itself.  In  this  section,  the 
important  dynamic  characteristics  of  test  specimens  will  be  discussed,  and  methods  will  be 
shown  for  coupling  the  acoustic  fields  with  the  structures  to  provide  an  estimate  of  the 
structural  response. 

The  structural  system  considered  here  is  assumed  to  be  a  linear  system,  such  that 
the  free  vibrations  of  the  system  can  bo  described  by  a  set  of  homogeneous  second-order 
linear  differential  equations,^ 

IMt  jwj  +  (C|  jwj  +  IK!  jwj  =0  (.2,19) 

If  the  damping  term  is  momentarily  neglected,  the  set  of  equations  can  be  solved  to 
give  the  undamped  natural  frequencies  of  the  structural  system,  and  the  corresponding  mode 
shapes.  Letting  (4>1  represent  the  mode  shapes,  we  cun  introduce  the  transformation  matrix 

Jwj  =  (<J»)  Jr; j  (2.20) 

where  r/  is  the  normal  coordinate.  Substitution  of  equation  (2.201  into  equation  (2.19)  and 
pre-multiplication  by  the  transpose  of  the  modal  matrix  permits  the  equations  of  .notion 
to  be  re-written  in  normal  coordinates  as 

|*|  *  j''j  +  J’lj  c°  U-2D 

Equation  (2.21)  includes  the  effect  of  proportional  damping,  and  displays  the  natural 
frequencies  of  the  structure. wir  Note  that  the  matrices  containing  £  and  to,,  are  diagonal. 
This  property  has  the  effect  of  uncoupling  each  equation  front  the  others.  Liter  in  this 
Section,  a  normalized  forcing  function  known  as  the  participation  factor  will  be  introduc¬ 
ed  that  permits  the  response  at  each  mode  to  be  estimated  by  using  the  uncoupled  equations 
of  motion  in  normal  coordinates. 


Returning  to  equation  (2.21 ),  we  consider  a  system  whose  normal  mode  frequencies  and 
damping  properties  are  known.  This  equation  describes  n  single-degree-of-freedom  systems 
whose  responses  combine  linearly  to  describe  the  behavior  of  the  total  system. 

The  response  of  the  system  at  each  no;-mal  mode  frequency  can  be  discussed  in  terms 
of  the  angular  frequency  and  the  damping.  The  angular  frequency  .  the  n^1  normal  mode 
is  denoted  by  con.  At  that  frequency,  the  effect  of  damping  in  the  system  is  described  by 
fn,  where  £  is  the  non-dimensional  damping  ratio.  If  a  smgle-degree-of-freedom  system 
characterized  by  o>n  and  £n  is  excited  by  a  unit  impulse  function  5(t),  the  response  of  the 
system  is  commonly  designated  by  h(t),  the  unit  impulse  response  function. 

From  h(t),  another  widely  used  property  of  dynamic  systems,  the  complex  frequency 
-esponse  H(a>),  can  be  defined  as 


H(«)  =  f  h(t)e'Jwtdt 

‘'no 


(2.22) 


The  form  of  the  integral  in  equation  (2.22)  indicates  the  frequency  response  function  is  the 
Fourier  transform  of  the  impulse  response  function,  li(t).  The  complex  frequency  response 
function  for  a  single  degree  of  freedom  system  in  non-dimensional  form  is 


H(<o)  = 


(2.23) 


For  typical  aerospace  structures,  equation  (2.23)  describes  the  system  at  each  of  the 
n  normal  mode  frequencies  of  the  distributed  structure. 

Another  property  of  (he  vibrating  structural  specimen  is  the  dynamic  mode  shape 
assumed  by  the  structure  at  a  norma)  mode  frequency.  For  the  present,  pane!  mode  shapes 
will  be  represented  by  characteristic  beam  functions.  This  simplification  will  clarify  the 
acoustic  fiold/struetural  response  coupling  that  is  of  primary  interest  here;  For  example, 
consider  a  panel  with  simply  supported  edges,  and  length  6.  Hie  normalized  mode  shape 
along  the  x -dimension  of  the  panel,  as  x  varies  from  0  to  C.  may  then  be  represented  by 


4>x  =  sin  ~ 


where 


(2.24) 


(n  *  1 .  2,  3, - .) 


(2.25) 


The  term  A,  represents  the  wavelength  of  the  structural  response  mode.  Note  that  for  nK! , 
the  mode  shape  function  is  a  simple  half-sine  wave,  with  maximum  deflection  at  the 
center  of  the  panel.  For  t\-2,  a  full  sine-wave  mode  shape  results, 

Tire  expressions  for  pressure  distribution  developed  in  Section  2.1  can  now  be  com¬ 
bined  with  the  mode  shape  functions  to  yield  a  parameter  known  as  the  participation  factor, 
I\  defined  for  one  dimensional  continuous  system  by 

i  r! 

rx«  -  /  PuM>*  (x)dx  (2,26) 
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The  two-dimensional  case,  which  includes  panels,  is  a  logical  extension  of  equation  (2.26). 
Since  both  the  pressure  distribution  and  the  mode  shape  are  separable  functions  of  x  and  y, 
the  participation  '"actor  for  panels  can  likewise  be  represented  by  a  product  of  the  participa¬ 
tion  fac  tors  in  the  x  and  y  directions. 


r  =  r  r 

x,  y  x  y 


(2.27) 


The  participation  factor  is  a  measure  of  the  extent  to  which  each  mode  shape  d* 
participates  in  the  response  of  the  structure  to  a  given  pressure  distribution,  p(x).  It  is  the 
mean  value  of  the  product  of  the  structural  mode  shape  and  the  pressure  distribution  along 
the  structure.  When  the  structural  mode  shape  and  the  pressure  distribution  have  the  same 
phase  relationship,  strong  coupling  exists  between  the  acoustic  Held  and  the  structure,  and 
significant  structural  response  will  result.  To  illustrate  the  val  .e  of  the  participation  factor 
in  estimating  the  structural  response,  equation  (2.21 )  can  be  rewritten  with  the  excitation 
expressed  in  terms  of  the  participation  factor  for  the  r^‘  mode  as 


+u>r~  'ir^  !>,n 


(2.28) 


The  solution  to  equation  (2. 28)  can  be  st  iiyd  in  terms  of  the  previously-developed 
complex  frequency  response  function  H(te).  f  or  a  simple  harmonic  excitation  at  a  fiequency 
uo,,  the  normalized  response  in  the  r!*'  mode  is  given  by 


b(u)0)  * 


p  r 

*  o  '  r 

T* .  fl(to)  F(w  ) 

Mr 


<2.2‘n 


This  result,  which  is  derived  in  Appendix  A,  is  readily  extended  to  any  periodic  function 
that  can  be  separated  into  its  harmonic  components  through  a  Fourier  scries  analysis  At 
each  of  the  in  harmonics  of  the  fundamental  frequency,  the  solution  becomes 


(>  r 
ro  '  r 

itfutu>0) c  H  iww0i  l- 


(150) 


The  response  of  the  structure  is  found  .by  making  the  transformation  from  normal 
coordinates  back  to  the  original  coordinate  system. 

|w(ui0)j  *  Hd  jr?tw0jj  t2.5l) 

For  random  forcing  .met ions,  the  excitation  and  response  s?f  the  system  will  he  expressed 
in  terms  of  power  spectral  density  functions  The  excitation  can  be  represented  by  j  Sp) .  a 
matrix  containing  the  auto-  and  cross- power  spectral  densities  of  the  excitation  in  normal 
coordinates.  Then,  with  the  complex  response  frequence  Utu?)  as  defined  earlier,  the  re¬ 
sponse  is 

l  t  t-  i 

(2  32) 

where  (Sf,j  is  the  matrix  of  power  spectral  density  functions  for  the  normalized  responses 


[h*(u?)1  is 

fw> 

\  <1 

P  t  N 

In  terms  of  the  participation  factors,  equation  (23 7 i  rs  rewritten  as 


IS  I  = 


H#(uf)  (F.Fjl  H(to)l  S>v 

nJ  '  sj  nn 


(2.33) 
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where  the  P  terms  are  the  participation  factors  for  the  i^  and  modes.  Spjpj  is  the  auto¬ 
power  spectral  density  of  the  excitation  pressure  at  some  arbitrary  point  in  the  .  oustie 
field. 

For  structural  panels,  the  dissipation  is  usually  low,  and  the  equations  relating  the 
response  to  the  excitation  can  be  simplified,  as  discussed  in  Appendix  A.  In  equation 
(2.33),  the  off-diagonal  terms  may  be  neglected,  yielding 


V  Vi 


(2.34) 


This  result  shows  that  the  response  in  normal  coordinates  is  reduced  to  a  function  of 
the  participation  factors  and  the  frequency  response  functions.  The  response  may  be  stated 
in  terms  of  the  original  coordinate  system  by  making  the  transformation  from  normal 
coordinates, 

iS^,!  =  ($1  tS?ji  l«M  1  (2.35) 

Continuing  with  the  participation  factor  concept,  we  can  now  show  how  the  pressure 
distributions  and  mode  shapes  are  related  to  produce  significant  excitation  on  a  test  speci¬ 
men.  As  an  example,  the  pressure  distribution  in  ;t  progressive  wavv  action  with  perfectly 
aneehoie  termination  was  given  in  equation  (2.2).  Neglecting  she  time- varying  component 
and  rewriting  equation  (2.2),  we  can  express  the  spatial  distribution  of  the  pressure  by 


pt  x  I  •*  cos 

*  t 


+ 1  sin  -r 


(2-Jtii 


liquation  (2.3b)  is  for  a  system  with  geometry  as  shown  in  Figure  1 .  The  addition  of  a 
constant  to  the  distance  satiable  x  in  tins  equation  would  cause  the  pressure  distribution 
to  be  currrxpondtngl  shifted  skmg  the  \  avis,  liquation  <  2  5<> )  is  derived  in  Appendix  A 

From  equation  (2,24).  we  can  now  select  a  mode  drape  which  we  wish  to  examine 
for  its  coupling  with  (he  above  pressure  distribution.  If  we  let  n*M.  the  participation  factor 
for  this  combination  of  prexyure  dtstribution  and  mode  shape  is  evaluated  as 


I  /*V  2 ex  2*x  „ 

'  \  /$  I  cos  i“-““  Mn  d  x 

I  L  X*  K 


t2.3?a> 


un  z.H.1  sin 


■cm} 


Tile  vubwnpU  R  and  1  arc  used  to  denote  real  and  imaginary  components,  to  account 
for  the  imaginary  operator  in  the  pi.\)  expression  Tin-  components  are  comhmed  to  give 
the  absolute  value  of  the  magnitude  of  the  p  srticipatton  factor  for  the  progressive  wave 
section  with  perfectly  ancehose  termination. 


r-tr 


■>  ,  ■>  i 

♦  r,-i 


(2.3S) 
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Participation  factors  for  other  terminations,  where  the  pressure  distribution  is  represented 
by  a  single  real  trigonometric  function,  are  found  in  a  similar  manner.  Only  a  real  value  for 
T  will  result,  as  in  equation  ( 2.37a )  For  a  progressive  wave  section  with  a  rigid  termination, 
the  pressure  distribution  takes  (he  form 


p(x) 


Po  ^os 


X, 


(2.34) 


where  the  exact  location  q"  the  distribution  within  an  enclosure  can  be  shifted  by  adding  a 
constant  to  the  x  variable.  The  participation  factor  for  n-1  mode  now  becomes 


r=  i/t  /  cos  ~ 


irx 


:n\ 


sin 


Xr 


dx 


(2.40) 


It  is  interesting  to  compare  the  progressive  wave  test  section  with  a  normal  incidence 
test  section.  In  normal  incidence  testing,  the  acoustic  waves ;  ravel  in  a  direction  perpendicular 
to  the  test  specimen  surface,  and  the  excitation  pressure  're  essentially  uniform  across  the 
specimens  At  low  values  of  wavelength  ratio,  where  the  structural  response  wavelength  is 
much  less  than  the  acoustic  excitation  wavelength,  the  pressures  are  also  uniform  along  the 
specimen.  Under  the  condition  of  low  wavelength  ratio,  a  consideration  of  the  specimen 
response  would  show  no  difference  between  the  progressive  wave  test  section  and  the  normal 
incidence  arrangement. 

The  participation  he  I  of  provides  some  insight  into  the  coupling  between  »he  potentially 
large  number  of  combinations  ot  structural  modes  and  excitation  frcqucitcn ».  In  the  above 
expressions  for  t\  the  term  X  represents  !•«.-  wavelength  of  the  structural  response,  as  re¬ 
lated  to  the  geometry  "f  the  structure  according  to  equation  *.2  2>  ).  Xp,  on  the  other  hand, 
is  the  wavelength  associate*'!  with  the  excitation  pressure  frequency ,  as  determined  from  the 
usual  relationship  between  wavelength,  velocity  ot  propagation,  and  frequency. 

.  \  «  J  •! 

4  »  ■"  V  i 

V 

Wc  Can  mm  define  a  new  parameter,  wive.,  -gth  r.ttn  .  X(.<\,,.  such  tlut 

X.  .Mu  V  3  , 


t  2  41  » 


wavelength  ratio  -<  ~L 


e  b 


oc 


losing  w-yvclcngth  ratio  as  an  independent  unable,  par (ictpat toil  (actors  can  be  evaluated 
fur  vuiKun  structure  modes  that  is,  for  a  given  value  of  it,  (hr  participation  factor  F  can 
be  evaluated  as  the  wavelength  ratio  ts  allowed  to  vary  f  igure  3  is  a  family  of  curves  for 
the  participation  factor,  as  given  to  equations  (2  3  ?i,  for  a  progressive  wave  section  with 
perfectly  anechorv  termination  The  participation  factors  are  presented  as  RMS  quantities 
to  facilitate  addition  of  comp  mentv  1  he  m- phase  component  of  participation  factor  has 
a  maximum  due  tor  the  tod  structural  mole  ' n-- i )  at  a  wavelength  ratio  of  zero  Tips 
represents  the  cave  where  the  wavelength  of  the  excitation  pressure  is  much  larger  than 
the  structural  response  wavelength,  so  that  the  structure  appears  to  ty  excited  by  a  pres¬ 
sure  distribution  that  is  umfot  a  across  the  structure  The  higher  order  modes  have  a 
maximum  value  of  part *eipa turn  factor  near  a  wavelength  ratio  of  one.  The  quadrat- v 
component  4  parti  t pat  ton  factor  always  has  a  max. mum  value  at  a  wavelength  ratio  oj 
one.  ahd  a  value  of  Zero  at  a  wavelength,  ratio  of  mo.  it  is  significant  that  the  participa¬ 
tion  factor  becomes  very  small  for  wavelength  ratios  greater  than  i\  >  At  larger  wavelength 
ratios,  acoustic  energy  becomes  very  inefficient  as  a  me  hamsm  for  producing  s»reas  in  a 
structural  specimen  In  Figure  4.  the  total  participation  factor  has  been  plotted,  according 
to  equation  (2.38). 
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Parallel  to  Wave  Propagation 


Tiie  above  discussions  have  considered  only  those  structural  modes  in  the  x  direction, 
the  direction  of  the  pressure  distribution.  In  a  progressive  wave  test  section,  the  pressure  is 
essentially  uniform  in  the  y  direction  (perpendicular  to  the  wave  propagation).  Therefore, 
if  we  consider  the  structural  modes  in  the  y  direction,  the  participation  factor  will  be  zero 
for  all  the  odd-number  modes.  This  result  is  plotted  in  Figure  4,  where  the  participation 
factor  has  been  calculated  from  an  expression  similar  to  equation  (2.26)  except  that  the 
variable  has  been  replaced  by  y.  With  uniform  pressure  distribution  in  the  y  direction,  p(y) 
becomes  a  constant,  and  the  points  shown  in  Figure  4  result  from  integrating  the  simple 
trigonometric  mode  shapes.  As  in  the  case  for  participation  factors  in  the  x  direction,  it  is 
apparent  that  the  acoustic-structural  coupling  is  very  poor  at  all  wavelength  ratios  except 
ratios  near  unity. 

As  examples  of  how  the  wavelength  ratio  and  participation  factor  concepts  can  be 
used  to  estimate  structural  response,  two  types  of  structural  configurations  -  beams  and 
plates  -  will  be  considered  here.  This  discussion  will  assume  that  the  excitation  frequency 
coincides  with  the  structural  resonance  frequency,  again  recognizing  that  the  acceptance 
of  energy  by  structural  members  occurs  principally  at  resonances.  For  rectangular  cross- 
section  beams,  the  bending  frequencies  are  determined  from 


h 

fr  =  ^ 


(2.43) 


where  h  and  £  are  the  thickness  and  length  of  the  beam,  respectively,  and  aj  is  a  constant 
determined  by  the  end  conditions  of  the  beam  and  the  order  of  the  mode  being  considered. 
For  many  commonly  used  engineering  materials,  the  E/p  ratio  does  not  change  widely,  and 
a  value  of  1 .02  x  10§  (inches)  was  used  here.  The  expression  for  the  resonant  frequency  of 
the  structure  may  then  be  written  as 
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(2.44) 


to  show  the  desired  dependence  on  thickness  and  length  only.  3f  we  now  constrain  the 
frequency  of  the  excitation  pressure,  fp,  to  be  equal  to  the  structural  resonance,  and  solve 
for  the  wavelength  of  the  pressure  from  equation  (2.41 ),  we  obtain 
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Finally,  using  the  expression  for  structural  response  wavelength  as  given  in  equation  (2.25) 
for  the  first  resonant  mode  of  the  beam  (n— 1 ),  the  wavelength  ratio  is  determined  as 
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where  the  new  constant,  a-',  includes  the  effect  of  boundary  conditions,  mode  number, 
speed  of  sound  and  physical  properties  of  the  beam  material.  The  values  of  wavelength  ratio 
for  three  end  conditions  for  beams,  as  a  function  of  h/fc,  are  plotted  in  Figure  5.  These 
curves,  then,  become  the  first  step  in  estimating  the  response  of  a  beam  to  an  acoustic  ex¬ 
citation.  From  the  h/£  ratio  and  end  conditions  of  the  beam,  the  wavelength  ratio  is 
determined  directly  from  the  curve.  The  wavelength  ratio  is  used  to  enter  a  curve  of 
participation  factor  for  the  type  of  acoustic  field  being  considered,  as  for  example.  Figure 
3.  The  relative  magnitude  of  the  participation  factor  gives  an  indication  of  the  effectiveness 
of  a  given  test  method  for  inducing  response  in  a  structure. 

Square  panel  structures  have  been  treated  in  a  manner  similar  to  the  beam  analysis, 
and  the  results  are  plotted  in  Figure  6.  Three  edge  conditions  are  considered  for  panels: 
fixed,  simply  supported,  and  free.  As  in  the  case  of  the  beams,  results  are  shown  only  for 
the  first  mode  resonances.  For  panels,  however,  it  is  important  to  show  how  the  wavelength 
ratio  changes  with  higher  order  modes.  The  curves  in  Figures  7, 8,  and  9  are  plots  of  nor¬ 
malized  wavelength  ratios  for  square  plates.  Values  of  mode-number  in  the  x-direction  are 
used  as  the  abscissa,  while  mode-numbers  in  the  y  direction  are  shown  as  constant  lines  on 
the  plot.  The  structural  response  wavlength  in  the  x  direction  (direction  of  acoustic  propa¬ 
gation)  has  been  used  in  the  numerator  of  the  wavelength  ratio  expression.  Each  wavelength 
ratio  is  normalized  relative  to  the  wavelength  ratio  for  the  first  mode, 
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Examination  of  Figures  7,  8,  and  9,  along  with  a  plot  of  the  participation  factor  for 
the  type  of  acoustic  enclosure  used  on  a  particular  test,  provides  some  insight  into  the  coup¬ 
ling  that  might  be  expected  between  a  specimen  and  an  acoustic  field. 

2.3  TEMPERATURE  EFFECTS  ON  ACOUSTIC  FIELDS 

Experiments  were  conducted  to  determine  the  effect,  of  thermal  energy  on  the 
acoustic  field.  The  experimental  apparatus  (described  in  more  detail  in  Section  4,  Figure 
26)  consisted  of  a  rectangular  cross-section  test  enclosure  with  a  bank  of  16  quartz  lamps, 
each  rated  at  1600  watts,  heating  a  1 5-inch-square  panel.  Tests  were  conducted  at  panel 
temperatures  ranging  from  1 50°F  to  475°F.  Measurements  were  made  to  determine  the 
power  required  to  hold  the  panel  center  temperature  at  various  levels,  both  in  still  air  and 
with  an  air  flow  over  the  heated  panel.  The  results  are  plotted  in  Figure  10,  where  the  air 
mass  flow  is  shown  per  unit  cross-section  area  of  the  enclosure.  The  increase  in  electrical 
power  generates  heat  tiat  is  lost  through  several  mechanisms  -  the  test  panel,  the  surround¬ 
ing  support  structure,  the  lamp-bank  reflector,  and  the  air  in  the  test  section.  If  it  is  assumed 
that  all  of  the  increase  in  electrical  power  is  converted  into  an  increase  in  the  gas  stagnation 
temperature,  then  the  increase  in  gas  temperature  can  be  computed  from  the  power  in¬ 
crease,  mass  flow;  and  specific  heat  from' 
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FIGURE  6 

FIRST  MODE  WAVE  LENGTH  RATIO 
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The  increase  in  power  at  a  panel  temperature  of  475°  was  approximately  2600  watts. 
This  power  increase,  combined  with  a  mass  flow  for  the  gas  of  2.5  pounds  per  second, 
corresponds  to  an  increase  of  4°F  in  the  gas  temperature.  This  small  change  in  gas  temper¬ 
ature  was  verified  by  the  experimental  work,  where  no  increase  in  gas  temperature  was 
observed.  Since  the  increase  in  gas  temperature  can  be  considered  negligible,  corresponding 
changes  in  gas  properties  (speed  of  sound  and  density)  can  also  be  neglected.  This  conclusion 
is  further  verified  by  comparing  the  transfer  functions  for  the  pressure  data  in  Figures  1 1  and 
1 2.  These  curves  represent  transfer  functions  between  two  microphones  at  two  temperature 
conditions.  One  microphone  was  located  in  the  thermal  field  near  the  center  of  the  panel, 
and  its  output  is  normalized  with  respect  to  a  reference  microphone  located  at  tire  input  to 
tho  test  section,  but  outside  of  the  thermal  field. 

The  effect  of  a  significant  increase  in  gas  temperature  on  the  acoustic  fields  would  be 
to  increase  c,  the  velocity  of  sound,  and  decrease  p0,  the  equilibrium  density  of  the  gas. 

These  changes  would  invalidate  the  basic  vvuw  equations,  which  are  derived  on  the  basis  of 
the  equilibrium  density  and  velocity  of  propagation  being  uniform  throughout  the  medium. 
Introduction  of  spatial  variations  in  c  and  pQ  would  requite  the  derivation  of  a  new  set  of 
wave  equations.  As  a  practical  matter,  however,  the  air  flows  that  are  normally  associated 
with  conventional  acoustic  noise  sources,  such  as  air-stream  modulators,  are  sufficient  to 
keep  the  air  at  a  near-ambient  temperature.  Consequently,  the  effect  of  air  temperature 
changes  on  the  acoustic  fields  can  be  neglected  in  elevated  temperature  sonic  fatigue  teso. 

7  4  TEST  TECHNIQUES  -  ACOUSTIC  ENVIRONMENT  SIMULATION 

Numerous  test  techniques  are  available  for  conducting  sonic  fatigue  tests  on  panel  struc* 
tores,  in  Sections  24  through  2.3,  the  acoustic  fields  for  some  fatigue  tests  were  described, 
end  the  coupling  of  the  acoustic  fields  with  structural  specimens  was  discussed.  In  this  section, 
earlier  results  will  be  used  in  a  numerical  example  to  illustrate  application  of  those  results 
to  an  experimental  situation, 

tteforc' proceed)^  with  tire  example,  it  diottld  Ire  noted  that  the  ultimate  objective  of 
a  sonit.  fatigue  test  is  to  induce  a  certain  dynamic  stress  in  the  test  specimen.  In  meeting 
this  ybjnetive,  it  might  not  Ire  necessary  to  simulate,  in  great  detail,  the  acUigUn-service 
acoustic  environment .  A  study  of  the  dynamic  properties  of  the  structure  tlunfid  reveal 
•thS  ma$t4ike*y  service  fmluie  mode  of  the  specimen,  mad  the  acoustic  field  should  tor 
tailored  to  ityhicvsr'this  failure  mode  in  the  labors  tory. 

Tire  lest  jnmcl  yH«c»fd  for  tire  present  example  is  a  square  with 

dimensions  df  J  2  %  1 2.x  i>.G32  inches. 'The  panel  is  considered  to  ho  mounted,  in  a  manner 
tha*  provides  simply  supported Tdges;.  -the  first  panel  cfisracicmtk  of  imereti  is  fire  tideferess- 
•Fd-Umgth  reito,  which  %  0.0326  Exttapcdatlon  of  tire  cum  fi  gmes  a; value  of 

0,073  for  the  first  mods  Wavelength  ratio,  Xr/XJV  Next,  we  Bgtire  0,  Which  shows  ; 
the  normalized  wavelength  ratio  for  tire  higher  o;der  modes.  Ore  values  of  wavelength  ratio 
in  Figure  9  are  rretmateed  with  respect  to, tire  wavelength  retie  st  ilk  tijrti; mq&Nso.  the . 
values  tekef  r*D*n  Figure  9  most  be  multiplied  byU-673  tt>  give-tire  actual  wayHeogtii -retiti;- . 
6»  this  specimen.  Table  f  presents  both  the  normalized  wavelength  ratio.  and  tire 
actual  value,  for  tresde  number  variations  -from  for  both  mjy  . 

direction!.  ..  —  . 
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It  wiii  now  be  assumed  that  the  panel  ill  be  installed  in  the  wall  of  a  progressive  wave 
test  section  vith  perfectly  anechoic  termination.  The  total  participation  factor  will  be  the 
product  of  me  participation  factors  in  the  x  and  y  directions,  rxrv,  according  to  equation 
(2.27).  The  x  and  y  participation  factors  can  be  taken  from  the  curves  in  Figure  4.  The 
values  of  participation  factor  are  tabulated  m  Table  2  for  the  same  variation  in  mode  numbers. 

An  inspection  of  the  values  of  FxFy  in  Table  2  shows  that  the  first  mode  (n=m=i )  is 
the  most  strongly  excited  mode.  The  excitation  of  the  next  strongest  mode  (n-1 ,  m=3)  is 
only  30%  of  the  first  mode,  and  the  other  modes  are  even  less  responsive. 

Use  of  the  wavelength  ratio  parameter  also  permits  us  to  compare  the  progiessive  wave 
test  section  section  with  a  normal  incidence  test  arranges .  \i.  We  reball  that  the  first  mode 
ratio  was  0.073  for  the  progressive  wave  setup,  giving  a  participation  factor  0.20.  For  normal 
incidence  testing,  the  pressure  across  the  specimen  would  be  uniform,  corresponding  to  a 
Xp  that  approaches  infinity,  and  a  corresponding  wavelength  ratio  of  zero.  The  zero  wave¬ 
length  ratio  for  normal  incidence  testing  gives  essentially  the  same  participation  factor  as 
that  for  the  progressive  wave  section. 

A  similar  argument  can  be  made  for  testing  this  panel  in  a  large  reverberation  chamber. 

In  this  case,  we  must  first  determine  if  the  acoustic  field  appears  te  be  directional  or  diffuse 
in  the  frequency  range  of  the  1,1  mode  of  the  pans!  specimen.  The  first  mode  frequency  - 
"an  be  computed  from  the  first  mode  wavelength  ratio  and  the  panel  length 


r . , * c  ;V'¥ 

M  to 


(2.49)' 


Equation  (2.49)  yields  a  first  inode  frequency  of  40.5  Hz  for  the  panel.  If  we  now  assume 
that  the  panel  is  lightiy  damped,  and  let  f  =  0.01  for  this  mode,  we  can  enter  Figure  2  and  ' 
see  that  a  room  volume  of  approximately  600,000  fP  would  be  required  in' order  to  obtain  a 
truly  diffuse  acoustic  field  at  the  first  mode  resonance.  From  equation  (2.16),  we  see  that 
a  room  volume  of  80,000  ftr  would  be  required  to  have  even  a  single  mode  excited  within 
the  bandwidth  of  the  structural  resonance.  The  acoustic  field  for  such  a  room  would  be 
considered  purely  directional  at  the  first  ..mode  resonance,  and.  the  same  participation  factor 
computed  for  the  progressive  wave  section  would  apply. 

This  example  indicates  that  the  choice  of  grazing  or  normal  incidence  field  may  not 
have  much  significance  from  a  consideration  of  the  response  induced  in  the  test  specimen. 
Enclosures  have  perfectly  anechoic  terminations  will  result  in  both  sine  and  cosine  terms,  .... 
in  the  pressure  distribution,  and  will  consequently  have  a  higher  probability  of  exciting 
structural  modes. 


SECTION  3 

THERMAL  ENVIRONMENT  SIMULATION 

laboratory  simulation  of  the  thermal  environment  is  an  important  aspect  cr  onic 
fatigue  testing  at  elevated  temperatures.  The  usual  problems  of  temperature  sinuu...ion  ate 
compounded  by  the  requirement  that  the  thermal  source  apparatus  be  capable  of  surviving 
in  the  high  intensity  acoustic  sound  field.  This  section  describes^  procedure  for  calculating 
the  temperature  distribution  in  a  heated  panel.  A  discussion  of  methods  for  implementing 
heat  sources  for  sonic  fatigue  testing  is  then  presented.  -  •'  "V 

3.1  THERMAL  FIELD  STUDIES 

Prediction  of  the  temperature  distribution  in  a  heated  sonic  fatigue  specimen  is  based 
primarily  on  the  intensity  of  the  heat  source  and  the  thermodynamic  properties  ofthe  test 
specimen.  A  secondary  consideration  is  the  effect  of  air  flow  associated  with  the  acoustic 
noise  sources  in  a  given  experimental  arrangement.  The  latter  effect  will  be  neglected  while 
the  procedure  for  predicting  temperatures  from  a  radiant  heat  source  is  discussed  .'The  analy 
sis  is  described  for  a  simple  panel,  but  can  be  modified  to  account  for  such  conditions  as  " 
nonuniform  heat  flux,  variable  edge  temperatures,  and  configurations  that  represent  more 
complex  aerospace  panels. 

The  analytical  model  used  is  shown  in  Figure  13.  The  panel  was  considered  to  be 
symmetrical  about  a  center-line  parallel  to  the  direction  of  air  flow,  with  the  i  and  j 
directions  as  shown.  Square  elements,  with  sides  equal  to  Ax,  were  chosen  for  simplicity 
in  calculation.  The  intersections  of  the  grid  lines  were  the  points  at  which  the  temperature 
was  to  be  computed. 

At  any  intersection  (i,  j)  a  heat  balance  may  be  written.  The  algebraic  sum  of  the  con¬ 
duction  from  each  of  the  four  neighboring  intersections,  plus  the  heat  generated  (absorbed) 
must  be  equal  to  zero  in  the  steady  state.  The  heat  generated  from  absorbed  radiation  is 
expressed  as  -  .  . 

%en  =  F  <Ax>2  (3.1) 

where  F  is  the  absorbed  radiant  heat  flux  in  BTU  per  second  per  unit  area.  Actual  values  of 
F  that  may  be  used  in  a  given  experimental  situation  will  be  discussed  in  Section  3,2. 

The  conduction  from  a  neighboring  intersection  may  be  written 
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qCond.  “  k  ^T-l,  j  ^j)  (3.3) 

where  qconcj  is  the  heat  conducted  from  the  (i — 1 ,  j)  intersection  to  the  (i,  j)  intersection. 
Similar  expressions  are  written  for  the  three  remaining  intersections,  and  these  are  com¬ 
bined  with  equation  (3.1 )  to  yield  the  heat  balance  equation 

k  Az j  +Tjt j+j  ♦TWiJ  +  Ty_, -4T|j!  *«W-0  (3.4) 

Solving  for  Tj  j  yields 
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FIGURE  13 

ANALYTICAL  MODEL  FOR  TEMPERATURE  PREDICTION 


ELEMENT  OF  THERMAL  MODEL 


Equation  (3.5)  must  be  evalvted  for  each  intersection  point  on  the  grid.  An  iterative 
technique  is  used,  in  which  an  initial  temperature  distribution  is  assumed.  A  “new”  tem¬ 
perature  is  calculated  for  each  point  using  the  “old”  temperature  for  the  neighboring  - 
points  in  the  calculation.  Each  “new”  temperature  then  becomes  the  “old”  temperature 
for  the  next  series  of  calculations,  and  this  procedure  is  repeated  as  the  difference  between 
the  “new”  and  “old”  temperature  becomes  vanishingly  small.  When  this  difference  is  less 
than  some  specified  value,  the  calculation  is  considered  complete. 

The  analysis  must  also  account  for  the  radiation  into  and  away  from  the  test  specimen. 

A  given  test  setup  will  include  an  array  of  quartz  lamp  radiant  heaters  backed  up  by  reflectors 
on  as  many  sides  as  possible,  consistent  with  the  air  flow  requirements  of  the  acoustic  noise 
sources.  To  enhance  the  radiant  heat  transfer  efficiency  of  the  setup,  the  lamp  array  and 
reflectors  would  be  placed  close  to  the  test  panel. 

The  thermal  model  for  the  radiative  heat  transfer  problem  considers  the  lamp 
reflectors  to  be  cooled,  high  reflectivity,  diffuse  reflecting  surfaces.  The  specimen  panel 
“sees”  above  itself  a  diffuse  reflection  of  itself;  i.e.;  another  panel  with  approximately 
"  the  same  emissivity  and  surface  temperature.  The  equation  expressing  the  net  radiant  power 
transferred  out  of  the  (i,  j)  intersection  in  such  a  configuration  is  written 

floutfi,  j)  ~  0  ^e^a  (l  i,  j  “  )  @.6) 

The  emissivity  factor,  Fe,  for  the  case  of  a  small  area  Aj  of  emissivity  e  j ,  surrounded 
by  a  larger  area  A2  of  emissivity  ej  can  be  expressed  as 

F  1 

e~1/ei  +  ^i/i  -  A  (3-7) 

A2  \f2  / 

When  A2  is  much  greater  then  Aj ,  which  is  the  case  for  the  assumed  model,  it  is  noted 
that 
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Similarly,  for  this  geometry,  the  view  factor  Fa  (ratio  of  the  solid  angle  subtended  by 
A2  from  A  j  to  the  solid  angle  of  hemispherical  space)  approaches  1 .0  as  A  j  /An  approaches 
zero.  Tr  is  the  effective  reradiative  temperature  produced  by  the  diffuse  reflector  reflecting 
the  specimen  panel  radiation.  The  reradiative' temperature  is  considered  to  be  spatially  uni¬ 
form,  and  equal  to  a  fourth  power  average  of  the  specimen  plate  temperature,  as  follows. 
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where  M  and  N'are  the  number  of  interior  grid  lines  an  the  panel  in  tile  i  and  j  directions, 
respectively. 


Equating  the  heat  balance  equation  (equation  (3.4))  to  the  radiant  exchange  (equa¬ 
tion  (3.6))  yields 
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(3,10) 
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which  is  the  basic  equation  to  be  solved  for  Ty.  The  presence  of  the  fourth  power  terms 
complicates  the  solution,  and  so  a  linearization  scheme  has  been  applied  to  the  equation. 
At  the  same  time,  an  algorithm  is  used  which  prevents  the  iterative  solution  from  diverg¬ 
ing.  In  linearizing  the  fourth  order  term,  the  “old”  temperatures  are  assigned  the  value 

while  the  “new”  temperature  remains  Tj  j.  Using  this  nomen- 
as 


j  Aj  j+{, . . .  etc.,  whil 
clature,  [Tjjj  is  written  as 


4  3  4 

TU  =  4AiJTU-3AiJ 

Then,  using  this  result  in  equation  (3.1 1 ),  the  basic  equation  becomes 
a  .  »  A  a  .  a  „  x  ,  f(Ax)2 

VlJ+AiJ+l+Ai+lJ  +  AiJ-l  -4Ti,j+-j^— 

kAz  kAz 

As  a  further  simplification,  two  new  constants  are  defined, 

F(Ax)2 


(3.12) 


(3.13) 


K»- 


kAz 


and 


K.  * 


aej  (Ax)^ 
kAz 


Finally,  solving  for  Tj  j  yields 


T.  .  »  - 

\  J  4 


1  (Ai-i,  j  +  Ajy+i  -+ Aj+i  j  +  +  K  +  3Kj  Aj>  j  +KjTr 


(3.14) 


(3,15) 


(3.16) 


i+KjAu 


This  equation,  then,  is  solved  in  exactly  the  same  manner  as  equation  (3.5)  for  the 
conduction  case  only.  A  computer  program  for  the  temperature  distribution,  as  defin¬ 
ed  in  equation  (3.16),  is  presented  in  Appendix  1). 


3.2  DESIGN  OF  HEAT  SOURCES 


Quartz  lamp  banks  and  their  associated  reflectors  were  used  as  heat  sources  in  this 
study  effort.  While  quartz  lamps  have  been  widely  used  in  standard  laboratory  conditions 
fov  a  number  of  years,  the  introduction  of  these  lamp  banks  into  a  high  intensity  acoustic 
environment  has  brought  new  problems.  In  the  following  paragraphs,  the  overall  considera¬ 
tions  for  a  lamp  bank  design  will  first  be  considered.  Then,  some  details  of  the  lamp  installs 
tion  will  be  discussed. 


3.2.1  Reflector  Design  Considerations 

The  purpose  of  a  quartz  lamp  bank/reflector  assembly  is  to  convert  electrical  energy 
into  thermal  energy  that  is  absorbed  on  a  test  specimen.  The  efficiency  with  which  this 
objective  is  reached  is  a  measure  of  the  success  of  a  given  design.  A  number  of  high  quality 
quartz  lamp  banks  are  available  commercially.  The  specified  performance  of  these  lamp 
banks  can  be  used  in  equation  (3.14)  as  part  of  die  heater  design  for  sonic  fatugue  test¬ 
ing.  The  basic  parameter  for  this  purpose  is  the  radiated  heat  flux  from  the  lamp  bank. 

For  a  representative  lamp  bank,  75  watts  per  inch  2  (approximately  10  BTU  per  ft  2-second) 
can  be  conservatively  obtained  at  the  face  of  the  lamp  bank.  This  output  is  fdr  each  lamp 
in  the  bank,  based  on  a  lamp  spacing  of  1  inch  between  centers,  and  electrical  power  of 
100  watts  per  inch  applied  to  the  lamps.  This  example  represents  an  electrical-to-thermal 
power  conversion  efficiency  of  75%.  Higher  flux  densities  may  be  obtained  by  a  com¬ 
bination  of  closer  lamp  spacing  and  higher  power  applied  to  the  lamps,  with  a  correspond¬ 
ing  decrease  in  the  useful  life  and  reliability  of  the  lamp  bank. 

TV  actual  value  of  heat  flux  absorbed  by  the  s^eimen  depends  on  a  number  of 
factors,  such  as  the  distance  of  the  lamp  bank  from  the  specimen,  the  absorptivity  of  the 
specimen,  and  losses  from  the  specimen.  Therefore,  the  value  for  F  used  in  equation 
(3.14)  should  be  lower  than  the  nominal  heat  flux  that  a  lamp  bank  can  generate.  A 
bright,  highly  polished  specimen,  for  example,  generally  absorbs  less  radiant  heat  energy 
than  a  dull,  darkened  specimen.  Likewise,  a  panel  with  insulation  on  the  side  away  from 
the  radiant  heater  could  be  raised  to  a  specified  temperature  with  less  applied  heat  than 
a  similar  uninsulated  panel.  These  factors  can  be  controlled  to  a  large  extent  in  sonic 
fatigue  tests. 

As  stated  previously,  the  acoustic  environment  introduces  new  problems  in  the  use 
of  quartz  lamp  banks.  Existing  commercial  equipment  is  designed  to  be  used  in  a  n blent 
conditions.  Therefore,  a  special  quartz  lamp  reflector  was  designed  for  use  in  the  pu  -ent 
study.  The  reflector  is  con  ventional  in  several  respects.  First,  a  sheet  of  0. 1  (Much- thick 
copper  is  used  as  the  basic  reflector.  The  reflector  was  drilled  to  permit  installation  of 
sixteen  T3/CL  qu.i.  tz  heat  lamps  at  a  lamp  spacing  of  1.0  inch  These  are  standard  16* 
inch  long,  1600-wait  quartz  lamps. 

A  unique  feature  of  the  reflector  assembly  is  the  reflective  surface  itself,  Most 
reflectors  are  designed  with  a  specular  finish.  The  reflector  for  this  study,  however,  was 
designed  to  provide  a  diffuse  surface,  while  retaining  a  high  degree  of  reflectivity.  Therefore, 
the  side  of  the  reflector  facing  the  quartz  lamps  was  surfaced  with  a  sheet  of  fine  emery 
paper,  onto  which  a  very  thin  layer  of  gold  had  been  vacuum  deposited.  The  resulting  sur¬ 
face  had  the  desired  diffusion  properties,  combined  with  a  reflectivity  of  approximately  95%. 

Another  departure  from  conventional  reflector  design  was  the  introduction  of  damp¬ 
ing  into  the  reflector  to  reduce  its  response  to  the  acoustic  field.  A  rigid  plate  (1,0-inch 
thick  aluminum)  was  bolted  to  the  back  of  the  copper  reflector  plate;  with  a  constrained 
layer  of  RTV  655  Silicone  Rubber  between  the  reflector  and  the  back-up  plate.  Although 
RTV  655  breaks  down  at  temperatures  above  S00°F,  the  material  was  not  permitted  to 
reach  that  temperature  is  this  installation.  The  relatively  low  temperatures  maintained  in  the 
R  l  V  055  resulted  from  a  combination  of  the  low  thermal  transmission  of  the  special  re¬ 
flective  surface  and  the  cooling  effect  of  the  water  circulated  through  the  copper  coils  on 
the  back  of  the  reflector,  lire  damping  effect  provided  by  the  RTV  655  rediwed  the  response 
of  tile  main  reflector  to  an  acceptable  level  throughout  tire  high  temperature  acoustic  tests. 


The  thermodynamic  performance  of  the  lamp  bank/reflector  assembly  was  evaluated 
in  a  sonic  fatigue  test  installation.  The  total  electrical  power  required  to  heat  the  test  panel 
to  a  given  center  temperature  was  recorded  with  no  air  flow,  and  then  with  an  air  flow 
representative  of  a  sonic  fatigue  test.  Since  commonly  used  acoustic  noise  sources  modulate 
a  stream  of  compressed  air,  some  air  flow  can  be  anticipated  in  elevated  temperature  sonic 
fatigue  tests.  Figure  1 0  shows  the  electrical  power  required  to  maintain  a  given  center 
temperature  on  the  test  panel  as  the  air  flow  is  increased  to  a  typical  value  for  the  Wyle 
WAS-3000  noise  source.  The  curve  is  given  in  units  of  air  mass  per  second  per  unit  cross 
sectional  area  of  the  test  enclosure.  This  data  indicates  that  the  additional  power  required 
to  maintain  specimen  temperature  in  the  presence  of  air  flow  does  not  exceed  the  capacity 
of 0  representative  lamp  bank. 

n.  ..Jdition  to  structural  dynamic  and  thermodynamic  influences  on  reflector  design, 
the  effect  of  the  reflector  on  the  acoustic  field  must  be  considered.  Three  types  of  acoustic 
test  arrangements  will  be  discussed:  progressive  wave,  normal  incidence,  and  large  reverber¬ 
ant  enclosures.  Tire  last  two  -  normal  incidence  and  reverberant  -  can  be  considered  at  the 
same  time,  since  in  both  arrangements  a  lamp  bank  reflector  will  effectively  mask  the  test 
specimen  from  the  acoustic  excitation  field.  If  the  lamp  bank  is  moved  away  from  the 
specimen  to  permit  the  acoustic  field  to  reach  the  specimen,  a  corresponding  reduction  in 
heat  flux  density  is  experienced.  This  method  might  prove  satisfactory  in  cases  where  a 
relatively  low  percent  of  the  available  energy  from  a  lamp  bank  is  needed  to  bring  the 
specimen  to  test  temperature.  Alternatively,  it  might  be  possible  to  heat  the  test  specimen 
from  the  side  away  from  the  acoustic  field.  In  the  ease  of  Haiti  panels  with  no  stiffeners, 
this  technique  should  cause  no  problem,  For  a  panel  with  structural  members  on  its  back 
face,  heating  from  the  back  would  result  in  hot  spots  on  the  structural  members,  while 
the  specimen  panel  itself  would  not  receive  uniform  heating,  the  structural  members  could 
be  insulated  from  the  heating  syst  .-.n  to  reduce  the  hot  spots,  but  this  would  not  alleviate 
the  non-uniform  heating  of  the  specimen  panel. 

For  tite  progressive  wave  test  sections,  the  lamp  bunk  has  far  less  effect  on  the  acoustic 
field.  The  lamp  bunk  can  be  oriented  so  that  only  the  edge  of  the  reflector  presents  Itself 
to  the  acoustic  field.  The  thickness  of  the  panel  is  small  compared  to  the  cross-section  dimen¬ 
sion  of  the  enclosure,  resulting  in  virtually  no  distortion  of  the  acoustic  field.  While  it  is 
possible  to  set  up  standing  waves  between  the  reflector  and  the  test  specimen .  the  frequen¬ 
cies  at  w  hich  standing  waves  would  appear  depend  on  the  spacing  between  the  reflector 
and  the  test  specimen,  A  typical  spacing  would  be  in  the  range  of  three  to  six  inches,  which 
would  support  standing  waves  at  2000  Hz  and  1000  Hz.  respectively.  These  frequencies  are 
above  the  range  in  which  fundamental  structural  resonances  normally  occur.  Under  these 
conditions,  the  effect  of  the  lamp  banks  on  the  acoustic  field  in  a  progressive  wave  test 
section  un  be  neglected. 

3.2.2  Ouariz  Lamp  Installations 

Normal  installations  of  quartz  lamps  in  reflectors  do  not  consider  the  dynamic 
motion  of  the  reflector.  The  high  intensity  acoustic  noise  environment,  however,  causes 
'Sine  motion  of  the  reflector,  even  in  the  highly-damped  reflector  described  above,  A  series 
of  experiments  was  conducted  to  determine  the  ability  of  quartz  lamps  to  withstand  dynam¬ 
ic  environments. 
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In  one  test,  an  unpowered  quartz  lamp  was  instrumented  with  strain  gages  to  detect 
bending  strain  in  the  quartz  envelope  of  the  lamp.  The  lamp  was  then  freely  suspended  in 
front  of  an  acoustic  noise  source,  and  subjected  to  increasing  sound  pressure  levels  up  to 
161  dB.  The  strain  in  the  quartz  envelope  was  recorded  at  several  sound  pressure  levels,  as 
plotted  in  Figure  14.  Extrapolation  of  this  data  indicates  that  the  quartz  lamps  will  sustain 
acoustically-induced  loads  at  sound  pressure  levels  in  excess  of  1 70  dB. 

Tins  conclusion  was  substantiated  by  conducting  vibration  tests  on  a  quartz  lamp.  For 
the  vibration  testing,  the  cuds  of  the  lamp  were  rigidly  bonded  to  a  vibration  fixture,  permit 
ting  excitation  to  be  applied  through  the  ends  of  the  lamp.  The  lamp  was  excited  at  input 
levels  up  to  30  g’s  at  the  first  bending  mode  frequency  of  the  lamp,  with  no  failure  of  the 
envelope.  Failure  of  a  filament  occurred  at  21  g’s,  at  which  level  the  stress  in  the  quartz 
envelope  was  approximately  3600  psi,  which  is  half  of  the  nominal  ultimate  strength  of 
the  quartz  material.  The  vibration  testing  was  considered  to  be  a  far  more  severe  test  of 
the  lamp  than  the  acoustic  testing,  since  no  filament  failures  were  experienced  during 
acoustic  testing.  In  addition  to  the  filament  lailure,  the  quartz  lamp  showed  signs  of 
“etching”  inside  Ure  lamp  as  a  result  of  motion  of  the  filament  support  discs  during  ' 
both  vibration  and  acoustic  testing. 

Quartz  lamps  were  also  subjected  to  testing  at  elevated  temperatures.  One  lamp,  which 
had  experienced  etching  as  a  result  of  earlier  acoustic  tests,  was  taken  to  full  rated  power 
without  failure,  indicating  that  the  etching  had  no  effect  oil  lamp  performance.  Another 
quartz  larno  was  operated  at  half  rated  power  and  subjected  to  sinusoidal  vibration  tests. 
The  vibration  input  was  applied  at  the  first  bending  frequency  of  the  lamp,  at  levels  up 
to  42  g’s.  No  failures  occurred  during  tins  test. 

The  tests  described  above  indicate  that  direct  excitation  of  the  quartz  limps  by  the 
acoustic  field  an  be  neglected.  This  observation  is  reasonable  from  a  consideration  of  the 
size  of  the  lamps  (3/8»inch  diameter)  relative  to  the  acoustic  wavelengths  (greater  than  I 
foot)  at  the  lamp  resonances.  At  these  wavetcagths,  the  pressure  gradient  acfo&s  the  lamp 
is  negligible,  so  Hint  no  significant  pressure  differential  appears  across  the  lamp  to  excite 
it.  On  the  other  Imidyvibratioual  to  ihe.Uuup  by  the  rerteo 

tor,  ■■ •  •  •’  •  . 

Tire  ifesgn  described  in  Section  3,2. 1  reduced  the  vibration  in  the  reflect  or.  ' 

Some  vjb.rattonfe  still  present  in  the  reflector,  howevcf,  and  is  transmitted  to  the  lamps 
through  the  lampholders.  Ttwcftm\  the  iampirdldcriisstaliation  itsdlf  was  investigated 
.'a*  aii.pwa  ibi  potential  imtwovemeni. 

Hie  original  quartz  iamp  installatidns  used  in  the  experimental  part  of  this  study  were 
made  with  stauferd  ceramic  LT-ri’K)  lampholders,  as  shown  hi  Figure  13,  An  excessive  num¬ 
ber  of  iihfp  fiiteeft  was  exiwrionced  with  this  installation  when  the  high  intensity  acoustk 
mitee  environment  exceeded  1 60  dB.  One  predominant  type  of  failure  was  a  fracture  of 
the  quartz  en vdopc  a t  the  end  .sal.  an  area  where  the  envelope  makes  a  transition  from  a 
circular  wsm-cbetkm  W  a- -fiat  tab.  It  is  through  this  tab  that  the lead  wire  is  brought  out. 
Another  tyj  tea)  failure  mode  was  fatiguing  of  the  stranded  lead  Wire  just  outside  tt  mve- 
lope.  Both  of  these  types  of  failures  result  from  motion  of  the  lamp  ends  within  the 
ceramic  holder,  since  the  hmps  fit  freely.  in  the  holders  in  a  normal  installation. 


FIGURE  15 

LT-100  LAMP  HOLDER 


Quartz  Lamp 


Modifications  to  the  normal  lamp  installation  in  the  LT-100  holders  were  tried  on 
several  test  runs.  One  modification  was  to  use  self-locking  nuts  on  the  lampholder,  to 
prevent  the  holder  from  becoming  loose  on  the  reflector.  Another  change  was  to  pack 
asbestos  cloth  or  fiberglass  cloth  around  the  lamp  end  seals  to  reduce  motion  of  the  lamps 
witiun  the  holders.  These  modifications  were  successful  in  reduemg  the  frequency  cf  lamt> 
failures.  The  cloth  materials,  however,  could  not  be  reliably  retained  in  the  holders  during 
high  intensity  acoustic  testing.  These  modifications  also  tended  to  change  the  predominant 
mode  of  lamp  failure  to  a  failure  within  the  quartz  envelope. 

A  different  type  of  lampholder  was  evaluated  on  some  lamps.  This  holder,  manufac¬ 
tured  by  Research  Incorporated,  is  identified  as  the  RI  8325-1  holder.  It  consists  of  a 
metal  slip  that  securely  clamps  the  flat  end  of  the  quartz  lamp,  as  shown  in  Figure  16. 

The  8325-1  holder  was  evaluated  in  its  “as-supplied”  condition,  as  well  as  with  asbestos 
or  fiberglass  cloth  wrapped  around  the  endseals,  Some  failures  were  still  experienced  with 
the  8325-1  holder,  live  addition  of  the  asbestos  or  fiberglass  cloth  on  the  endseals  did  not 
improve  the  failure  rate  with  the  8325-1  holder.  In  general,  it  is  concluded  that  the  8325-1 
holder,  retained  by  the  self-locking  nuts,  and  without  any  cloth  wrapping  on  the  endseals, 
provided  the  lowest  lamp  failure  rate.  This  lamp  mounting  method,  combined  with  a  high¬ 
ly  damped  reflector,  is  considered  to  be  the  best  of  the  heater  arrangements  evaluated 
in  this  study. 

Further  experimentation  in  techniques  for  holding  the  lamps  is  still  needed.  One 
approach  that  might  be  pursued  would  be  to  mechanically  uncouple  the  lampholders 
from  the  reflector  This  could  be  achieved  with  vibration  isolators  that  would  isolate 
the  lamps  in  the  frequency  range  where  filament  and  envelope  resonances  exist. 
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SECTION  4 

SPECIMEN  MOUNTING  EFFECTS 

Specimen  mounting  is  one  of  the  most  important  considerations  in  combined  thermal/ 
acoustic  environment  testing.  Even  when  the  acoustic  and  thermal  environments  are  accurate¬ 
ly  controlled,  improper  boundary  simulation  can  iead  to  incorrect  loads  in  the  test  specimen. 
The  specimen  mounting  restrains  one  or  more  of  the  degrees  of  freedom  at  tin  points  of 
attachment,  and  thereby  significantly  affects  the  specimen  response. 

In  the  following  section,  particular  attention  will  be  paid  to  panel  structures.  This  is 
done  for  ease  of  analysis  and  experiment;  however,  the  methods  developed  here  may  be 
extended  to  the  analysis  of  more  complicated  structures.  It  will  also  be  shown  that  in  many 
cases  specimen  response  to  a  combined  thermal/acoustic  forcing  function  may  be  determin¬ 
ed  from  the  superposition  of  responses  to  the  individual  forcing  functions.  Where  super¬ 
position  is  valid,  the  boundary  conditions  in  the  combined  environment  may  also  be  treated 
separately. 

Requirements  for  a  combined  thermal/acoustic  test  generally  fall  in  one  of  two  cate¬ 
gories.  First  one  may  be  required  to  test  a  specimen  to  a  prescribed  response  level.  Second, 
an  environment  (service  or  modeled!  may  be  specified,  in  either  case,  prediction  of  the 
magnitude  and  attribution  of  specimen  stress  is  essential  to  proper  boundary  design  and 
instrumentation  location. 

It  i%  the  purpose  of  this  work  to  examine -the  effects  of  specimen  mounting  on  panel 
state  of  stress  in  a  combined  environment  Both  analytical  and  experimental  techniques 
ore  discussed  and  information  on  free  and  fixed  panels  is  presented.  Following  the  state¬ 
ment  of  the  governing  equations,  response  Jo  thermal  and  dynamic  environments  will  be 
discussed  separately,  Consideration  of  speed  ten  mounting  in  the  combined  environment 
will  then  be  discussed  along  with  recommendations  for  extending  the  methods  developed 
here. 

4.1  GOVERNING  EQUATIONS 

the  analysis  of  the  general  thermal  /dynamic  stress  problem  can  be  accomplished  by 
obtaining  the  solution  to  the  eleven  governing  equations.  These  equations  are  the  continuity 
equation,  the  equations- of  motion,  and  the  energy  equation,  and  the  constitutive  equations; 
winch  are  presented  tti  Appendix  C  as  Equations  (<  i ).  (C2h  (0>.  and  fC'5).  respectively. 

For  purposes  of  later  analysis  and  j  eneraluation  it  will  be  conventient  to  write  the 
governing  equations  in  dimensionless  form.  The  loUowing  substitutions  will  be  made 
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As  explained  in  Appendix  C,  the  assumption  of  negligible  thermal  rates  due  to  me¬ 
chanical  motion  and  negligible  strain  rates  due  to  temperature  changes  leads  to  the  quasi¬ 
static  thermoe*astic  theory.  Under  the  above  assumptions  we  find  that  the  temperature 
distribution  in  a  body  can  be  calculated  from  the  Poisson  equation  for  a  steady  state  sys¬ 
tem.  Letting  q  represent  the  heat  generated  per  unit  volume  per,unit  time  internally  in 
the  body,  the  total  heat  How  to  a  specimen  including  net  absorbed  radiation  and  eon- 
vect’ve  losses  me  be  treated  as  the  amount  of  heat  flow  required  (q)  to  produce  a  given 
temperature  distribution.  The  temperature  is  then  computed  subject  to  the  given  bound¬ 
ary  conditions.  For  the  case  of  the  thin  panels  considered  here,  the  7  term  is  assumed 
to  be  zero.  Thus  the  Poisson  equation  can  be  written  as 
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The  analysis  of  the  state  of  stress  for  a  panel  subjected  to  a  temperature  distribution 
calculated  from  equation  4.2  requires  the  solution  of  the  equilibrium  and  compatibility 
equations.  With  the  usual  assumptions  of  small  deflections  of  a  linea>‘,  elastic  material, 
two  methods  may  be  employeu  to  reduce  the  required  number  of  partial  differential 
equations  to  one.  The  two  functions  which  may  be  employed  aye  the  .Airy’s  stress 
function  (biharmonic)  and  the  strain  (or  displacement)  potential  function  (harmonic).^ 
For  the  case  of  plane  stress  the  state  of  stress  in  terms  of  the  stress  function  is  specified  by 


axxxx  ^  +  2  3xwv  £  +  dyvw  $  =  aE(9v*T  +  dvvT) 


°xx  3yy 


°yy  3xx  &  °xy  ~  '"3xy  ^ 


(4.3a) 

(4.3b) 


The  state  of  stress  in  terms  of  the  strain  potential  for  plane  stress  is  given  in  dimension¬ 
less  form  by 

2 

^*il)  Syy*“|)(I+")«(TcT+TB)+C  (4A) 

It  is  assumed,  in  equation  (4.4a)  that  V*  \p  and  T  are  differentiable  functions  of  x  and  y.  The 
strain  potential  and  stress  relations  are  then  given  by 


b  a  — 

exx  “  ~  3xx  €yy  “  7  3yy  exy  ~  ^axy  ^ 


(4.4b) 


and 
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(4.5a) 

(4.5b) 

(4.5c) 


it  can  be  shown,  by  applying  the  substitutions  indicated  in  equation  (4.1 )  to  (4.3a),  that 
equation  (4.3a)  is  the  lapacian  of  equation  (4.4a).  For  ease  of  numerical  solution  in  the 
present  discussion  solutions  will  be  obtained  using  the  strain  potential. 


The  in-plane  forces  on  a  panel  may  now  be.determined  by  substitution  of  equations 
(4.5)  into  the  relations 


v  v 

r2  ri  n 

Sex  ^  axx  az>  ^yy  j  °yy  az>  ^xy  ~J^C  a> 


Hie  dynamic  plate  equation  (equation  C2  i )  may  be  written  in  dimensionless  form  as 
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where  the  \p  values  are  known  from  equation  (4.4a). 

On  the  basis  of  the  preceding  equations,  (4.1 )  through  (4.7),  the  solutions  to  the 
thermal  and  dynamic  stress  problems  may  be  determined  for  a  given  dynamic  forcing 
function  and  damping  value.  Due  to  the  complexity  of  the  problem  numerical .  "  itions 
to  the  thermal  stress  problem  were  obtained  principally  for  free  and  fixed  boundary  condi¬ 
tions  on  a  square,  steel  plate.  Experimental  solutions  were  obtained  for  a  boundary  condition 
between  free  and  fixed.  A  determination  of  the  effect  of  in-plane  forces  on  the  resonant 
frequencies  of  a  panel  was  obtained  experimentally.  Methods  for  extending  this  information 
to  other  boundary  conditions,  shapes,  and  materials  are  presented  in  sections  4.2  through  4.4. 

4.2  THERMAL  EFFECTS 

The  governing  equation  used  in  the  study  of  thermal  stress  in  panels  is  equation  (4.4). 
Since  this  equation  is  derived  by  the  integration  of  two  equilibrium  equations  an  arbitrary 
constant  is  included.  If  the  potential  and  its  derivatives  are  taken  to  be  zero  at  infinity  the 
constant  is  zero.  There  are,  however,  classes  of  practical  problems  which  can  be  solved  con¬ 
veniently  by  choosing  a  reference  for  the  potential  function  at  a  point  in  space  other  than 
infinity.  In  these  cases  it  will  be  necessary  to  retain  and  evaluate  the  constant  of  integration, 

C,  in  equation  (4.4). 

Evaluation  of  the  constant,  C,  requires  the  evaluation  of  the  expression 
/a'\2  /a\ 

3*5*  +7  avv^“  7  (l+»>)aT=C  (4.8) 
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at  some  point  x0,  yQ.  To  do  this  equations  (4.4b)  are  substituted  into  equation  (4.4a)  and 
evaluated  at  xn,  yn 


e*s|x=x0  *eyv  |x = x0  ^<1+l')“Tjs=x0  =c 
y=y0  y  =  yo  ?  =  y0 


The  following  discussion  will  be  limited  to  square  free  and  fixed  panels  with  a  constant 
boundary  temperature.  It  will  be  shown  later  that  equation  (4.8)  may  be  applied  to  more 
general  problems  but  s^me  additional  calculation  may  be  required.  For  the  fixed  panel  it  is 
known  that  at  the  comers  exx  =  Cyy  =  0  and  the  temperature  is  Tg.  Substituting  these 
values  into  equation  (4.9)  yields 


C=  —  (l-b')aTg 


(4.10) 


For  the  free  panel  exx  and  are  both  equal  to  oTg  at  the  corner  and  the  temperature  is 
Tg.  Equation  (4.9)  then  becomes 


C-  (l-i>)  aTg 


(4.11) 


The  second  part  of  the  boundary  value  problem  requires  the  evaluation  of  y;  at  the 
boundaries.  The  general  form  of  ^  (x,  y)  is 


0  (x,  y)  =  f j  (x)  +  f2  (y)  +  f3  (x)  f4  •?)  +  C 


(4.12) 


For  the  free  panel  we  have  the  requirement  that  the  shear  stresses,  anu  thus  the  shear  strains, 
are  zero  at  the  boundaries  (equation  (4.4b)). 


hy*  -  0 


(4.13) 


Substituting  equation  (4.12)  into  the  above  expression  and  evaluating  this  at  the  boundaries 
1x1  =  I  and  lyl  "  1  we  have  the  condition 


aXy  i//  =f'3(i)f,4(y)  =  f'3(x)f'4a)=o 

B 


(4.14a) 


which  gives 


f'3(x)  =0,  f'4(y)  =0 

J  lyi  =  1  1x1  =  1 

-1  <  Ixl  <  I  -1  <  lyK  1 


(4.14b) 


t  =fj  (x)  +  f2(y)vC 
B 


(4.14c) 


At  the  boundary  the  normal  stresses  must  be  zero.  In  terms  of  the  potential  function  equa¬ 
tions  (4.5)  become 


7  hx$  Z  hy  $  -(Hv)aT(l.y)  =  0 

b  1x1  =  1  b  yy  1x1  =  1  • 


(4.15a) 


-  any  *  +»  “  hx  *  -•(l+t»)aT(x,  1)  =  G 

b  yy  lyl  =  1  a  XX  lyl'  =  l 


(4.15b) 


Substituting  equation  (4.14c)  into  the  above  we  have 


£  f"l  0)  +  »~  f"2(y)  -  (Hv)  «T  (1,  y)  (4.16a) 

“f'^tD  +  y-f'^Cx^d+^fiTCx,  I)  (4.16b) 


Solving  equation  (4.16)  for  f  j  and  >  at  the  boundaries  we  have  for  the  constant  boundary 
temperature  case 


fi  (x)  = 

f2(y)= 


a(l-n^) 
2  hv 

b(l-H^) 

2av 


ttT  (x,  1)- 

«T  (1,  y)  - 


ff'd) 

2v  . 
f/d) 
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x  +  C|  x  +  C2 

y  +  C3  y + c4 


(4.17a) 


(4.17b) 


where  fj"  (i),  f2  "  (1 ),  Cj ,  and  Co  are  unknown  constants  to  be  evaluated  from  the. 
boundary  conditions  for  the  known  temperature  distributions.  C2  and  C4  are  arbitrary 
since  the  potential  is  only  known  to  within  a  constant.  Solving  equation  (4.17)  for  the  case 
of  a  square  panel  with  a  constant  boundary  temperature,  Tg,  we  have  upon  substitution  of 
the  second  derivatives  into  equation  (4.15) 

f"1d)  =  f"2(l)  =  d-^)aTB  (4.18) 


Applying  the  condition  of  zero  rigid  body  motion  we  have  Cj  -  C3  =  0  in  equations  (4  17). 
Thus  the  potential  is 


ip  (x,y) 


=  aTB  (x2  +y2) 

B 


+  c2  +  c4 

b 


(4.19) 


Since  \j/  (x,y)  can  only  be  known  to  within  an  arbitrary  constant  we  can  assign  a  value 
to  a  point..  In  this  case  we  will  assign  a  zero  value  to  the  midpoint  on  a  boundary.  This 
results  in  a  value  of  -aTg/2  for  (C4  +  C2). 


For  the  fixed  panel  the  boundary  conditions  are  that  the  displacements  at  the  bound¬ 
ary  are  zero. 


95  * 


lyl  =  1 
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|xi=l 


(4.20) 


Applying  these  conditions  to  equation  (4.12)  we  have 

f'j  (x)  +f'3  (x)  f4  (1)  =  0  (4.21a) 


f'2(y)+f3(l)f'4  (y)  -  0 


(4.21b) 


Integrating  and  combining  like  functions  results  in 

fj  (x)+f4(l)f3(x)-Cl  =0  i5  (x)=  Cj  (4.22a) 

f2  (y) + f3  0 )  f4  (y>-  C2  =  0  .’.  r6  (y)=  C2  (4.22b) 


stsBangfete 
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where  f3  (x)=  fj  (x)+  f4  (1)  f3  (x)and  similarly  for  fg  (y).  Since  Cj  and  C2  are  arbitrary 
constants  they  can  be  given  zero  value  yielding 


(x,  y) 


(4.23) 


For  the  free  and  fixed  square  panel  with  constant  boundary  temperature  we  can  now  deter¬ 
mine  a  function  ip  which  satisfies  equation  (4.4a)  and  the  boundary  conditions  given  by 
equations  (4. 10)  and  (4.23)  or  (4.11)  and  (4. 1 9).  From  potential  theory,  it  can  be  shown 
that  with  p  known  on  the  boundary,  the  function  obtained  unde/  the  above  conditions 
is  the  unique  solution.^  There  are  several  methods  for  the  solution  of  equation  (4.4a).  For 
simple  temperature  distributions  and  classical  boundary  conditions,  closed  form  analytical 
solutions  are  available.  For  general  temperature  distributions  and  various  boundary  condi¬ 
tions  two  other  methods  of  solution  seem  preferable.  One  method  is  the  numerical  analysis 
of  equation  (4.4a).  The  other  method  consists  of  a  combination  of  experimental,  numerical, 
and  statistical  techniques. 


Three  techniques  may  be  employed  in  the  numerical  analysis  of  equation  (4.4a): 
numerical  integration,  finite  elements,  or  finite  differences.  In  the  present  study  the  finite 
difference  technique  was  employed  due  to  its  compatibility  with  the  low  storage  time  shar¬ 
ing  computers  and  the  existence  of  several  applicable  algorithms  for  the  solution  of  elliptic 
equations  by  finite  differences.  As  shown  in  the  literature,^  the  solution  to  the  finite  dif¬ 
ference  representation  of  the  second  order,  elliptic  equation  converges  uniquely  to  the 
solution  of  the  partial  differential  equation.  The  finite  difference  solution  is  the  super¬ 
position  of  the  homogeneous  and  particular  solutions.  A  description  of  the  computer 
programs  used  in  solving  this  equation  by  the  finite  difference  technique  is  given  in 
Appendix  D.  Once  the  program  has  been  developed  solutions  may  be  obtained  quickly 
and  economically;  however,  each  case  requires  an  individual  solution.  It  is  often  desirable 
to  have  an  indication  of  the  variation  of  the  stress  for  a  range  of  independent  parameters. 

In  this  case  a  few  computer  and  experimental  solutions  may  be  combined  with  any  of 
several  statistical  techniques  to  yield  design  curves  giving  a  reasonable  approximation  to 
the  actual  stress  response.  *  * 

The  form  of  T  for  the  nonhomogeneuus  portion  of  equation  (4.4a)  will  be  taken  as 


T (x,  y)  =  Aj  +  A2lx|  +  A3lx2|  +  A4Ix3|  +A5Ix41 

-1.  <x<  1.  -  jx|  <y  <  |x| 

(4.24a) 

T  (x,  y)  =  Aj  +  A2|  y|  +  A3  |y2|  +  A4  jy3|  +  A5  |y4| 

-1.  <y<  1.  -  | y  1  <x<jy| 

(4.24b) 

This  form  of  temperature  distribution  is  symmetric  with  respect  to  the  x  and  y  axes  and 
has  a  constant  boundary  temperature.  Forming  the  dimensionless  temperature  ratio  T 
equation  (4.24a)  becomes 

T  (X,  y)  =  A2  ( x  -1.)  +  A3  ( -1 .)  +  A4  (|xj3-I)  +  A5  (|X|i-L) 

-l.<  x  <1.  - |x«y<  | S|  (4.25) 

where  A2  =  A2/Aj  ,  A3  =  A3/Aj  ,  etc.  A  similar  expression  for  equation  (4.24b)  can  be 
written.  The  normalized  center  temperature  is  then 

Tc  =  ~^2  -  A3  -  A4  -  A5  (4.26) 
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Since  stress  is  the  response  parameter  of  interest  we  would  like  to  find  o  =  a  (x,  y,  T).  Mak¬ 
ing  the  substitutions  indicated  in  equation  (4.1 )  in  the  expression  for  stress  (equation  (4.5)) 
we  define  a  dimensionless  response  parameter  as 


cxx 


xx 


b  -  a 
_  a  dxx  ^  +  v°  3yy  I jl 
oEX  (1  +  v)  o:T 
\-v 


- 1 


(4.27) 


with  similar  expressions  for  3pp.  Since  \J/  depends  on  the  coefficients  of  the  non-homogeneous 
portion  of  equation  (4.4a)  and  the  boundary  conditions,  equation  (4.27)  can  be  written  in 
parametric  form  as 

—  a 

a  =  a  (A;-,  v,  boundary)  (4.28) 

1  b 


Particular  attention  will  be  paid  to  the  center  stress  parameter,  oq  as  a  function  of  Aj  for 
a/b  =  1 .0,  v  =  0.27,  for  free  and  fixed  boundaries. 

The  initial  intent  of  this  study  was  to  apply  the  statistical  technique  known  as  Latin 
Squares  to  determine  the  relation  between  a  and  T,  ^  However,  a  few  solutions  to  the 
thermal  stress  problem  indicated  that  for  the  above  type  of  symmetric  temperature  distri¬ 
butions  and  boundary  conditions,  oq  varies  linearly  with  each  coefficient,  A:,  in  this  case 
the  use  of  a  tool  as  powerful  as  Latin  Squares  is  not  required.  For  problems  where  the 
functional  relationships  are  more  complicated  the  Latin  Squares  technique  should  be 
employed.  The  value  of  oq  for  a  temperature  distribution  which  is  a  function  of  more 
than  one  coefficient  is  found  by  the  superposition  of  the  solutions  to  the  various  single 
coefficient  problems.  Computer  solutions  of  Oq  for  various  values  Aj  for  free  and  fixed 
panels  are  presented  in  Tables  3  and  4,  respectively.  The  values  and  locations  of  the 
maximum  stresses  are  also  presented.  For  each  case  where  the  temperature  distribution 
was  purely  linear,  quadratic,  cubic,  or  quartic  a  least  squares  polynomial  fit  of  oq  vs  Aj 
was  performed.  In  each  case  the  form  of  the  polynomial  was  oq  -  CjAj  +  Dj.  For  the  case 
of  a  fixed  plate  the  Cj  are  all  0.626  and  the  Dj  are  -1.0.  For  the  free  plate,  Cj  =  -0.372 
and  Dj  =  0.  As  noted  in  equation  (4.26),  f  c  =  Ej~Aj  so  each  plot  for  T  having  oniy  one 
non-zero  coefficient  is  also  a  plot  of  oq  vs  T.  lire  oq  vs  fc  cum  for  both  the  free  and 
fixed  panel  are  presented  in  Figure  17.  It  should  be  noted  that  for  T  ®  0,  the  constant 
temperature  case,  there  is  zero  stress  in  the  free  panel  and  oq  ~  •••  oriu/l  -v  for  the  fixed 
panel.  When  T  =  1.0,  Ty  -  ambient,  the  panel  responds  as  a  fixed  panel  regardless  of 
physical  bounuary  conditions  and  oq  is  --  (I  ~v)J 2.  Cases  i  1  and  1 2,  Table  4,  for  the 
flxedjpanei  and  cases  14  through  16,  Table  3  for  the  free  panel  verify  that  applying  T  - 
E  -  Aj  and  the  equations  for  d  vs  f  yields  the  same  3  value  as  the  solution  of  equation 
(4.5). 


The  results  of  the  analytical  studies  shown  in  Tables  3  and  4  indicated  that  for  tempera¬ 
ture  distributions  of  quadratic  and  higher  order  the  maximum  stress  is  not  in  the  center 
Temperature  and  stress  distributions  for  the  fixed  and  free  panels  for  first  through  fourth 
order  temperature  distributions  are  presented  in  Figures  1 8  through  25,  The  shape  of  the 
stress  distribution  shown  for  the  higher  order  distributions  was  also  observed  in  the  ex¬ 
perimental  investigations. 

In  addition  to  the  analytical  studies  a  thermal  test  was  conducted  on  a  15-inch-square, 
stainless  steel  panel.  The  te*  t  setups  used  for  the  simulated  free  and  fixed  boundary  tests  are 
shown  in  Figures  26  and  27,  respectively.  Heat  was  applied  to  the  panel  by  the  quartz  lamp 
and  reflector  arrangement  shown  in  the  figures. 
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FIGURE  19 

THERMAL  STRESS  DISTRIBUTION 
FIXED  BOUNDARY  CONDITION 
2nd  Order  Temperature  Distribution 
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FIGURE;  22 

THERMAL  STRESS  DISTRIBUTION 
FREE-FRL*.  BOUNDARY  CONDITION 
1st  Order  Temperature  Distribution 


Thermal  Distribution 


. 


FIGURE  24 

THERMAL  STRESS  DISTRIBUTION 
FREE-FREE  BOUNDARY  CONDITION 
3rd  Order  Temperature  Distribution 
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FIGURE  25 

THERMAL  STRESS  DISTRIBUTION 
FREE-FREE  BOUNDARY  CONDITION 
4th  Order  Temperature  Distribution 
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Temperature  and  strain  data  were  obtained  at  several  points  on  the  fixed-edge  test  panel. 
A  least  squares  polynomial  fit  was  performed  on  the  centerline  temperatures  and  the  coeffi¬ 
cients  were  normalized  to  obtain  the  Aj  values  described  in  equation  (4.26).  The  value 
calculated  in  this  manner  is  compared  with  T^  obtained  as  (T^  -  Tg)/T^  using  an  average 
boundary  temperature.  The  results  are  summarized  in  Table  5.  During  the  initial  low  tempera¬ 
ture  thermal  tests  the  stress  data  for  the  experimental  panel  was  consistent  with  the  theoretical 
values  shown  in  Figure  17.  During  the  high  temperature  tests  the  experimental  boundary  was 
substantially  less  stiff  than  a  fixed  condition  as  shown  by  the  fact  that  the  strain  at  the 
boundary  was  approximately  aTg/2.  The  boundary  strain  for  a  fixed  panel  would  be  zero 
while  the  boundary  strain  for  a  free  panel  would  be  aTg.  A  least  squares  analysis  of 
Qq  vs  Tq  for  the  high  temperature  tests  on  .the  experimental  panel  is  presented  in  Table  6. 

The  results  of  this  analysis  are  shown  in  Figure  28. 

The  techniques  used  in  this  section  for  determining  the  stress  distribution  due  to  thermal 
effects  can  be  powerful  tools  in  analyzing  actual  laboratory  cases.  As  shown  previously,  for 
simple  experimental  cases,  a  good  estimate  of  the  actual  stress  experienced  by  the  test  panel 
can  be  obtained  by  utilizing  normalized  polynomial  temperature  distributions.  A  few  experi¬ 
mental  cases  may  then  be  run  to  determine  Oq  vs  T q.  For  variable  boundary  temperatures 
and  difficult  boundary  conditions  equations  (4.4a)  and  (4.9)  must  be  solved  for  C  and 
’/'boundary*  respectively  and  these  values  used  in  the  computer  programs  shown  in  Appendix 
D.  If,  however,  the  boundary  can  be  assumed  to  be  close  to  one  of  the  conditions  studied  here, 
Figure  28  may  be  used  to  give  an  approximate  center  stress.  Tablfes  3  and  4  also  give  an  indi¬ 
cation  of  how  maximum  stress  varies  with  increased  order  of  temperature  distribution. 

To  verify  the  accuracy  of  the  technique  described  in  this  section  a  comparison  was 
made  between  the  closed  form  solution  of  the  rectangular  plate  problem  (shown  in  Figure 
29  and  described  in  reference  14)  and  the  solution  obtained  by  the  finite  difference  solution 
of  the  strain  potential  equation.  The  plate  in  this  case  is  free  to  expand  in  the  y  direction. 

The  boundaries  at  x  =  +  b/2  are  placed  in  the  “free”  condition  by  the  application  of  a  stress 
in  the  x  direction  of  2aET0/3,  according  to  Timoshenko. 

At  the  corners  we  have  the  conditions 


°xx  =  2«ET0/3,  Oyy  =  0.  ,Tcorner  =  0. 


(4.29) 


Solving  equations  (4.5a)  and  (4.5b)  for  exx  and  eyy  and  substituting  into  (4.9)  the  constant, 
C,  is  found  to  be  2(1  ~v)<xTJ3.  We  assume  to  be  of  the  form  given  in  equation  (4. 1 2). 

The  values  at  the  boundary  then  become 


0 


x  =  +^- 


+  f2(y)  +  f3  ~  f4(y) 


(4.30a) 


* 


=  +a  =fj  (x)  +  f2(“)+  f3  (x)  14(^) 


(4.30b) 


The  boundary  conditions  are 

axx  =  2aET0/3,  ayy  =  0 


at 


b 

x  =  ±  ~,y  =  ± 


a 

2 


(4.31) 
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TABLE  5 

EXPERIMENTAL  EVALUATION  OF  Tc 
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TABLE  6 

EXPERIMENTAL  THERMAL  STRESS  DATA  ANALYSIS 
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- 

7 
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x 

*  T 

MEAN  VALUE  2F  Y 

= 

- .  36428  6 

STD  ERR0R  0F  Y 

3.41258E-C2 

Y 

■  O' 

P0LYFIT  0F 

DEGREE  l  INDEX  OF  DETER U  - 
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1 

• 1 .4126 1£-C£ 

X-ACTUAL 
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ulFF 

PCT-D1FF 

*354000 

.  -.340000  36.2032 
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- . 40  t  000  -.362597 
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.406000 

-.301000  -.362766 

6. 17661  E-OS 

-17.0264 

.453000 

-.390000  -.363430 

-2  .65699E-02 

7.21 088 

*460000 

-.383000  -.363529 

-1 .94710E-02 

5.35612 

.528000 

-.370000  -.364490 

-5.51 O47E-03 

1.51183 

s 

i 

-.365000  -.371157 

6. 15707E-03 

-1.65889 

STD  ERR0R  Z-  ESTIMATE  F0R  Y  *  3.72253E-02 
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Substitution  of  equations  (4.30)  and  (4.3 1 )  into  (4.5a)  and  (4.5b)  will  allow  solution  for 
the  functions  fj ,  {2,  fj,  and  f^  to  within  a  constant.  Applying  the  requirement  that  the 
values  of  \p  given  by  equations  (4.30a)  and  (4.30b)  be  the  same  at  the  corners  allows 
evaluation  of  the  arbitrary  constants.  The  values  of  ip  at  the  boundaries  are  then  given  by 


^  >■  (-$) 


va?o 

3 


(4.32a) 


5(l4^)aT0  ya:T0 
48  +  I2~ 


* 


aTox2 

3 


»T0b^ 

12 


(4.32b) 


Solution  of  equations  (4.4)  through  (4.5)  subject  to  the  constant  in  (4.4)  being 
2  (1  -P)al0l3  was  accomplished  by  use  of  the  computer  programs  presented  in  Appendix 
D.  Table  7  lists  the  results  of  this  analysis  compared  to  the  closed  form  solution. 


In  this  case  TQ  was  taken  as  5387. 9°F,  thus  for  a  steel  panel  a  ET0  becomes  1 .00  x  I  (A 
The  error  in  the  values  of  oxx  are  obviously  due  to  truncation  of  the  decimal  values  of 
improper  fractions.  The  error  in  is  due  to  the  method  of  computing  strain  which  re¬ 
sults  in  taking  the  difference  of  large  numbers.  The  accuracy,  if  required,  can  be  improved 
by  using  double  precision  arithmetic,  increasing  the  number  of  iterations  in  the  potential 
program,  and  using  a  finer  mesh  in  the  finite  difference  routine. 


Substituting  the  closed  form  solution  (4.33)  into  equations  (4.5a)  and  (4.5b)  and 
solving  for  strain  we  find  that  the  sum  of  the  strains  minus  (l  +p)  atTat  any  point  in  the 
plate  is  2(1  -v)«T0/3.  This  is  the  constant  of  integration  evaluated  from  equation  (4.9) 
and  the  boundary  conditions  at  the  corners. 


For  commonly  encountered  boundary  conditions  and  materials  it  is  recommended 
that  the  techniques  developed  previously  be  used  with  a  statistical  method  such  as  Latin 
Squares  to  extend  the  design  information  presented  here.  For  details  on  the  use  of  Latin 
Squares  the  reader  is  referred  to  references  1 1  and  13.  It  is  suggested  that  a  four  level 
Latin  Square  plan  be  used  to  determine  0 ^  and  5^a>£  as  a  function  of  t .  a/b.  and  bound¬ 
ary  stiffness.  The  variation  of  o  with  the  parameter  v  need  not  be  considered  since  it  is 
approximately  the  same  for  most  engineering  materials  and  E/l  is  contained  in  the 
parameter  a,  Tire  boundary  stiffness  may  be  determined  either  qualitatively  or  quantitatively. 
'Hie  stiffness  used  in  the  analysis  must  be  n^n-dirnensionalized  to  an  appropriate  value.  An 
example  of  a  four  level  Latin  Square  design  for  determining  a  as  a  function  of  boundary. 

T,  and  a/b  is  shown  in  Figure  30.  This  analysis  requires  4 2  or  16  experiments  as  compared 
to  4^  or  64  experiments  for  a  conventional  analysis  of  variance. 


The  values  of  Oj  j ,  5j  3.  and  oj^  are  obtained  from  Figure  28.  The  values  of  3->|  .031. 

1  •  0^4.  034.  and  044  may  be  obtained  by  using  computer  programs  presented  in  Appendix 
D.  The  remaining  o1(  may  be  determined  by  experiment  or  numerical  solution  as  described 
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in  references  U  a-',.  (3.  The  partial  responses  oc  =  f  (T),  oc  =  f  (K-r),  oQ  =  f  (a/b),  and  the 
total  response.  oc  -  f  (T,  Kj,  a/b)  can  be  determined  from  statistical  techniques.  If  it  is  de¬ 
sired  to  find  5jy?ax  as  a  function  of  T,  Kj,  and  a/b  the  form  of  T  may  be  another  independ¬ 
ent  parameter  in  which  case  a  Greco-Latin  Square  would  be  used.  This  would  require  4-^  or 
64  experiments  instead  ot  c.5o  tor  the  conventional  analy^.s  of  variance. 

Tlie  results  of  the  above  analysis  would  provide  the  design  curves  necessary  for  thv 
selection  of  boundary  stilfness  to  induce  the  proper  stress  level.  For  a  given  b^’  r.uary  an 
estimate  ot  the  center  and  maximum  stress  for  a  given  temperature  distribution  could  also 
be  obtained  from  the  design  curves. 

4.3  DYNAMIC  EFFECTS 

Hie  solution  ot  the  dynamic  stress  problem  requires  the  solution  of  the  following 
equations: 


“rx “ :,T7i  K’’ 


(4.34a) 


=  __  (  v->  0V-  ws  r  r  i  -)  a-  w 

2(1  r)  h  }y  -i  Xx  51  j 


(4.34b) 


°xr2(';^V  ;ixy  *s 

Ihus  tor  known  material  parpertios  the  stress  can  be  determined  when  the  deUec 
is  known. 


(4.34c) 


tiou,  w 


Die  first  step  in  obtaining  the  deflection.  vy.,  requires  the  solution  of  the  dynamic  plate 
equation  (4.  / 1.  The  values  of  N'xx.  N^..  Nxv  are  given  by  equations  (4:5)  and  (4,n),  h  t>. 
evident  that  the  thermal  stress  problem  must  be  solved  prior  to  solving  tire  combined, 
dynamic  and  thermal  problems.  From  equation  (4.7)  the  parameters  that  affect  (he  undamp¬ 
ed  mode  shape  tor  free  vibrations,  w,  arc  boundary  conditions,  a  ■  b ,  a  c,  a  T  i*.  and  ipj  ^  u,>t  -  j 
igt;  C*>.  Hus  last  parameter  is  of  spevu!  interest  die  to  the  w  ‘  (  nn  In  general  w.  2  is  " 
ns>t  known  atm  itt  fact  may  be  considered  .vs  a  dependent  variable  which  is  a  function  of  the 
mode  shape.  m- plane  forers.  boundary  conditions.  geometric  and  material  properties.  It 
wdl  he  convenient  to  multiply  the  dimensionless  parameter  containing  w  2  by  I  2i  1  r2i  to 
yield  u>?l  igl.  C  I  -•  1  ■"" )  f>  a  .)  Hie  term  in  brackets  van  be  considered  as  the  natural 
frequency  of  the  panel  at  ambient  temperature  divided  by  a  constant.  Thus  we  write  this 
dimensionless  parameter  as  ^“‘Ow,,2  ambient.  Hie  value  of  C2  for  various  boundary 
conditions,  materials,  and  a/b  ratios  >s  readily  available  m  Uie  literature.  •  ^6 

Eor  zero  in-plane  forces  <on  ts  the  frequency  Risen  by  equation  (Cl 9 1  and  for  increas¬ 
ing  compressive  loads  u.-n  is  decreased  to  zero  3t  She  Orxf  critical  buckling  load.  Since  equa¬ 
tion  (4.7)  is  a  homogeneous  equation  the  previously  employed  finite  difference  technique 
yields  only  the  trivial  solution  w  -  0  However,  a  number  of  computer  program.*  arc  avail¬ 
able  in  the  industry  for  computing  u:n  and  w  for  various  >;nicnirai  mo  lev  Most  of  these 
programs  require  only  the  geometry  and  material  properties  as  input  A  lumped  parameter 
nvcsx  and  stiffness  matrix  is  then  generated  using  techniques  such  as  the  finite  dement 
technique.  The  values  of  w  arc  then  obtained  as  (he  eigenvectors  and  u?n  as  the  eigenvalues 
s>f  the  dynamic  matrix  by  one  of  several  numerical  techniques. 


J t: 


The  main  problem  with  this  type  of  solution  is  accounting  for  the  thermal  loads. 
Equation  (C20),  Appendix  C,  rewritten  in  the  form  Mdtt|X,+KxpO  becomes 
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yyi-yy  -  j) 
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From  ihis  it  is  apparent  that  the  expression  in  braces  is  a  stiffness  term.  Thus  it  would  be 
necessary  to  modify  finite  dement  programs  such  as  the  McDonnell  Douglas  Automation 
Company  dynamic  analysis  program  to  provide  for  a  reduction  of  the  terms  in  the  stiffness 
matrix  !  .  a  matrix  determined  from  t**  .us  involving  N-:  3y  w.  It  may  also  be  possible 
to  use  a  modified  form  of  the  current.  puiar  NASTR.4N  program. ' ' 

Although  the  above  modified  programs  would  provide  a  solution  to  Urn  dynamic 
equation  including  in-plane  loads,  it  would  *,t«!l  be  necessary  to  use  a  large  number  of 
elements  to  obtain  reasonable  accuracy.  .  a  the  experimental  1 5  inch  square  panel,  if 
elements  are  taken  every  inch  (fairly  large  increments)  the  stiffness  and  mass  matricies  are 
1 5  x  1 5  for  only  one  degree  of  freedom.  The  size  of  cue  matrix  is  increased  by  !  5  per  side 
for  each  addition?'  degree  of  freedom.  It  is  obvious  that  the  storage  limitations  and  cost 
soon  become  limiting  factors.  Single  runs  of  the  above  mentioned  programs,  not  modified 
to  include  thermal  loading,  may  cost  several  hundred  dollai.,. 

The  cost  involved  in  the  numerical  solution  and  the  lack  of  closed  form  solutions  for 
equation  (4.7)  indicates  that  either  experimental  solutions  or  combined  experimental-and 
numerical  solutions  should  be  employed.  Mapping  mode  shapes  in  the  combined  thermal/ 
acoustic  environment  can  be  an  expensive  and  time-consuming  process,  thereby  making 
experimental  determination  of  w  difficult.  One  possible,  low  cost  method  of  obtaining  de¬ 
sign  curves  for 


0  9 

n  /CZwri  vs  independent  parameters  is  to  experimentally 


determine  the  values  of  to 


for  various  high  temperature  cases.  In  some  cases, 


such  as  simply  supported  polygonal  plates  with  constant  in-plane  forces,  the  mode  shape 
is  independent  of  the  magnitude  of  the  in-plane  load.  Thus  a  computer  solution  or  experi¬ 
mental  mode  mapping  is  not  necessary.  In  cases  where  the  mode  shape  is  not  known  from 
the  literature  it  may  be  necessary  to  perform  a  ground  vibration  test  to  determine  the 
modal  response.  Normalized  mode  shapes  can  then  be  computed  from  the  measured  data 
and  the  mass  matrix.  Equations  (A4)  and  (A5)  relate  the  normalized  mode  shapes  and,  the 
natural  frequencies  to  the  mass  and  stiffness  matricies.  It  is  these  matricies  which  must  be 
known  for  the  elevated  temperatuio  case.  The  generalized  mass  and  stiffness  matricies  are 
then  given  by 

T-1 

(M)  =  10] 1  [0]  (4.36a) 

IK]  «[01T-1r«n?.J  [0]_1  ’-(4.36b) 

It  should  be  noted  here  that  two  important  assumptions  have  been  made;  first,  it  is  assum¬ 
ed  that  the  method  of  excitafion  and  measurement  does  not  effect  the  normal  modes; 
second,  the  damping  must  be  small  so  that  the  m<“!e  shape  computed  from  tho  measured 
data  is  the  same  as  the  undamped  normal  mode.  It  is  recommended  that  measured  modes 
be  checked  for  orthogonality  to  test  the  accuracy  of  the  above  assumptions. 
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Once  the  system  characteristics  have  been  determined  either  by  one  of  the  finite 
element  computer  programs  or  the  experimental  methods  described  above  the  deflection 
can  be  determined  for  a  given  forcing  function  and  damping  value  by  the  methods  present¬ 
ed  in  section  2.2.  This  section  presents  methods  for  determining  the  deflections  which  must 
be  substituted  into  equations  (4.34)  to  obtain  stress  in  terms  of  the  coupling  between  the 
structure  and  the  given  acoustic  field. 

In  the  present  experimental  study  a  plate  with  both  free  and  fixed  boundary  condi¬ 
tions  was  subjected  to  acoustic  loads  at  ambient  temperature  apd  at  a  nominal  450°F  plate 
center  temperature.^  The  experimental  test  setups  for  the  free  and  fixed  panfels  are  the 
same  as  those  for  the  static  thermal  tests  (Figures  26  and  27).  Mathematical  models  used 
in  the  ambient  temperature  modal  analysis  for  the  free  and  fixed  panels  are  shown  in 
Figure  31 .  The  mode  shapes  and  natural  frequencies  from  the  analysis  and  experiment 
at  ambient  temperature  are  presented. in  Figure  32  for  the  fixed  and  free  panel.  The 
experimental  frequencies  shown  are  the  resuits  of  vibration  excitation  since  some  modes 
could  not  be  excited  acoustically  and  mode  mapping  was  not  convenient  in  the  acoustic 
environment.  However,  these  results  give  an  indication  of  the  quality  of  the  dynamic 
boundary  conditions.  Frequency  shifts  during  acoustic  testing  from  ambient  to  450°F 
center  temperature  are  shown  in  Figures  33  and  34  for  the  free  panel  and  Figures  35  and 
36  for  the  fixed  panel.  Table  8  summarizes  the  frequency  changes  and  the  percentage 
frequency  shifts.  As  shown  in  the  table  the  average  reduction  in  natural  frequencies  was 
3.4%  for  the  free  panel  and  9.4%  for  the  fixed  panel.  A  reduction  of  5.0%  in  Young’s 
Modulus  due  to  temperature  would  account  for  a  2.2%  reduction  in  natural  frequency. 

Thus  a  1 .2%  and  7.2%  decrease  in  frequency  in  the  free  and  fixed  panels  can  be  attributed 
to  the  effect  of  in-plane  forces. 

4.4  THERMAL/ACOUSTIC  MOUNTING 

The  preceding  sections  have  shown  how  thermal  stresses  and  dynamic  stresses  can  be 
analytically  evaluated  in  a  sonic  fatigue  test  specimen.  The  importance  of  the  support 
structure  for  combined  thermai/acourric  loading  is  apparent'  from  the  governing  equa¬ 
tions.  In  designing  a  mounting  for  testing  panels  in  the  combined  environment,  one 
must  consider  both  the  thermal  conditions  at  the  boundary  of  the  panel,  and  the 
mechanical  restraint  imposed  on  the  edges  of  the  panel.  Theremainder  of  this  section 
will  discuss  specimen  mounting  designs,  starting  with  the  thermal  aspect  of  the 
problem.  After  the  method  for  controlling  the  boundary  temperature  has  been 
established,  one  may  proceed  to  the  design  of  a  mechanical  boundary  that  will  produce 
the  desired  thermally-induced  in-plane  loads.  Finally,  attention  may  be  directed  to  the 
design  of  panel  edge  repaints  that  will  result  in  the  proper  panel  frequencies  and  mode 
shapes. 

The  first  part  of  the  boundary  synthesis  problem  is  the  design  of  the  thermal  boundary. 
The  specification  of  boundary  temperature  and  qwill  determine  the  temperature  distribu¬ 
tion.  Maintaining  the  proper  boundary  temperature  may  be  as  important  as  the  physical 
stiffness  of  the  boundary.  As  described  in  section  4.2,  even  a  free  panel  with  boundaries 
maintained  at  ambient  temperature  will  have  the  same  thermal  stress  as  a  fixed  panel. 
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FIGURE  31 

FREE*PANEL  CONFIGURATION  FOR  ANALYTICAL  STUDIES 


FIGURE  32 

FIXED  PANEL  MODAL  CHARACTERISTICS 
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FREE  PANEL  MODAL  CHARACTERISTICS 
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FIGURE  33 

ACCELERATION  RESPONSE 
FREE-PANEL 


Temperature  -  Ambient 


FIGURE  34 

ACCELERATION  RESPONSE 
FREE-PANEL 


Temperature  -  450°F 


FIGURE  35 

ACCELERATION  RESPONSE 
FIXED  PANEL 

Temperature  -  Ambient 


FIGURE  36 

ACCELERATION  RESPONSE 
FIXED  PANEL 

Temperature  -  480°F 


TABLE  8 

CHANGE  IN  NATURAL  FREQUENCY  OF  PANEL  DUE  TO  THERMAL  EFFECTS 
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9.4 

There  are  several  methods  available  for  controlling  boundary  temperatures.  A  massive 
fixture  with  good  thermal  conductivity  will  act  as  a  heat  sink  which  maintains  ambient 
boundary  temperatures.  This  method  of  controlling  boundary  temperature  is  particularly 
applicable  when  dynamically  fixed  boundaries  are  required  to  be  maintained  at  ambient 
temperature.  Convective  cooling  or  heating,  where  a  gas  or  liquid  is  forced  to  flow  over 
the  boundary,  may  also  be  used.  This  method  requires  a  good  analysis  of  the  convective 
heat  transfer  problem  to  determine  proper  flow  rates,  contact  areas  and  fluid  temperature. 
Care  must  also  be  exercised  in  the  attachment  of  hardware  to  control  the  fluid  flow  since 
the  in-piane  and  bending  stiffness  of  the  boundary  may  be  affected.  A  third  method  of  con¬ 
trolling  boundary  temperature  is  electrical  resistance  heating  of  the  boundary.  In  this  case 
the  specimen  must  be  electrically  isolated  from  the  boundary  but  have  good  thermal  contact. 
Several  materials  with  these  properties  are  available  commercially,  such  as  cupric  oxide  ce¬ 
ments.  In  each  case,  boundary  heating  due  to  radiation  from  the  specimen,  lamp  banks,  and 
reflectors  must  be  accounted  for  or  the  boundary  must  be  insulated  against  radiation  from 
these  sources. 

Once  T(x,  y)  hasjbeen  established  with  the  boundary  controlled  at  a  constant  tempera¬ 
ture  Tg,  the  value  of  Tr  (Figure  37)  may  be  computed.  As  shown  in  section  4.2,  if  small 
gradients  exist  along  the  boundary,  use  of  an  average  boundary  temperature  yields  satisfac¬ 
tory  results.  The  next  problem  is  to  design  a  mounting  such  that  the  in-plane  stiffness  at 
the  boundaries  produces  the  desired  static  stress.  To  examine  the  possible  stresses  which 
may  be  developed  in  the  panel  we  observe  the  stress  limits  shown  in  Figure  37.  As  mention¬ 
ed  in  section  4.2,  when  the  boundary  temperature  is  ambient  (Tq  =  1 .0),  the  center  stress 
is  independent  of  the  in-plane  boundary  loads.  The  value  of  the  center  stress  for  the  ambient 
boundary  temperature  case  is  given  by  -{1  ~-v)/2.  Thus  for  a  given  value  of  v,  a  plot  such  as 
the  one  in  Figure  37  can  be  constructed.  The  points  on  the  Oq  axis  have  urn  range  -1 .0< 

<j(j  <0.  These  points  represent  the  possible  in-plane  boundary  stiffnesses  from  infinite  to 
zero  (fixed  to  free).  As  shown  in  the  figure  for  values  of  oq  >  -  (1  -v)/2  a  gradient 
O’c  >  0)  will  increase  the  magnitude  of  the  center  stress.  For  values  of  Gq  <  -  ( 1  -u)/2  a 
gradient  will  result  in  a  decrease  of  the  center  stress.  As  shown  in  Figures  1 8  through  25 
the  stress  at  points  other  than  the  center  is  a  function  of  the  order  of  T(x,  y).  Hie  curves 
in  these  figures  show  two  significant  characteristics:  the  form  ot  the  stress  distribution  is 
independent  of  the  boundary  condition,  and  the  magnitude  of  the  stress  at  a  point  is 
related  to  the  stress  at  the  same  point  in  the  free-boundary  case  by  an  additive  constant. 

This  constant  depends  on  the  in-plane  stiffness  of  the  boundary. 

Figure  37  can  be  used  to  design  the  mechanical  boundary  to  produce  a  given  thermal 
stress  in  a  panel.  The  value  of  as  a  function  of  for  a  given  boundary  is  determined  by 
a  straight  line  which  lies  between  live  lines  Qq  =--((1  ~v)(2)Tq  and  Gq  -  --  1 .0  +  ((|-H>)/2)  T^ 
Since  each  line  must  pass  through  the  common  point  (o^  =  -  ( 1  -n)/2,  ■=  1 .0)  the  line  for 

a  given  mechanical  boundary  is  completely  determined  by  the  specification  of  5^  at  Y q  =  0. 
The  value  at  Tp  =  0  establishes  the  center  stress  for  the  constant  temperature  case.  This 
stress,  however,  is  the  stress  at  every  point  in  tiie  panel  including  the  boundary.  Thus,  if 
we  are  given  Y^  =  0.5,  a  desired  center  stress  of  -0.48  aBTq/(  1  -v),  and  wish  to  design 
a  boundary  which  will  produce  this  result,  we  proceed  as  follows:  A  line  is  drawn  through 
the  points  ?£  =  0,5,  3 q  -  -0.48,  and  ?£  =  1.0,  vq  =  -  (l-v)/2.  Tim  point  where 


this  line  intersects  the  Tq  -  0  axis  then  found.  In  tins  case  the  intercept  has  the  value  7?/-  - 
-0.6.  Substitution  of  oXK  =  ayy  Boundary  =  -0.6  aHT^/O  -t>)  into  the  generalised  Hooke’; 
law  (equations  4.5)  yields  a  value  of  0.4  a  T(-  for  the  strain  at  the  boundary.  Thus  60%  of 
the  free  boundary  strain,  a Tp,  must  be  suppressed.  To  do  this  a  load  per  unit  length.  F, 
must  be  applied  to  each  boundary  where 
C 

P  =J a  dz  =  .6  «ETC  c/(l  ~v)  '  (437) 

_C 
2 

One  method  of  accomplishing  this  loading  is  to  spring  load  the  panel  edges!  Keeping 
the  spring  at  ambient  temperature  the  spring  constant  may  be  determined  from  the  ’load  and 
deflection.  In  this  case  the  deflections,  u  and  v,  for  zero  panel  rotation  or  tnv'xiatiou.  ?re 
0.2  txT^a.  The  spring  constant  per  unit  length  must  then  be  0.31ic/a(l  tJ).  The  problems 
involved  in  designing  such  springs  for  combined  environment  testing  will  be  discussed  later 
in  this  section,  if  the  temperature  distributions  arc  not  of  the.  form  discussed  in  section  4,2 
it  is  suggested  that  the  computer  analysis  described  in  that  section  be  used  to  check  the 
stress  distribution  tor  the  designed  boundary. 

The  problem  of  designing  a  mechanical  boundary  to  control  the  d>  namic  response  of 
a  panel  is  not  as  straightforward  as  the  static  case,  While  the  boundary  problem  is  one  of 
producing  a  given  bending  stiffness  at  the  boundary,  the  dynamic  response  is  vie  pen  lent  oti 
the  in-plane  thermal  loading,  damping,  ami  coupling  with  the  acoustic  field.  The  eerigner 
has  the  option  of  controlling  the  panel  mounting,  the  acoustic  field  or  both  to  produce 
the  desired  dynamic  stress. 


The  basic  problem  in  direct  boundary  synthesis  (and  analysis)  is  the  determination  of 
stress  as  a  function  of  mode  shape  and  frequency.  Thus  control  of  the  acoustic  Held  is 
generally  the  easier  method  of  producing  a  desired  dynamic  stress  level.  If  the  acoustic 
field  is  io  be  the  controlled  variable,  a  dynamic  boundary  may  be  selected  on  tire  basis  of 
convenience  and  case  of  analysis. 


In  the  present  study,  fined  and  free  boundaries  were  designed.  The  design  of  the  fixtsa 
boundary  was  based  on  the  information  presented  in  Reference  19.  The  tree  boundary'  was 
simulated  by  laying  the  panel  on  thermal  insulating  material  with  low  damping  and  low 
translational  stiffness.  The  panel  was  kept  from  rotating  or  translating  as  a  rigid  body  by 
the  application  of  fibergkr?  tape  to  the  panel  edges.  The  tape  also  prevent  the  thermal 
insulation  from  being  disturbed  by  the  air  flow  from  the  acoustic  source,  in  both  case?,  thi 
experimental  boundaries  resulted  in  frequencies  and  mode  shapes  closely  approximating 
the  analytical  predictions. 


Another  mounting  technique  which  proved  successful  is  Ule -application of  Tiyfeimg  re¬ 
straint  (dynamic  boundaries)  after  the  test  specimen  Is  at  the  desired  temperature.  The 
specimen  in  this  case  is  allowed  to  expand  on  rollers  placed  at  0©  boundary.  After  the 
desired  temperature  is  reached  the  edges  may  be  restrained  by  %hfep.ing  the  bolt  ton  a 
picture  frame  type  mounting.  A  mounting  of  this  type  is  shown  in- Figwre 13$.  Some  con¬ 


trol  over  the  dynamic  response  can  be  achieved  by  using  the  variable  dynamic . ■restraint 


(amount  of  torque  on  picture  frame  bolts/  to  tune  the  dumped.- mdurai  frequerteiet  at  the 


higli  temperature  condition.  These  frequencies  would  be  turvedTb  produce,  the  acoustic 


coupling  required  for  the  given  stress  level. 


Although  the  design  of  specimen  mountings  has  t  een  approached  as  the  superposition 
of  three  independent  boundary  designs  it  is  obvious  that  the  three  boundaries  do  interact. 

In  some  cases  the  interaction  of  the  three  boundary  conditions  can  be  used  to  the  designers 
advantage.  For  example  if  the  boundary  is  to  be  at  ambient  temperature  and  both  the  bend¬ 
ing  and  in-plane  stiffness  is  to  be  large,  one  massive  mounting  fixture  would  satisfy  all  three 
conditions,  in  other  cases  it  is  difficult  to  obtain  certain  combinations  of  the  three  bound¬ 
ary  conditions,  it  is,  for  example,  difficult  to  obtain  a  statically  fixed  and  dynamically  free 
boundary  since  physical  methods  of  in-plane  loading  usually  introduce  some  bending 
stiffness.  Two  design  considerations  chou!d  be  included  In  most  combined  environment 
mountings.  First  the  static  and  dynamic  boundaries  should  either  be  isolated  from  the 
thermal  environment  or  maintained  at  a  constant,  known  temperature.  Second,  where 
interaction  between  the  static  and  dynamic  boundaries  is  a  problem,  priority  should  be  giv¬ 
en  to  proper  static  boundary  simulation  since  the  dynamic  stress  can  often  be  controlled  by 
control  of  the  acoustic  field. 


SECTION  5 


SONIC  FATIGUE  INSTRUMENTATION  AND  MEASUREMENTS 

Simulation  of  the  thennal/acousiic  environment  and  proper  mounting  of  the  specimen 
within  this  environment  are  predominant  factors  in  elevated  temperature  sonic  fatigue  test¬ 
ing.  Measurement  of  the  test  conditions  and  monitoring  the  response  of  the  test  specimen 
are  matters  of  no  less  importance.  Consequently,  one  phase  of  this  program  was  to  study 
methods  to  measure  the  parameters  of  interest  in  combined  environment  tests.  The  study 
consisted  of  determining  the  availability  of  instruments  to  make  the  various  measurements, 
and  indicating  the  advantages  and  limitations  of  different  types  of  devices.  The  following 
parts  of  this  section  discuss  sonic  fatigue  environment  and  response  measurements,  and 
methods  for  detecting  sonic  fatigue  cracks. 

5.1  SONIC  FATIGUE  MEASUREMENTS 

Elevated  temperature  sonic  fatigue  tests  require  a  number  of  measurements.  The  most 
basic  parameter  is  the  acoustic  field  itself.  Of  almost  equal  importance  are  the  temperature 
conditions  and  the  response  of  the  specimen  undergoing  test.  Surveys  of  devices  to  make 
these  measurements  were  conducted  as  part  of  this  study  program.  The  results  of  these 
surveys  are  summarized  below. 

5.1.1  Acoustic  Field  Measurements  ' 

A  number  of  good  transducers  are  available  to  measure  high  intensity  acoustic  noise 
Fields.  The  transducers  operate  on  a  variety  of  principles,  including  piezoelectric,  strain 
gage,  capacitive,  inductive,  magnetostrictive,  and  pressure-sensitive  semi-conductors.  Each 
of  these  devices  has  a  metallic  diaphragm  upon  which  the  acoustic  pressure  impinges,  and 
the  resultant  force  is  converted  into  an  •Tectrical  signal  through  one  of  the  above  principles. 
In  selecting  a  pressure  transducer  for  a  specific  application,  the  user  must  recognize  that 
high  sensitivity  and  wide  frequency  range  will  be  difficult  to  obtain  in  a  single  instrument. 
High  sensitivity  results  from  a  low  stiffness  diaphragm,  but  low  stiffness  also  contributes 
to  a  low  resonant  frequency.  Therefore,  a  trade-off  between  upper  frequency  range  and 
sensitivity  must  be  made  in  all  types  of  pressure  transducers. 

Piezoelectric  and  capacitive  transducers  are  generally  usable  to  temperatures  in  the 
500°F  range  without  external  cooling.  The  change  of  sensitivity  with  temperature  is 
ieliably  predicted  in  this  temperature  range,  so  the  transducers  can  be  used  with 
confidence.  Above  500°F,  external  cooling  of‘fhe  transducer  is  generally  required. 

Some  devices  are  being  developed  with  materials  that  can  sustain  higher  temperatures, 
but  the  use  of  cooled  transducers  is  more  reliable  at  the  present  time. 

Several  methods  of  cooling  pressure  transducers  are  employed.  The  most  common 
method  requires  tiiat  a  coolant  fluid  be  circulated  within  tin;  transducer  to  maintain  the 
sensing  element  at  some  specified  temperature.  Another  scheme  is  to  build  the  transducer 
into  a  heat  sink  which  is  then  cooled.  Some  transducers  arc  mechanically  isolated  from 
mechanical  members  that  experience  high  temperatures,  and  are  cooled  by  a  flow  of 
air  or  some  other  gas  between  the  transducer  and  the  high  temperature  environ,  umt.  When 
a  cooled  transducer  is  used,  special  attention  should  be  given  to  the  noise  generated  by  the 
cooling  system.  Likewise,  the  amplifier  used  with  the  pressure  transducers  must  have  a  noise 
level  that  is  sufficiently  low  to  permit  measurement  of  the  required  acoustic  levels. 


Another  concept  evaluated  during  this  study  was  that  of  coupling  the  acoustic  field  to 
a  pressure  transducer  through  a  small  tube.  In  the  case  of  elevated  temperature  sonic  fatigue 
tests,  this  method  would  permit  the  pressure  transducer  to  be  located  away  from  the  ele¬ 
vated  temperature  environment.  A  reliable  room  temperature  transducer  could  then  be  used 
to  measure  the  high  temperature  acoustic  field.  Some  cooling  of  the  coupling  tube  might  be 
required,  but  this  could  be  easily  accomplished. 

The  limiting  factor  in  the  use  of  a  coupling  tube  is  the  frequency  response  of  the  tube. 
The  behavior  of  the  tube  can  be  analyzed  by  considering  the  pressure  transducer  to  be  an 
infinitely  stiff  termination  of  the  tube,  as  shown  in  Figure  39.  Such  an  arrangement  is 
similar  to  the  closed-end  tube  discussed  in  Section  2.1.  The  pressure  at  any  point  x  in  the 
tube  (with  the  x=£  point  taken  as  the  “input”  to  the  tube)  is  given  by 

p(x)  =  - -g—  cos  kx  (5.1) 

sin  k£ 


Although  equation  (5.1 )  neglects  losses,  it  is  adequate  ,j  show  the  limitations  of  the 
coupling  tube.  The  equation  is  readily  evaluated  at  the  termination  and  input  points,  x=0 
and  x=£  respectively.  Finally,  the  gain  of  me  tuoe  is  found  by  taking  the  ratio  of  the  pres¬ 
sure  at  the  transducer  io  the  pressure  at  the  input,  yielding 


P(o)  _  1 

p(S)  cos  k? 
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(5.2) 


(5.3) 


At  sufficiently  small  values  of  kC,  cos  k?  *»  ! ,  and  the  transducer  pressure  is  approxi¬ 
mately  equal  to  the  input  pressure.  These  conditions  are  only  met,  however,  for  freauencies 
at  which  the  tube  length,  8,  is  much  smaller  than  the  acoustic  wavelength.  The  variation  of 
gain  in  the  tube  is  plotted  in  Figure  39  for  two  representative  tube  lengths.  At  frequencies 
where  kC  equals  odd  multiples  of  90°,  equation  (5.2)  indicates  that  the  ratio  approaches 
infinity.  The  losses  in  an  actual  tube,  however,  will  always  limit  the  pressure  ratio  to  a 
finite  value.  Even  with  some  losses  in  the  tube,  it  is  apparent  that  the  use  of  coupling 
tubes  should  V  limited  to  low  frequencies.  This  restriction  reduces  the  usefulness  of 
coupling  tubes  for  sonic  fatigue  tests,  since  signals  up  to  several  kilohertz  must  ordinarily 
be  measured. 


Calibration  of  pressure  transducers  for  use  on  elevated  temperature  sonic  fatigue  tests 
can  be  accomplished  in  several  ways,  The  reciprocity  technique  yields  the  most  accurate  re¬ 
sults,  because  it  requires  only  the  measurement  of  electrical  quantities  and  distances.  This 
method  can  also  be  used  at  elevated  temperatures,  since  the  variation  in  the  physical  proper¬ 
ties  of  air  can  be  determined  for  elevated  temperatures.  Comparison  calibrations,  using  a 
known  “reference”  transducer,  are  also  suitable  and  widely  used.  A  less  common  calibration 
method  is  the  use  of  a  shock  tube,  in  which  a  pressure  step-function  is  applied  to  the  trans¬ 
ducer.  For  those  transducers  whose  bandwidth  extends  down  to  zero  frequency,  a  static 
pressure  calibration  can  be  used  to  establish  the  transducer’s  sensitivity.  Knowledge  of  the 
transducer’s  resonant  frequency  and  damping  are  then  needed  to  extend  the  frequency 
range  of  the  calibration. 
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Measurement  of  the  acoustic  field  in  elevated  temperature  sonic  fatigue  tests  can  be  ac¬ 
complished  with  currently  available  transducers.  At  temperatures  above  500°F,  transducers 
with  external  cooling  should  be  employed.  Properly  selected  and  controlled  co-Jing  methods 
will  result  in  reliable  acoustic  field  measurements  throughout  the  temperature  range  of  sonic 
fatigue  tests.  As  transducer  manufacturers  continue  to  develop  materials  that  will  sustain 
higher  temperatures,  pressure  transducers  will  become  available  that  can  be  used,  uncooled, 
beyond  the  temperature  limits  of  current  instruments. 

5.1.2  Sonic  Fatigue  Strain  Measurements 

A  knowledge  of  the  state  of  strain  in  a  sonic  fatigue  test  specimen  is  a  basic  factor  in 
evaluating  test  results.  The  measurement  of  such  data  is  difficult  in  the  combined  high 
temperature/acoustic  environment  because  this  environment  imposes  detrimental  elevated 
temperatures  and  dynamic  acceleration  loads  on  the  strain  gages.  To  determine  the  avail¬ 
ability  of  strain  gages  for  sonic  fatigue  applications,  a  survey  was  made  of  strain  gages 
currently  on  the  market.  Significant  findings  regarding  the  properties  and  applications  of 
strain  gages  are  summarized  below. 

The  use  of  strain  gages  at  elevated  temperatures  is  subject  to  errors  from  several 
sources.  As  one  example,  the  resistance  of  the  strain  gage  material  may  change  as  a 
function  of  temperature.  In  addition,  there  is  usually  a  differe;  between  the  thermal 
coefficients  of  expansion  cf  the  strain  gage  material  and  the  man.  !  which  is  being  meas¬ 
ured.  The  elastic  moduli  of  the  test  material  and  the  strain  gage  material  ore  also  subject 
to  change  as  a  function  of  temperature.  To  account  for  errors  from  these  and  other  factors, 
strain  gages  must  be  compensated  for  temperature  variations,  according  to  the  specific 
materials  on  which  the  gage  is  installed.  An  alternative  to  compensation  is  calibration  of 
the  strain  sensing  system  throughout  the  load  and  temperature  range  of  intended  usage, 

Several  methods  of  temperature  compensation  may  Ihv  employed.  A  common  method 
employs  a  second  strain  gage  (dummy  gage)  which  is  exposed  to  the  same  temperature 
environment  as  the  gage  of  interest  (active  gage)  without  applying  the  actual  strain  condi¬ 
tions  to  the  dummy  gage.  The  temperature-dependent  output  of  the  dummy  gage  can  thus 
be  used  to  cancel  out  the  apparent  strain  indicated  by  the  active  gage.  For  sonic  fatigue 
applications,  this  method  is  limited  by  the  difficulties  in  locating  the  dummy  gage  in 
precisely  the  same  temperature  as  the  active  gage,  and  yet  assuring  that  the  dummy  gage 
does  not  experience  any  load.  f  . 

Another  method  of  temperature  compensation  is  the  use  of  a  3-wire  system.  This 
method  requires  that  an  additional  lead  be  attached  to  one  of  the  strain  gage  terminals,  and 
be  routed  out  along  with  the  two  normal  gage  leads.  The  additional  lead  is  used  to  identify 
thermally-induced  changes  in  lead  resistance,  which  might  otherwise  be  mistaken  for  changes 
in  gage  resistance  and  read  out  as  strain.. 

Self-compensating  gages  are  also  available  for  temperatures  up  to  800°F.  These  gages 
are  manufactured  with  known  temperature  characteristics  that  can  be  used  to  correct  the 
indicated  strain  levels.  Corrections  depend  on  an  accurate  knowledge  of  the  temperature  at 
the  gage,  and  the  gage  manufacturer’s  compensation  data  for  the  type  i r  material  being 
tested.  Some  gages  are  supplied  with  the  gage  factor  variation  as  a  function  of  temperature. 
This  information  is  useful  when  the  dynamic  component  of  strain  is  being  measured,  since 
static  drift  or  pparent  strain  can  be  zeroed,  and  the  time-varying  strain  may  be  determined 
by  using  the  correct  gage-factor  value  at  a  given  temperature. 


The  method  of  attachment  of  the  gage  to  the  specimen  must  also  be  considered  in 
selecting  gages  for  sonic  fatigue  tests.  The  most  common  method  is  bonding,  with  the  choice 
of  bonding  agent  depending  primarily  on  the  temperature  range  of  the  gage.  Bonding  materials 
range  from  low  temperature  epoxy  cements  to  high  temperature  ceramic  cements  which  are 
use  at  temperature  of  1200  to  1400°F.  For  even  higher  temperatures,  flame  spraying  is 
used  to  attach  strain  gages  to  the  material  being  measured.  Some  gages  are  encapsulated  in 
metal  sheaths  which  arc  spot  welded  directly  to  the  test  specimen. 

The  selection  of  a  gage  attachment  method  must  also  consider  the  compatibility  of 
the  method  with  the  material  being  tested.  Spot  welding,  for  example  is  likely  to  reduce 
the  fatigue  life  of  some  materials,  and  therefore  may  not  be  suitable  for  sonic  fatigue  test 
applications.  Bonding  materials  have  different  thermal  conductivities  than  the  test  specimen 
materials  and  may  cause  sharp  thermal  gradients  when  the  specimen  is  heated. 

In  sonic  fatigue  testing,  a  strain  gage  attached  to  a  panel  experiences  high  dynamic 
accelerations  as  the  panel  is  tested.  This  motion  is  perpendicular  to  the  strain  gage  surface 
and  lead  wires,  and  tends  to  throw  the  wiring  and  strain  gage  connections  away  from  the 
test  panel.  Published  strain  gage  data  dealing  with  “dynamic"  loads  is  derived  from  repeated 
loading  as  applied  by  fatigue  testing  machines.  A  high  speed  machine  of  this  type  may 
operate  at  10  to  20  Hz,  frequencies  that  are  low  compared  to  acoustic  excitation  frequencies. 
The  resultant  low  accelerations  do  not  present  the  lead  retention  problem  that  is  experienced 
on  sonic  fatigue  tests,  where  it  has  become  necessary  to  bond  or  otherwise  attach  the  loads 
securely  to  the  specimen. 

in  summary,  it  is  concluded  that  strain  measurement  methods  for  high  temperature 
sonic  fatigue  tests  require  much  additional  development.  Problem  areas  include  the  follow¬ 
ing: 

a.  Accuracy  of  gages  at  elevated  temperatures  is  questionable. 

b.  Specimens  undergoing  sonic  fatigue  tests  experience  high  stress  h-vels,  and  re¬ 
quire  gages  with  good  fatigue  life  at  high  strains. 

c.  Attachment  of  gages,  whether  by  bonding,  flame  spraying,  or  welding,  can  cause 
localized  thermal  gradients  in  the  specimen.  The  same  problem  applies  to  the 
attachment  of  load-out  wires. 

Manufacturers  of  strain  gages  are  continually  improving  the  temperature  characteristics 
of  their  products.  However,  as  long  as  it  is  necessary  to  physically  attach  a  device  to  an  ele¬ 
vated  temperature/sonic  fatigue  specimen,  there  will  be  uncertainties  related  to  the  thermo¬ 
dynamic,  chemical,  and  dynamical  effects  of  the  attached  device.  Therefore,  other  means  of 
estimating  the  stress  level  in  a  sonic  fatigue  panel  should  be  pursued.  Concepts  that  have 
potential  for  this  purpose  include  optical  methods,  and  accurate  specimen  deflecrion  meas¬ 
urements  from  which  the  stress  can  be  calculated. 

5.1.3  Displacement  Measurements 

Measurement  of  the  dynamic  displacement  of  a  specimen  undergoing  sonic  fatigue  * 
testing  provides  a  good  indication  of  the  dynamic  properties  of  the  specimen.  The  wide 
range  of  amplitudes  and  broad  frequency  range  of  sonic  fatigue  specimens  place  severe  re¬ 
quirements  on  devices  used  to  make  displacement  measurements.  As  part  of  this  program, 
a  survey  was  made  to  identify  devices  that  might  be  used  to  make  displacement  measure¬ 
ments  on  sonic  .fatigue  specimens  at  elevated  temperatures.  Several  types  of  devices,  along 
with  their  limitations,  are  described  in  the  following  p  "  raplis. 


87 


Devices  to  measure  displacement  were  considered  to  fail  in  two  broad  categones:  de¬ 
vices  which  must  be  attached  to  the  specimen  being  tested,  and  devices  which  do  not  make 
physical  contact  with  the  specimen.  The  former  type  includes  accelerometers,  potentiomet- 
ric  devices,  and  linear  variable  differential  transformers  CLVDT’s).  The  disadvantage  of 
contacting  devices  is  that  they  add  mass  to  the  test  specimen.  In  the  case  of  thin  sonic 
fatigue  test  specimens,  even  light-weight  accelerometers  can  affect  the  response  characteris¬ 
tics  of  the  tc>  i  specimen.  Another  disadvantage  is  that,  in  the  combinedthermal/acoustie 
environment,  devices  attached  to  the  specimen  produce  undesirable  temperature  gradients 
in  the  area  of  the  attachment. 

The  frequency  range  of  accelerometers  is  generally  adequate  for  sonic  fatigue  work. 
Other  types  of  contacting  devices  have  more  limited  frequency  response.  For  example, 
LVDT’s  depend  on  modulation  of  the  signal  applied  to  the  fixed  windings  of  the  device. 
Since  the  excitation  frequencies  are  generally  a  few  kilohertz,  the  frequency  of  the  motion 
being  measured  is  limiie  '  fo  several  hundred  Hertz.  This  type  of  device  is  not  suitable  for 

measuring  motion  at  the  higi.  -r  frequencies  encountered  in  sonic  fatigue  testing, 

» 

Potentiometric  devices  are  limited  because  or  the  mechanical  problems  associated  with 
their  operation.  Coupling  of  the  specimen  motion  to  the  resistive  elements  is  cumbersome, 
and  requires  special  precautions  to  avoid  free-play.  If  the  potentiometric  device  incorporates 
a  wire-wound  resist1' -  element,  resolution  is  limited  by  the  size  of  the  wire.  Continuous  re¬ 
sistive  elements  have  their  resolution  limited  by  the  wiping  element  in  the  device.  All 
potentiometric  devices  are  subject  to  wear  and  noise  problems,  and  are  not  recommended 
for  sonic  fatigue  specimen  displacement  measurements. 


Noncontacting  displacement  measuring  devices  operate  on  one  of  several  principles. 
The  noncontacting  transducers  may  be  subdivided  into  two  general  categories:  electrical 
and  optical  devices. 


The  electrical  devices  depend  on  either  capacitive  or  inductive  measurements  to  deter¬ 
mine  the  motion  of  a  specimen.  In  both  types,  a  probe  is  moved  to  within  a  specified  dis¬ 
tance  of  the  test  specimen,  which  is  referred  to  as  the  “target”.  Then,  motion  of  the  target 
relative  to  the  fixed  probe  is  read  out  us  a  change  in  capacitance  or  inductance,  or  as  a 
change  in  capacitive  or  inductive  coupling  between  elements  of  the  measuring  system. 
Displacement  can  therefore  be  measured  without  physically  contacting  the  moving  sped- 
ment.  The  frequency  range  and  resolution  of  these  systems  are  adequate  for  sonic  fatigue 
test  requirements.  A  basic  problem  with  systems  of  this  type  U  one  of  linearity.  Since  the 
electrical  quantities  being  measured  are  not  truly  linear  functions  of  the  distance  between 
the  probe  and  target,  it  is  necessary  that  the  range  of  operation  be  limited.  The  device  can, 
of  course,  be  calibrated  over  a  larger  range  than  the  linear  range.  High  temperature  opera¬ 
tion  of  these  devices  is  usually  limited  by  the  temperature  characteristics  of  the  probe  unit. 
Since  the  measurements  involve  some  elecMcal  property  of  the  target,  the  resistivity, 
permeability,  and/or  permittivity  of  the  target  material  can  affect  the  accuracy  at  high 
temperature. 


An  increasingly  important  type  of  noncontacting  transducer  is  the  optical  device. 

Some  optical  devices  focus  a  beam  of  light  on  a  prescribed  spot  on  the  target  materials,  and 
track  the  spot  as  the  target  moves.  Other  types  include  laser  interferometers,  which  produce 
alternate  light  and  dark  lines  on  a  detector  as  the  target  specimen  moves  through  small 
distances.  Optical  devices  have  a  common  advantage  with  the  electrical  honcontaeting  de¬ 
vices,  in  that  neither  type  mass-loads  the  specimen  or  affects  the  thermal  properties  in  any 
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way.  The  optical  devices  may  be  located  farther  from  the  material  than  electrical  devices  - 
typically,  several  feet  rather  than  several  inches.  For  sonic  fatigue  tests  at  elevated  tempera¬ 
tures,  however,  optical  methods  presently  have  practical  limitations.  For  one,  the  location 
of  the  test  specimen  in  a  progressive  wave  tube  or  other  enclosure  may  restrict  the  optical 
path  between  the  specimen  and  the  measuring  device.  If  the  test  is  being  conducted  at 
elevated  tempera'  ..re,  the  heat  lamps  and  their  reflectors  tend  to  further  complicate  the 
optical  path  problem.  Another  problem  present  at  elevated  temperature  is  the  radiation 
from  the  test  specimen.  Optical  displacement  measuring  devices  have  their  own  light 
sources,  and  specimen  radiation  would  appear  as  noise  to  the  optical  detector. 

In  summary,  several  types  of  devices  to  measure  the  motion  of  sonic  fatigue  specimens 
are  available.  Those  devices  which  must  be  physically  attached  to  the  specimen  add  mass, 
and  therebv  distort  the  thermal  and  dynamic  properties  of  the  specimen  To  circumvent  ’ 
these  disadvantages,  noncontacting  electrical  and  optical  devices  are  more  desirable.  The 
use  of  noncontacting  transducers  requires  recognition  of  the  physical  limitations  of  the 
devices,  particularly  at  elevated  temperatures.  Thermal  effects  related  to  the  sonic  fatigue 
test  specimen,  the  measuring  device,  and  the  medium  between  the  specimen  and  measur¬ 
ing  device  can  all  cause  errors  in  the  use  of  noncontacting  instruments. 

Additional  development  work  is  needed  on  noncontacting  electrical  and  optical  dis¬ 
placement  measuring  devices.  The  use  of  these  devices  in  the  thermal/acoustic  environment 
should  be  fully  evau/ited.  Noncontacting  devices  are  the  most  promising  instruments  for 
measuring  sonic  fatigue  nocunen  displacements. 

5.1.4  Temperature  Measurement 

Temperature  measurements  in  the  range  of  sonic  fatigue  tests  -  up  to  about  2000°h  - 
area  well  established  procedure  under  static  conditions.  In  the  high  intensity  acoustic 
nnise.  environment,  however,  the  motion  of  the  test  specimen  causes  difficulties  with 
those  temperature  transducers  which  must  make  phy  •'cal  contact  with  the  specimen. 

It  is  fortunate  that  noncontacting  devices  arc  becoming  available  which  promise  to 
make  sonic  fatigue  temperature  measurements  a  more  routine  operation. 

The  most  widely  used  type  ol  tempomture  .transducer  is  tire  acurate,  low  cost,  thermo¬ 
couple.  this  device  operates  on  the  temperature-dependent  change  in  contact  potential 
between  two  dissimilar  metals.  A  number  of  materials  are  commonly  used  for  thermo¬ 
couples  -  copper,  constantan.ehromel,  iron,  platinum,  and  rhodium  are  examples.  'Hie 
proper  combinations  of  materials  am'  alloys  provide  thermocouples  .or  use  at  tempera¬ 
tures  well  in  excess  o'"  3000°F. 


Tile  difficulty  in  the  use  ol  thermocouples  for  sonic  fatigue  applications  results  liom 
the  requirement  that  the  thermocouple  junction  make  direct  contact  with  the  test  specimen. 
Common  practice  is  to  spot-weld  the  thermoc.  uple  junction  in  place.  The  thermocouple 
wires  must  also  be  attached  to  the  specimen,  eidicr  with  cement  or  with  metal  tabs  which 
are  spot  welded  in  piaec.  Hie  thermocouple  must  then  sustain  the  large  dynamic  m.  lions 
that  the  test  specimen  experiences  when  it  responds  to  the  acoustic  excitation.  The  use 
of  thermocouples  on  sonic  fatigue  specimens  thus  becomes  a  question  of  how  well  the 
user  can  physically  attach  the  thermocouple  to  the  test  specimen.  Another  deleterious 
result  of  a  thermocouple,  and  its  associated  attachment  mechanism,  is  that  thermocouples 
rcprtsctU  materials  with  diltcient  thermal  properties  than  the  test  specimen.  Undesirable 
thermal  gradients  in  the  test  specimen  are  a  consequence  of  this  situation. 


Radiation  transducers  largely  circumvent  the  difficulties  inherent  in  thermocouples 
because  radiation-sensitive  devices  do  not  make  ph,  sical  contact  with  the  device  being 
measured.  The  transducer  can  be  located  remote  from  the  specimen  at  a  distance  ranging 
from  several  inches  to  many  feet,  depending  on  the  individual  device. 

Some  radiation-detecting  devices,  generally  known  as  optical  pyrometers,  measure 
the  brightness  temperature  of  a  specimen.  Such  instruments  are  useful  only  for  measuring 
temperature  in  excess  of  about  1G00°F,  since  the  measured  object  must  be  at  a  sufficiently 
high  temperature  that  its  brightness  can  be  measured. 

Other  types  of  radiation  detectors  operate  over  a  wider  range  of  temperatures.  The 
thermal  radiation  from  a  small  area  of  the  test  specimen  can  be  focused  on  a  heat-sensitive 
detector  whose  output  is  calibrated  to  indicate  the  temperature  of  the  radiating  surface. 
Radiation-sensing  devices  can  be  used  with  specimens  whose  temperature  ranges  from  100°F 
to  several  thousand  degrees. 

A  highly  advanced  form  of  radiation  detector  is  the  infrared  scanning  system.  This 
system  continuously  scans  the  surface  being  measured,  detects  the  infrared  radiation  from 
the  surface,  and  presents  a  video  dis»  lay  of  the  thermal  gradient  of  the  surface.  When  the 
temperature  of  a  single  point  on  the  surface  is  known,  the  thermal  gradient  display  is 
calibrated,  and  the  temperature  distribution  is  completely  defined.  The  advantages  of  this 
system  are  immediately  apparent.  First,  the  infrared  scanning  system  does  not  have  any 
physical  contact  with  the  specimen.  Also,  this  single  instrument  provides  complete  tempera¬ 
ture  information  for  the  entire  surface,  eliminating  the  need  to  measure  the  temperature  at 
a  number  of  points,  and  then  estimating  the  temperatuie  between  measured  locations. 

The  use  of  radiation  detectors  in  high  intensity  ucoustr  environments  is  subject  to 
some  restrictions.  The  detectors  must  first  be  able  to  view  the  specimen.  The  detection 
devices  must  also  be  sufficiently  rugged  to  survive  the  acoustic  environment.  Both  of  these 
limitations  are  minimized  by  the  small  size  of  some  devices.  The  infrared  scanning  system 
requires  a  satisfactory  optical  path  between  the  measured  surface  and  the  scanning  camera. 
Since  the  system  operates  at  infrared  wavelengths,  special  optical  components  may  be 
required  to  provide  the  required  viewing  conditions. 

Measurement  of  test  specimen  temperatures  during  sonic  fatigue  tests  can  be 
accomplished  with  either  thermocouples  or  radiation-sensing  devices.  Tin’  familiar  thermo¬ 
couples  have  several  problems  when  used  for  measuring  specimens  undergoing  dynamic  mo¬ 
tions,  as  discussed  above.  Modern  radiation  detectors  have  been  developed  for  temperature 
measurements  without  nuking  any  physical  contact  with  the  specimen.  Further  evaluation 
of  rudiation  detectors  under  actual  sonic  fatigue  test  conditions  is  needed,  to  establish 
special  techniques  for  operation  in  the  acoustic  environment.  Tire  future  for  temperature 
measurement  of  sonic  fatigue  specimens  clearly  lies  iii  the  use  of  noncontacting  radiation 
detectors. 

5.2  SONIC  FATIGUE  CRACK  DETECTION 

Hie  detection  of  cracks  in  sonic  fatigue  test  specimens  presents  a  number  of  challenges 
to  the  experimentalist.  The  complicated  configuration  of  a  representative  aerospace  test 
structure  usually  does  not  suggest  that  any  single  area  is  most  likeiy  to  fail.  Theretorc,  the 
entire  structure  should  be  monitored  for  failure.  The  locations  of  crack-detection  devices 
would  have  to  cover  all  zones  where  failure  is  considered  likr'y  to  occur.  The  addition  of 
elevated  temperature  conditions  further  complicates  the  problem. 
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Several  methods  of  sonic  fatigue  crack  detection  were  evaluated  as  part  of  this  pro¬ 
gram,  Direct  crack  detection  methods,  using  devices  whose  electrical  continuity  was 
interrupted  when  a  crack  occurred,  were  considered  first  because  of  the  success  with  these 
devices  on  static  and  low-frequency  tests. 

One  type  of  crack  detection  material  that  was  evaluated  was  a  conductive  copper  tape 
with  a  dielectric  adhesive  backing.  This  tape  is  available  commercially  in  a  nominal  thickness 
of  0.002  inches,  in  widths  as  small  as  0.01 5  inches.  Use  of  the  tape  as  a  crack-detecting 
device  requires  that  the  tape  be  applied  to  the  specimen  in  an  area  where  a  crack  is  likely  to 
occur.  Tire  tape  is  then  wired  into  an  electrical  circuit,  and  the  continuity  of  this  circuit  is 
monitored.  When  the  specimen  cracks  in  the  designated  area,  the  crack  will  propagate 
through  the  tape,  and  the  resulting  discontinuity  is  detected. 

Evaluation  of  the  conductive  tape  on  a  styiizcd  aluminum  panel  disclosed  that  the 
copper  tape  was  cracking  before  the  test  panel,  giving  false  indications  of  failures.  The  test 
tape  was  made  from  relatively  pure  copper,  which  has  poor  fatigue  characterise  To  be 
useful  as  a  fatigue  detection  device,  tape  with  better  fatigue  properties  would  be  required. 
Material  such  as  beryllium  copper  might  be  suitable.  Present  dielectric  adhesives  are  stable 
to  about  400°F,  and  so  a  different  adhesive  material  would  be  needed  at  higher  temperatures. 

Another  type  of  crack  detecting  material  is  conductive  paint.  This  paint  consists  of 
finely  powdered  silver  suspended  in  a  fast  drying  easier.  The  paint  is  applied  to  the  test  speci¬ 
men  in  narrow  strips  along  areas  most  likely  to  be  cracked  during  sonic  fatigue  tests.  The 
paint  must  be  insulated  from  specimens  which  are  themselves  conductive.  Copper  foil  solder 
tabs  are  located  at  the  ends  of  the  stripes  for  convenience  in  attaching  leads  to  the  conduc¬ 
tive  paint.  Continuity  of  the  conductive  stripe  is  monitored,  and  a  break  in  continuity  is 
interpreted  as  an  indication  of  a  crack  at  some  point  along  the  painted  stripe. 

A  special  panel  was  designed  to  evaluate  this  fatigue  crack  detection  method.  The 
panel,  as  shown  in  Figure  40,  is  clamped  at  each  end  and  subjected  to  vibration  vxctUuon, 
The  notches  in  the  panel  result  in  stress  concentrations  at  the  notches,  and  cracks  are 
initiated  in  the  panel  after  only  a  few  minutes  of  excitation. 

Figure  41  shows  the  display  that  results  when  a  crack  propagates  through  the  con¬ 
ductive  paint.  During  that  part  of  a  vibration  cycle  when  the  painted  surface  is  in  tension, 
the  conductive  path  opens,  and,  when  the  surface  is  in  compression,  the  conductive  path  is 
closed.  It  was  also  noted  that  continuity  was  restored  when  the  vibration  was  stopped, 
indicating  that  the  cracked  edges  of  the  painted  stripe  maintain  good  electrical  continuity 
when  the  panel  is  undellected. 

The  carrier  material  used  in  the  conductive  paint  evaluated  here  is  suitable  for  use  to 
300° F.  Within  this  temperature  range,  these  tests  indicate  that  conductive  paints  are  a  satis¬ 
factory  method  of  detecting  sonic  fatigue  failure  in  specimens  where  the  failure  location 
can  be  predicted.  Development  of  carrier  materials  that  are  stable  at  higher  temperatures 
would  extend  the  use  of  this  method. 

There  are  other  schemes  lor  observing  sonic  fatigue  failures  that  do  not  detect  the 
failure  directly,  but  rather  monitor  some  parameter  which  is,  in  turn,  influenced  by  the 
existence  ol  a  crack.  One  such  parameter  is  the  strain  at  some  point  on  the  panel.  A  strain 
gage  location  shoo'd  be  selected  based  on  the  predominant  normal  mode  shapes  of  the 
panel  under  test.  A  change  in  stress  di.  Tibulion  in  the  panel,  resulting  from  a  fatigue  crack, 
will  then  be  detected  by  that  strain  gage. 


91 


FIGURE  40 

CRACK  DETECTION  PANEL  DESIGN 

(All  Dimensions  in  Inches) 
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An  experiment  to  demonstrate  the  above  effect  was  conducted  as  part  of  this  program. 
A  multipie-bay  panel  was  tested,  with  a  strain  gage  attached  to  the  panel  near  one  of  the 
rows  of  fasteners  on  the  panel.  Since  the  excitation  was  random  acoustic  noise,  the  output 
of  the  strain  gage  was  normalized  with  respect  to  the  output  of  the  microphone  in  the 
acoustic  field.  By  observing  the  ratio  of  the  two  signais,  small  changes  in  strain  gage  output 
resulting  from  small  changes  in  acoustic  excitation  would  not  be  mistaken  for  a  change  in 
specimen  characteristics.  When  an  abrupt  change  in  the  ratio  was  observed,  the  test  was 
stopped  and  the  panel  was  inspected.  A  section  of  the  panel  was  found  to  have  buckled, 
causing  the  reduction  in  the  output  of  the  observed  strain  gage.  While  this  metfi'M  is 
subject  to  the  usual  limitations  of  applying  strain  gages  to  high  temperature  sonic  fatigue 
specimens,  it  has  the  sensitivity  to  detect  panel  changes  caused  by  buckling  or  cracking. 

Another  means  of  detecting  failure-induced  changes  in  the  dynamic  properties  of  a 
panel  is  by  use  of  a  microphone.  In  this  application,  a  panel  undergoing  tests  was  enclosed 
on  the  side  that  was  not  exposed  to  the  acoustic  excitation.  The  enclosure  consisted  of  a 
rigid  walled  box,  lined  with  material  that  had  high  :  oustic  absorption.  A  microphone  was 
placed  inside  the  enclosure,  and  its  output  was  dominated  by  the  fundamental  response 
frequency  of  the  panel.  As  the  panel  edges  cracked,  small  changes  were  observed  in  both 
the  amplitude  and  frequency  of  the  microphone  signal.  The  use  of  a  microphone  does  not 
physically  contact  the  specimen.  High  temperatures  on  the  panel  would  cause  no  problems, 
because  the  temperature  of  the  microphone  could  be  controlled  independent  of  the  speci¬ 
men.  In  addition,  the  microphone  lias  no  effect  on  either  the  thermodynamic  or  structural 
dynamic,  properties  of  the  test  panel. 

A  more  sophisticated  method  ot  monitoring  sonic  fatigue  cracks  would  be  through 
the  use  of  the  infrared  scanning  device  described  in  Section  5.1 .4.  liven  a  small  crack  in  a 
heated  specimen  would  represent  a  sharp  thermal  discontinuity  in  the-  test  specimen-  A 
discontinuity  of  this  type  would  be  ck  ^ly  and  immediately  evident  on  the  video  display 
of  the  infrared  scanning  system. 

Based  on  the  results  of  tho  study,  the  use  of  a  microphone  to  monitor  the  response 
of  a  panel  is  recommended  as  a  reliable,  low  cost  technique  for  detecting  specimen  failure., 
litis  method  is  most  suitable  for  simple,  single-bay  pamT  which  the  response  is  dominated 
by  one  or  two  panel  resonances,  for  more  complex  structures,  the  use  of  conductive  paint 
shows  promise  for  monitoring  several  potential  failure  areas  simultaneously.  Additional 
development  work  is  needed  to  identify  paint  materials  usable  at  temperature*  in  cxeers 
of  300nF. 

Only  one  method  •  the  infrared  scanning  system  •  Iras  the  ability  to  eonUmioudy 
monitor  a  full  structural  surface  for  fatigue  cracks.  The  high  cost  and  optical  path  re¬ 
quirements  of  this  system  may  limit  its  application.  Nevertheless,  the  feature  of  this 
method  w  ill  mein  ns  vonsidcr.it ton  on  tonic  fatigue  tests  where  a  precise  knowledge  of 
fatigue  crack  initiation  o  of  prime  importance. 
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SECTION  6 


SUMMARY  AND  RECOMMENDATIONS 

Tasks  described  in  this  report  covered  four  major  areas  related  to  high  temperature 
sonic  fatigue  tests.  These  areas  are  the  acoustic  environment,  thermal  environment,  effects 
of  specimen  mounting  on  thermal  and  dynamic  response,  and  instrumentation  and  measure¬ 
ments.  The  results  in  each  area  are  briefly  summarized,  and  recommendations  for  further 
study  are  indicated. 

The  studies  of  acoustic  fields  have  shown  that,  in  cases  where  the  acoustic  field  is 
not  diffuse,  simple  trigonometric  functions  can  be  used  to  describe  pressure  distributions. 
Although  these  functions  do  not  epresent  exact  equations  for  the  pressure  field,  they 
are  useful  in  estimating  the  response  of  a  structure  by  the  participation  factor  concept. 

It  was  also  determined  that  the  air  flow  from  air  stream  modulator  noise  sources  is  suffi¬ 
cient  to  preclude  sharp  rises  in  the  air  temperature  during  elevated  temperature  tests. 

The  properties  of  air  are  such  that  under  the  flow  conditions  described  here,  there  is 
low  heat  transfer  to  the  air  from  the  heat  sources  and  the  test  specimen.  Therefore,  the 
usual  ambient  values  of  air  density  an  eed  of  sound  can  be  used. 

Additional  experimental  studies  would  be  desirable  to  show  the  effect  of  acoustic 
field  directional  properties  on  structural  response  As  a  first  tep,  an  existing  progressive 
wave  and  reverberant  test  enclosure  should  be  mapped  to  determine  the  pressure  distribu¬ 
tions  in  these  enclosures.  Then,  a  stylized  panel  could  be  designed  with  a  high  participation 
factor  in  the  p"\gressive  wave  test  section.  The  panel  would  be  tested  in  both  enc’osures, 
to  show  the  difference  in  response  level  for  different  excitation  fields. 

For  thermal  environment  simulation,  a  computer  program  has  been  written  to  calculate 
the  two-dimensional  temperature  distribution  in  a  panel  that  is  heated  by  a  quartz  lamp- 
bank.  Work  has  also  been  reported  on  lampbank  design  considerations  for  a  combined 
thermal/aeerstic  environment.  Methods  for  introducing  damping  into  the  reflectors  and 
for  improving  the  quartz  lamp  supports  were  discussed. 

A  needed  extension  to  the  thermal  environment  studies  is  the  addition  of  temperature 
gradients  along  the  thickness  of  the  test  panel.  This  would  permit  prediction  of  temperature 
in  honeycomb  structures  mid  other  sonic  fatigue  specimens  where  temperature  variations 
normal  to  the  surface  may  be  expected.  Additional  work  on  quartz  lamp  holders  is  also  need¬ 
ed  to  improve  the  reliability  of  quartz  lamps  in  the  high  intensity  acoustic  noise  environ¬ 
ment.  Still  another  area  for  work  would  be  a  detailed  study  of  the  mechanism  by  which 
quart;  mips  deliver  heat  to  a  test  specimen.  This  effort  would  require  an  experimental 
determination  of  the  variation  of  the  heat  transfer  process  at  different  temperatures, 
and  identification  or  al;  thermal  losses  that  occur  in  heating  sonic  fatigue  specimens.  Empir¬ 
ical  data  of  this  kind  are  needed  10  provide  a  realistic  estimate  of  what  fraction  of  the  power 
in  quartz  lamps  is  available  for  heating  test  specimens. 

Considerable  work  has  been  done  on  the  thermal  and  dynamic  response  of  a  test 
specimen  as  influenced  by  the  specimen  mounting.  In  a  large  number  of  cases,  the  thermal 
and  dynamic  problems  can  be  treated  separately  and  then  combined,  according  to  the 
principles  of  superposition.  Computer  programs  were  developed  for  the  plane  stress  response 
to  a  thermal  field.  Basie  equations  were  also  presented  to  show  how  thermal  in-plane  loads 
enter  into  the  dynamic  equations.  Experimental  data  was  presented  to  illustrate  the  thermal 
effect  on  the  dynamic  properties  of  a  panel. 
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A  logical  continuation  of  the  panel  thermal  stress  studies  would  be  to  extend  the  stress 
analysis  to  the  three-dimensional  problem.  Even  in  the  plane  stress  analysis,  however,  other 
length-to- width  ratios,  boundary  conditions,  and  material  properties  could  be  evaluated  to 
supplement  the  temperature-stress  curves  given  in  this  report. 

The  determination  of  stress  in  a  panel  from  a  knowledge  of  the  panel  geometry, 
material,  and  deflection  is  indicated  by  the  dynamic  equations.  Evaluation  of  the  stress 
would  first  require  a  response  test  on  the  panel  to  identify  and  map  the  first  few  normal 
modes  of  the  panel.  Then,  the  panel  should  be  subjected  to  a  specific  thermal/acoustic 
environment  and  the  deflection  should  be  measured  at  two  or  three  points  on  the  panel. 
From  these  deflection  measurements,  the  modes  in  which  the  panel  is  responding  can  be 
determined  and  amplitudes  assigned  to  the  mode  shapes.  The  deflection  of  the  panel  can 
also  be  calculated  from  the  measured  mode  shapes  if  the  coupling  with  the  acoustic  field 
is  known.  Determination  of  bending  stresses  in  the  panel  can  then  be  obtained  by  carry¬ 
ing  out  the  differentiation  operations  in  the  dynamic  equations. 

A  survey  of  devices  to  make  high  temperature/sonic  fatigue  measurements  was  made. 
Currently  available  instruments  to  measure  the  acoustic  pressures  and  specimen  ten  ,;  ..ia- 
tures  appear  to  be  adequate  for  the  ranges  currently  encountered  in  thermal/acoustic 
testing.  Measurement  of  specimen  deflections  and  strains,  however,  are  definite  problem 
areas  in  the  combined  environment.  Because  of  inherent  shortcomings  in  strain-measuring 
devices  that  must  be  physically  attached  to  the  specimen,  indirect  methods  for  measuring 
strain  should  be  pursued,  as  discussed  above,  Deflection  measurements  with  noncontacting 
optical  and  electrical  devices  are  potentially  suitable  for  combined  thermal/acoustic  test¬ 
ing,  although  individual  devices  should  be  thoroughly  evaluated  to  determine  corrections 
that  must  be  added  in  the  thermal/acoustic  environment. 


APPENDIX  A 

STRUCTURAL  COUPLING  AND  ACOUSTIC  FIELD 


A.i  STRUCTURAL  RESPONSE  -  PERIODIC  FORCING  FUNCTION 


The  structural  response  to  periodic  forcin'?  functions  will  be  c  ‘ermined  and  the 
results  will  be  presented  in  terms  of  the  participation  factor  and  the  single-degree-of- 
freedom  system  functions. 

Tire  coupled  dynamic  equations  of  equilibrium  will  be  written  and  the  equations  will 
be  decoupled  by  rotation  into  normal  coordinates.  The  Fourier  Transform  of  the  response 
in  normal  coordinates  will  be  performed  to  show  the  relation  between  the  response,  the 
participation  factor,  and  single-degree-of-freedom  system  functions. 


I, 


Dynamic  Equation: 

[Ml  jwj  +  [Cl  jwj+iR]  jwj  -  jp(x,t)j 
Dynamic  Equation  -  Undamped  Case  Homogeneous  Solutions 

!o! 


[Ml  jwj+t,Kl  jwj  -  joj 
(ID)  -w2  [ID  |w  =  jo; 


(Al) 

(A2) 

(A3) 


Eigenvalue' Solution  and  Modal  Matrix 

o.'  =tonj  i  =  1 , 2,  3, . natural  frequencies 

1 4> )  Modal  matiix  -  each  column  represents  an  eigenvector.  The  matrix  is  also 
normalized  such  that 


(4>]T  l M)  [4*3  =  [11 
[4>1t  [K]  14»1  =  r«nlj 
Dynamic  Equation  -  Dampened  Cases 

14,1  Hr-  hi  ... 

[<frlT[M][$]  «  +1*]1  [C)  [4>1  j5|  +  [4>]rlK|  m  jsj  = 
jp(x,  t)J 

Let  [C]  be  diagonal  or  proportional  to  [Ml  and/or  [K] 

{sJ+^Wn-J  jsj  +  rcon2_J  jsj  =  [4>]Tjp(x,t)J 

a 

pvu 

P0[$]T* 


[*1 


T 


Let  p(x,  t)j  =  P0  |p(x)j  p(t) 

js j  +  r 2?  wn.J  js!  +r-cen2^j  j$ 


P(x)  P(t) 


(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 


Let 


jrj  =  [4>3T  jp(x)j  (A10) 

[0(x,t)j  =p0  jrj  p(t)  (All) 

Participation  Factor  is  F.  The  Participation  factor  depends  on  the  product  of  the  spatial 
distribution  of  pressure,  and  the  structural  mode  shape  of  the  panel. 
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The  Fourier  transform  of  Equation  (A9)  results  in  a  solution  of  the  equation  in  terms 
of  its  system  characteristics  (system. function)  and  the  forcing  function  (Participation 
Factor). 

Let  the  Fourier  transform  of  6(t)  for  the  ith  equation  be2® 


1  rTo12  e  -jnuo.t 
Ai(mF0)=  — -  j  5(t)e  0  dt 

lo  J~To/2 


(AI2) 


The  Fourier  transform  of  (A9)  is 

j  [o^2  -  (mto0)2]  +  j  2£  coni  mco0|  Aj  (mf0)  =  P0Fj  P(mfQ) 
Po  Fi  P(mfo) 


A;(mf0)  = 

m 


-/mcoof  +)2?  (mwo\ 
'wni '  l^ni  / 


(A13) 

(A14) 


where  the  system  function  (Hj(mf0))  is  expressed  as 
Hj  (mf0)  = 


(A1S) 


i  -  ra  +j2r  ps> 


,  CO, 


111 


mw. 


CO 


111/ J 


The  response  is  expressed  in  terms  of  the  participation  factor  and  the  single-degree-of- 
freedom  system  function. 


Aj  (mf0)  =  Hj  (mf0)rj 


PoP(mf0) 


<o, 


m 


(A  16) 


A. 2  STRUCTURAL  RESPONSE  ■  RANDOM  FORCING  FUNCTION. 

The  cross-correlation  of  the  pressure  field  in  normal  coordinates  will  be  expressed 
in  terms  of  the  cross-correlation  of  the  pressure  field  in  the  original  coordinate  system. 
The  cross-power  spectral  density  for  the  pressure  field  in  the  normal  coordinates  will  be 
expressed  in  terms  of  the  cross-power  spectral  density  and  the  modal  matrix.  This  expres¬ 
sion  will  permit  the  transfer  of  any  known  cross-power  spectra  for  the  pressure  field  from 
the  original  coordinate  system  to  normal  coordinates,  where  structural  response  can  be 
determined.  An  expression  which  transfers  the  response  cross-power  spectra  from  normal 
coordinates  to  the  original  coordinate  system  also  is  presented. 

The  cross-correlation  function  is  expressed  as2^ 

Vi(t)‘tIT«  ct/t V>Pj«*>d*'  (A17) 

Rotation  of  the  pressure  p.(t)  into  normal  coordinates  can  be  accomplished  by 
f,  N  j 

Pn(t)=  2  $ni  pj(t)  (A18) 

i=l 
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(A19) 


N  'p 

Pm({)=  2  *nij  P/^ 


j=l 


RPmPm^  2  2  $ni  ^mj  T  2T  /  pi^pj  (t+r)dt 
i=1  i=l  J—T 


N  N  .T  .T  Lim  1_  S 


(A20) 


i-1  j~l 


RPmPm  =  S  4*,Ij  RpiPKO 

1=1  J=1 


(A21) 


The  Fourier  transform  of  the  cross-correlation  is  the  cross-power  spectral  density. 


•>  Spnpm <*>*/  npnPm^^5  ^ 


(A22  )’ 


N  N 


T  T 


Spnpm^=  jfj  *ni  $mj  Spipj(f) 
In  matrix  form  this  equation  is 
l  j1l  =  t<J>lTlSp]  [4*] 


(A23) 


(A24) 


.Equation  (A24)  transforms  the  cross-power  spectral  density  of  a  known  pressure  field  into 
normal  coordinates.  The  cross-correlation  of  response  in  normal  coordinates  is 


Lim  1  r 1 

R6  6  (T)  ”  t  _  tv  I  5n(t)  5m  (t+r)dt 

°n°m  T-^oo  2TJ_T  n  10 


(A25) 


The  transfer  of  a  response  parameter  from  normal  coordinates  to  the  original  coordinate 
system  can  be  accomplished  by 


N 

wi(t)=  2  5n(t) 

r*=l  11 

Therefore,  the  cross-correlation  of  the  response  is 


(A26) 


N  N  i  jjj«  i  -i 

Vi(r)=S,Z  nl  Wn(^r)dt  (A27) 

1  J  n=l  m=l  1  t 


T 

T 


N  N 

( r )  “2  2  <hjm  R§  5m  (f) 

1  J  n=l  m=l  11 


(A28) 


The  Fourier  transform  of  the  cross-correlation  is  the  cross-power  spectral  density. 

■/  N  N 

.^WjWi^X*8  2,  2  ,  ^iti^jm  ^6n5m^^ 


1  J  n®!  m=l 


(A29) 
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In  matrix  form  this  expression  is 


[Sw]  ^  C^l  [S5]  [^] 1  (A30) 

Equation  A3G  transfers  the  structural  response  parameter  from  normal  coordinates  to  the 
original  coordinate  system. 

A.3  POWER  SPECTRA  FOR  TYPICAL  ACOUSTIC  ENCLOSURES 

Normal  acoustic  test  enclosures  have  progressive  wave  or  reverberant  fields.  Both  of 
these  cases  can  normally  be  considered  single  input  multiple  output  systems.  The  acoustic 
enclosure,  in  a  large  majority  of  cases,  is  excited  by  a  single  or  a  group  of  sources  which 
have  the  same  pressure  amplitude  and  phase.  For  this  type  of  system  the  source  can  be 
considered  single.  The  pressures  at  the  structural  test  panel  will  have  a  variation  in 
amplitude  and  phase  which  is  distributed  spatially  across  the  panel.  The  pressures  at  these 
points  (pDj)  for  this  study  will  be  considered  to  be  multiple  outputs  from  the  single  pressure 
source  (pjj). 

We  will  investigate  first  multiple  input-output  systems.  The  multiple  output  of  a  multi¬ 
ple  input  system  can  be  described  as  a  function  of  the  conjugate  system  function,  input 
cross-power  spectral  density  and  the  system  function.  In  matrix  notation  +he  expression  Is 


IS.  }  =[H*1  IS-]  (HI  1 

po  Pi 

The  cross-power  spectral  density  between  the  ith  and  jth  location  is^ 

N  N 

Sp  .p  .(f)-  2.  2,  Hin(OHjm(f)S 
PoiPqj  n=l  m=l  111  pinpim 


(A31) 


(A3  2) 


If  we  assume  that  we  have  a  single  input  system,  that  is,  the  enclosure  is  excited  by  a  single 
source  or  a  group  of  sources  which  can  be  considered  as  a  single  source,  the  cross-power  i 
spectrum  for  the  pressure  between  the  ith  and  jth  location  is 

sPip/0=  H*u  Hji  CAM)  ' 

,  ‘  I 

Let  the  normalize  \  cross-power  spectrum  be  . 

Cpj  Pj  (f)  *  H*U  Hjj  (0  (A34) 

Equation  (A34)  gives  the  normalized  cross-power  spectral  density  in  terms  of  system 
functions  which  are  characteristics  of  the  acoustic  enclosure.  If  these  system  functions  for 
various  types  of  acoustic  fields  can  be  determined,  then  forcing  functions  in  normal  coordi¬ 
nates  can  be  estimated  by  Equation  (A24). 

A.4  SYSTEM  FUNCTION  -  ANECHOIC  TERMINATION 

In  Section  A.3  we  have  shown  that  the  cross-power  spectrum  for  a  pressure  field  is  a 
function  of  the  system  functions  for  the  enclosure.  In  this  section  the  system  functions  for 
an  enclosure  with  an  anechoic  termination  will  be  developed.  A  system  function  by  defini¬ 
tion  is  the  Fourier  transform  of  the  unit  impulse  response.  The  response  of  a  one-dimension¬ 
al  pressure  field  with  anechoic  termination  to  a  unit  impulse  is 


p(x,t)«  P0fi(ct-x) 


(A35) 
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Ins  Fourier  transform  of  Equation  (A35)  is 


/co  t 

p(x,  t)e  JW  dt 

OO 


COX 


*(? 


'2nx 
P 


(A36) 

(A37) 

(A38) 

(A3  9) 


P(x,f)  =  P0e 

'  ^  (2ts\\  .  /2tcx\ 

_ "  P(x, f)  =  p0 (cos  -)sm(^-) 

The  normalized  system  function  (H(x,  f))  is 

;  m  /2irx\  (2irx\ 

H»,f)  =  coS^-} 

The  system  function  of  Equation  (A39)  describes  the  pressure  field  as  a  function  of 
frequency  and  distance  from  the  source.  Equation  (A39)  is  Equation  2.2,  the  pressure 
distribution  in  a  tube  with  anechoic  termination. 

A.5  RESPONSE  TO  RANDOM  ACOUSTIC  FIELD-ANECHOIC  THERMINATION 

The  normalized  cross-power  spectra  for  the  ith  md  jth  location  in  an  acoustic  field  with 
anechoic  termination  is 

!  (f)Hj,  (f) 


W * (cos (v)+jsin(?))  (£0S © ”jsi"  ("?)) 

Transforming  the  above  equation  into  normal  coordinates  gives 


(A40) 

(A41) 


C0S,'X 


(A42) 


’  Equation  (A42)  is  the  normalized  cross-power  spectrum  in  normal  coordinates  and  it  is 
separable  In  i  and  J/  Equation  (A42)  can  be  rewritten  in  terms  of  the  participation  factor. 

cWf,r'r«r">  (M3) 

The  normalized  power  spectrum  for  the  case  n  equals  m  in  normal  coordinates  is 


eg  g  (0*  in 

vr. 


2 


(A44) 


Equations  ( A43 )  and  ( A44 )  indicate  that  the  cross-power  spectrum  in  normal  coordi¬ 
nates  for  onessure  tMd  with  an  anechoic  termination  is  a  function  of  the  participation 
factors  which  were  developed  for  periodic  functions.  The  normalized  cross-power  spectra 
for  the  pressure  field  in  matrix  form  is 

(A4S) 
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The  response  of  structures  in  normal  coordinates  is 

[S6]  =  r~Htj  [  V*  rm3  r-H-j  Sp_p,  (f)  (A46) 

The  system  function  H  in  Equation  (A46)  is  the  single  degree-of-freedom-system 
function  in  normal  coordinates  for  the  structure.  For  structural  damping,  which  is  normal^ 
ly  less  than  5  percent,  the  following  approximation  can  be  written  at  the  normal  mode 
frequencies. 


lHu(f)  P»H*iHjj 

tS6l  »tx|[12'|H|2s|  S  (0 

11 


(A47) 


(A48) 


The  above  expression  is  a  diagonal  matrix  and  is  approximately  true  for  cases  when  the 
modal  frequencies  are  separated  by  at  least  1 0  percent. 

Equation  (A48)  indicates  that  the  response  to  a  random  acoustic  field  in  an  enclosure 
with  an  anechoic  termination  is  a  function  of  the  participation  factors,  Tj.  The  response  to  a 
random  acoustic  field  is  similar  to  the  response  to  periodic  functions  since  the  forcing  func¬ 
tion  in  normal  coordinates  is  the  participation  factor.  The  effect  of  participation  factor 
on  response  wilt  therefore  be  the  same  for  bc*h  random  and  periodic  acoustic  pressure. 

A.6  RESPONSE  OF  STRUCTURES  TO  RANDOM  ACOUSTIC  FIELD-ABRUPT 
TERMINATION  . 

The  response  ~  •>  structure  to  a  random  acoustic  field  will  depend  on  the  cross-power 
spectra  of  the  pressure  within  the  enclosure.  Equation  (A45)  indicated  that  the  normalized 
cross-po  ver  is  a  function  of  the  system  function  for  the  enclosure.  The  system  functions  for 
an  enclosure  with  a  particular  type  of  an  abrupt  termination  will  be  determined.  The  results 
of  this  analysis  will  be  extended  to  the  other  types  of  terminations  and  to  fields  with  three- 
dimensional  properties.  The  system  functionjis  the  Fourier  transform  of  the  unit  impulse 
response.  The  termination  to  be  examined  will  be  an  open  tube.  The  unit  impulse  response 
will  be  found  by  applying  boundary  conditions  to  the  wave  equation.  The  wave  equation 
for  particle  displacement  in  terms  of  unit  impulses  is 


0(x 


0(x.  t)=  0+  6  (ct-x)  +  0„  §(ct+x) 

The  Fourier  transform  ■  :  Equation  (A49)  is 

-jw't 

j3(x,t)e  dt 

’—00 

_iwx 

0(x,o*0+c  V  +fL  o  c 

The  boundary  conditions  for  the  open  tube  are 
x  =  0,  jj  (o,  t)  =  0O  5(t) 

x«  L,  p(L,  t)  =  0 

The  relation  between  particle  displacement  and  pressure  is 


■p(x,t)»  ~P0c2  0(x,t) 


P(x,f)»-jp0c 


.2/" 


-cT“?+ e 
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-j  w> 


(A49) 

(A50) 

(A51) 

(A52) 

(A53) 

(AS4) 

{A55) 


Applying  boundary  conditions  of  (A52)  and  (A53) 

* 

/30=/3+  +  p_ 

j  coL  -j  toL 

O  =  0_.e  'r_  -/3+e  x“ 


0+  = 


P-  e 


i  coL 

r~c~ 


-j  wL 


(A56) 

(A57) 

(A58) 


0-  =  0o 


-j  to  I 
a  tr 


j  wL  -  j  coL 

e  +  e 


p_.  £°f 


-jwL 

T“ 


2co^t) 


jtoL 


0oe 


“C~ 


P+  =  ^ 


2  cos 

The  particle  displacement  for  an  open  tube  is 

/to 


P(x,f)  = 


Po  lcos\c 


(L-x)l 


cos 


/wL\ 


Hie  pressure  is 


C™\T) 


( toL' 


P(x,  f)  “  0opo  ete  tan  \~~~J  cos 


sin 


(A59) 


(A60) 


(A61) 


(A62) 


(A63) 


(A64) 


Hie  normal  mode  frequency  of  the  tube,  to|V,  is  g,v~n  by  ton  -  2jt(iic/4L).  Then,  tan  (wnL/c)s 
tan  (mr/2)  -  although  the  small  losses  present  i;;  the  tube  will  limit  the  pressure  amplitude 
to  a  finite  value.  The  pressure  distribution  in  the  tube  can  then  be  approximated  by 


P(x,f)  =  P0cos 


<A6S) 
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Equation  (A65)  is  the  system  function  for  an  open  tube  and  corresponds  to.  Equa¬ 
tion  2.3  in  Section  2.1 .  The  same  procedure  used  in  Section  A.5  will  result  in  a  similar  ex¬ 
pression  for  the  response.  That  is.  the  response  will  be  a  function  of  the  single-degree-of- 
freedom  system  functions  of  the  structure  and  the  participation  factor,  Tj. 

For  a  three-dimensiona  structure  the  system  function  for  a  large  room  is 


H  (x.  y,  z,  f)  =  [cos 


I2**) 

1  rrvcl 

cos/ 

V 

u  ( 

7 


i 


(A66) 


Sections  A. 2  through  A.6  indicate  that  the  acoustic-structural  coupling  for  random 
functions  can  be  estimated  in  a  manner  similar  to  the  acoustic  structural  coupling  for 
periodic  functions.  The  response  will  be  a  function  of  the  structural  single-degree-of-freedom 
system  functions  and  the  participation  factors. 
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APPENDIX  B 

TEMPERATURE  PREDICTION  COMPUTER  PROGRAM 

The  program  described  here  is  used  to  solve  for  the  temperature  distribution  T-  ■  as 
defined  by  Equation  (3.16)  in  Section  3.1.  This  example  is  for  a  stainless  steel  panel0 

ab°Ut  3  centerlme>  1110  half'Panel  has  been  divided  into  5(  segments  (10x5) 
i  this  division  appeanng  as  1 1  x6  in  the  dimension  statements.  The  1 1  x6  dimension 

accounts  for  the  temperature  heme  calculated  at  corners  rather  than  centers  of  the 
segments. 

.  ™!lqUantituy  f is  the  original  assumed  temperature,  from  which  T  is  calculated  B 
s  the  difference  between  A  and  T  for  each  iteration.  Initially,  10  iterations  are  program- 
The  user  may  request  additional  groups  of  1 0  iterations,  as  indicated  in  line  440. 

t  pr°^m.’  the  c°nstant  R  (line  60)  is  the  same  as  K  in  Equation  (3.14)  Like¬ 

wise,  Rj  (line  130)  is  the  same  as  Kj  in  Equation  (3.1 5), 

The  constant  S  (line  70)  is  initially  assigned  a  value  of  zero.  This  terra  is  later  com¬ 
puted  as  the  value  of  TR  (Equation  (3.9)). 

puter  wSt’em^mfrrdTr116'6/8  f°r  USe  °n  a  conversational,  time-sharing  com¬ 

puter  system.  Therefore,  different  computers  will  require  different  input  (output)  formats. 


PLATEM 


10  DIMENSION  ACJl*6)*TCll*6>*BCli,6> 

20  DATA  B/66*0 •/ 

30  DATA  T/66*0 • / 

40  DATA  A/!2*80.*V*400.*SO.,80.,9*400.*80.*80.*9*400.*80.*80. 
5049*400 *>80 **80 .>9*400. *80./ 

60  READ*  R 
70  S*0. 

80  X*0. 

90  D0  1  1*1*1  1 

100  D0  l  J«l,6  . 

110  AC I*  0)®A( I*0)+460 . 

120  1  CONTINUE 
130  Rl«8.8E-3 
140  R2=Rl/300. 

150  25  D0  2  1*2*10 
160  D0  3  0*2*5 

170  TC1*U>*<AC1,U-1>+A<I-1*U>+ACI*U+1  >+Atl*l* J)*R)/4. 

180  Ta*0>*C4.*T(I*  J)+100.*R2*S*R!*CAU#  J)/100.  >**4> 

190  3  TCI*U)-T(1* J)/C4t*Cl .♦R2*CA(I* J)/10C.)**3>) 

200  0*6 

210  TCI* 0>*c2.*AU,U-»  >+A(i~i*0>+Am>* 0>>R)/4. 

220  T ( I* J)*(4.*T < I*0>  +  1 00.*R2*S+R1 *<A< I* J J/l 00. )**4) 

230  2  TCI*0)*T<I*0)/<4«*Cl . +R2* (AC  I  * J ) / 1 00 . )*#3) > 

240  D0  4  1*2*10 
250  D0  4  0*2*6 

260  4  Ba*0)*<T<I*0)-A(I*0))/4. 

270  D0  8  1*2*10 
280  D0  9  0*2*5 

290  9  TCI*0>*T<  I*0>*B(I-|*U>+B{UI*0>+BCI*U-1  >*B<I*0+1  ) 

300  0*6 

310  8  T<I*0)-TU#0>*Ba-t*0)<B<I*l,0)*2.*B<l*J-l) 

320  1FCX.LT.I0.>  00  TO  63 
330  2*0. 

340  DO  5  1*2*10 

350  D0  5  0*2*6  ♦ 

360  IFCA8S  (BCl*J)).L£.2>  GO  T0  5 
370  2*ABS(Bd*0>  )  . 

380  5  CONTINUE 
390  81*4. *2 
400  PRINT  100*21 

410  100  FORMAT {"THE  MAXIMUM  DIFFERENCE  IS"*F7.3>"  DEG.  F.  M) 
420  PRINT, 

430  PRINT* 

440  PRINT, "IF  REITERATION  DESIRED*  TYPE  11  0  T0  PRINT  OUT" 

450  READ*  22 

46C  IFC22.EQ.0.)  GO  TO  57 
470  G0  te  62 
480  57  CONTINUE 
490  PRINT* 

500  PRINT. 


PLATEM  CONTINUED 


510  D0  6  1=2, 10 

520  D0  6  J=2,6 

530  6  T(I,J)*T(I,J)-460. 

540  PRINT20Q,<TCI.»2>»TCI,3),T(I,4>.>T(I,5>>T(1,6>,I »2,  10) 
550  200  F0RMAT<5F14.3> 

560  G0  T0  73 
570  62  X=0. 

580  63  CONTINUE 
590  X*X+S . 

600  S«0. 

610  D0  7  1*2, 10 
620  D0  7  J*2, 6 
630  S*S+(TCI«  J')/1'00* )*#4 
640  7  A(LJ)»T(I,J) 

650  S*S/45. 

660  60  T0  25 
670  73  CONTINUE 
680  END 
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APPENDIX  C 


DERIVATION  OF  GOVERNING  EQUATIONS 


The  complete  solution  to  the  problem  of  specimen  response  to  a  combined  thermal 
and  dynamic  loading  requires  the  evaluation  of  eleven  parameters  (in  Cartesian  coordinates). 
These  parameters  are  T,  p,  Vx,  Vy,  V2,  oxx,  ay y,  o22,  oX2  and  oyz  (due  to  symmetry  of  the 
stress  tensor  oxy  =  uyx>  etc.)  2? 

To  determine  the  above  parameters  eleven  equations  must  be  solved.  The  solution  of 
these  equations  with  the  associated  boundary  conditions  is  usually  quit.:  difficult  and  often 
can  only  be  accomplished  by  approximate  numerical  methods.  Fortunately  in  many  engineer¬ 
ing  problems  tire  number  of  equations  and  the  complexity  of  th„  equations  may  be  reduced 
by  the  use  of  certain  simplifying  assumptions.  Tire  following  are  the  eleven  basic  equations 
and  the  simplifying  assumptions  which  can  often  be  applied. 

The  first  of  the  eleven  equations  is  the  equation  of  continuity.  This  equation  expresses 
the  conservation  of  mass  and  can  be  written 


P  +  3X  (p  Vx)  +  3y  (p  Vy)  -  d2  (>  V,)  =  0 


(Cl) 


Under  the  assumption  of  a  constant  and  known  density  we  can  reduce  the  number  of  requir- 
e  1  equations  and  unknowns  by  one. 

The  next  three  equations  can  be  derived  from  conservation  of  momentum,  i  hesw 
equations  state  that  the  acceleration  of  a  body  is  equal  to  the  sum  of  the  forces  and  the 
specific  surface  forces.  This  can  be  written  as: 


p  (d(  Vs  ♦  vx  vx  <*  Vy  ds  Vx  +  V?  Vx)s 

p  l*s  +  uxx  +  dy  oXy  *  u»  Ux, 

p  by  +  dj.  -O^y  ♦  dy  Uyy  *  0^  Oy£ 

P  ldt  v£  *  vx  V,  4-  V,.  a.  v,  +  V,  K  V,)  - 


y  vi  't  tvz  1 


e  **2  *  h  Vva*  °u 


(C2a) 


(C2b) 


(C2c) 


lire  fifth  equation  to  be  dealt  with  is  the  energy  equation.  This  equation  states  that 
tire  fate  oi  change  of  kinetic  energy  plus  tire  rate  of  change  of  internal  energy  must  equal 
the  tate  of  noroeelunical  energy  loss  plus  the  (rower  of  the  surface  and  body  forces  acting 
on  the  continuum. 


(C3) 


[3tP  +  3xpVx  +  3ypVy  +  3  zpVz]  (1/2VXVX  +  Zz  Vy  Vy  +  >/2VzVz  +  e)  + 

[d*e  +  Vx3xe  +  Vy3ye  +  Vz3ze]  + 

|yxPt3tVx  +  VX3XVX  +  Vy3yVx  +  v  „3  ,VX]  + 

V  o!  3*V  +V3V  +V3V  +V3  v  u 
vyPi°t  y  x  x  y  yy  y  z  z  yJ 

Vzp[3tVz  +  VX3X  Vz  +  Vy3yVz  +  VZ3ZVZ]  j  - 
jVx[pFx  +  3X<?XX  +  3yoXy  +  3zaxz]  + 

VytpFy  +  3XCJXy  +  3yOyy  +  3^]  + 

Vz['oFz  +  3x0xz  +  3y0yz  +  3z0zzi|  r 

axx^x^x  +  axy^x^y  +  °xz^x^z  +  °xy**y^x  +  °yy^y'vy  + 
oyzdyVz  +  oxz3zVx  +  ayz3zVy  +  azz3zVz 

_  ^x^x  “  ‘V'fy  “  ^z^z 

From  Equation  (Cl),  the  first  term  in  brackets  on  the  left  side  of  (C3)  is  zero.  From 
Equations  (C2)  it  can  be  shown  that  the  expressions  in  braces  on  the  left  and  right  side  of 
(C3)  cancel  term  f o>  term.  By  definition  3jVj  =  t'y.  Equation  (C3)  can  now  be  written  as 

p( 3^e  +  Vx3xe  +  Vy3ye  +  Vz3ze)  - 
°xx^xx  +  °yy^yy  +  °zz^zz  +  “(°xyGxy  + 

°xzGxz  +  °yzGyz^  “  ^xqx  ~  ^y^y  ~  ^z%  ^ 

The  last  six  of  the  eleven  required  equations  are  the  constitutive  equations,  better  known 
in  this  case  as  the  generalized  Hooke’s  law. 

E 


XX  - X  i-Jv)(  1-2 !>)  K  '-^exx  +  ^eyy  +  ezz>  ~ 

(C5a) 

E 

fyy  =  ( i+t^X l-2r>)  f ( *“^fcyy  +  ^exx  +  czz^  ~  (  1  +i^)cvT' ] 

(C5b) 

E 

rzz  =  ( 1  +!>)(  1  -  2v)  l 1 r !  -^ezz  +  ^XX  +  eyy )  -  (  1 

(C5c) 

xy  ~  ^  ^xy 

(C5d) 

xz  ~  ^  7xz 

(C5c) 

ryz  =  C  Tyz 

(C50 

The  term  in  parentheses  on  the  left  of  Equation  (C4)  is  the  total  time  derivative  of  (he 
nonmechanical  internal  energy,  e.  This  energy  is  JcpT  where  cp  is  the  specific  heat  at  con¬ 
stant  deforestation.  We  a"piy  Equation  (CS)  to  replace  o  in  Equation  (C4)  in  terms  of 
strain  and  note  that  if  we  consider  zero  input  forces,  7xy>  TXZ)  and  7yz  are  zero  since  they 
are  not  a  function  of  temperature.  Further,  we  assume  a  purely  conductive  heat  transfer 
mechanism.  Equation  (C4)  then  becomes: 
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(C6) 


E 

PcFt  *  <exx  +  eyy  +  €zz>  + 

kOxxT  +  ayyT  +  azzT> 

The  term  involving  T  represents  the  rate  of  change  of  temperature  due  to  the  motion  of 
the  body.  If  this  temperature  rise  is  small  (as  it  is  in  most  cases)  then  T  **  0.  Similarly  the 
terms  involving  represent  the  mechanical  vibrations  due  to  thermal  changes  and  can  in 
general  be  neglected.  24  The  energy  equation  in  this  case  degenerates  to  the  well  known 
Laplace  Equation 25 

3xxT  +  dyyT+azzT---°  (C7) 

Equation  (C7)  shows  that  the  temperature  can  be  calculated  without  mechanical  con¬ 
siderations.  This  represents  the  quasi-static  thermoelastic  theory.  Before  turning  to  the  cal¬ 
culations  of  the  veiocity  components  from,  the  equations  of  motion  and  the  stress  compon¬ 
ents  from  Equat;ons  (C5),  we  will  examine  the  calculation  of  the  temperature  fleld  in  more 
detail. 

Equation  (C7)  is  the  general  expression  for  the  three  dimensional  steady  state  temp¬ 
erature  distribution  with  zero  heat  sources  and  losses.  In  the  current  discussion  the  specimen 
will  be  a  thin  structure  (9ZT  «  0)  that  has  reached  steady  state.  It  is  assumed  that  the  heat 
is  applied  from  a  quartz  lamp  reflector  system  in  such  a  way  that  we  can  assume  uniform 
heat  generation  within  the  specimen.  The  absorbed  radiation  is  of  such  a  magnitude  that  all 
convective  and  radiative  heat  losses  are  accounted  for  and  the  required  temperature  distribu¬ 
tion  is  achieved.  Equation  (C7)  is  then  written  as 

k.(dxxT  +  ayyT)  =  q  (C8a) 

and 

A  =  ^radiation  loss  +  ^convective  “  9  absorbed  (C8b) 

i 

A  description  of  the  computer  program  used  to  calculate  the  temperature  distribution 
is  presented  in  Appendix  B.  For  the  remainder  of  this  discussion  it  will  be  assumed  that  the 
temperature  distribution  is  known. 

The  analysis  of  the  thermal-dynamic  problem  now  requires  only  the  solution  of  the 
equations  of  motion  and  the  constitutive  equations.  This  is  still  not  an  easy  task,  but 
through  the  use  of  certain  assumptions  further  simplifications  may  be  made.  We  first  make 
use  of  the  theory  of  superposition  and  assume  that  all  deflections  are  small.  Two  problems 
are  then  treated  separately,  one  thermal  problem  and  one  dynamic  problem.  The  compon¬ 
ents  of  the  static  and  dynamic  stress  tensors  add  linearly.  We  will  consider  the  thermal 
problem  first. 

The  thermal  stress  problem  is  a  static  problem.  The  left-hand  sides  of  Equations  (C2) 
represent  the  total  time  derivatives  of  the  velocity  components  which  for  the  static  case  are 
zero.  In  this  analysis,  since  T(z)  is  constant  and  we  deal  with  thin  panel  structures,  the  case 
of  plane  stress  may  be  assumed.  Plane  stress  implies  that  azz,  oxz  and  ayz  are  zero.  From 
Equation  (C2c)  we  note  Fz  must  also  be  zero.  The  equations  of  motion  then  become 

axCTxx  +  ayffxy  =  “  '  (C9a) 


3x^xy  ay°yy 
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(C9b) 


The  body  forces,  Fx  and  Fy,  will  be  neglected  in  the  static  equilibrium: Equations  (C9)> 
For  the  plane  stress  case,  the  constitutive  equations  become 


°xx  ^exx  +  ^eyy  ~0+y)aTJ 

■  V  .  tblOa) 

E  - 

ayy  =  t  v2  feyy  +  wxx  ~  -  ":..v  .. 

<CK*b)  : 

<7Xy  ~  G  t'Xy 

(cioc) 

To  obtain  the  equilibrium  equations  ihTerms  of  displacements  we  substitute  Equations 
(CIO)  into  Equations  (C9)  and  make  use  of  the  following  strain-displacement  relation  for  . 
small  displacements: 

exx“9xU; 

(Cl la)  ; 

eyy=ayv 

(Cl  lb) 

exy  =  9xv+9yu  .. 

Equations  (C9)  then  become:  • 

(Cl  1c) 

j_„2  ~  W  +  udyv  3  +  2(  ^  9y*9yu  +  V  1  =  9xT 

(Cl  2a) 

"  3yl  V  +  +  ^  V  V  +  V  1  =  j§8yT 

(Cl  2b) 

In  order  to  eliminate  the  need  for  the  solution  of  two  simultaneous  partial  differential 
equations  we  introduce  a  function  \p  such  that 

3xi p  =  u 

(Cl  3a) 

Sy^“V 

(Cl  5b) 

This  function,  \p,  is  known  as  the  strain  potential  function 7-^  Substituting  Equations 
(Cl 3)  into  Equations  (Cl  2)  we  have: 

V^x^^yy^  +(1-^)  3y3xy^  = 
dx(9xx^  +  3vy^)  “(!+*')  «3XT 

3y[9yy^  +  ^3XXl/']  +  (l-*0  ~ 

3y[9XX^  +  3yy^]  =  (I +*»)<*  3yT 

(C 1 4a)  ‘ 

Collecting  like  terms, 1  ’  tuations  (C 1 4)  become 

ax^xxl^  +  9yy^"^I+^Q:T^ 

(Cl  5a) 

3y{3xx^  +  ayy^-(i+^)aT]  =0 

(Cl  5  b) 
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Integrating  Equations  (CIS),  we  have 

•  '  3xx^-3yy^-(l+,)aT=fi(y)  +  Ci  (Cl 6a) 

- ;V  hx*  +  3yy  V  -  0  +»)  aT  =  f2(k)  +  C2  (Cl  6b) 

•.  From  Equations  (Cl 6)  we  can  see  that  f(x)  =  f(y)  «  0  and  the  basic  equations  become: 
n._.  +  C  (C17) 

•The  constant  C  is  evaluated  from  a  boundary  condition  on  9XX$  +  Ky'P  ~  ( 1  +n)n T  or 
its  derivatives-..  Thus  the  solution  of  the  thermal  stress  problem  rests  with  the  solution  of 
Equation  (Cl  7)  for  the  one  variable;  ip.  -  'Mv 

The  solution  of  the  dynamic  problem  can  also  be  simplified  so  the  complete' solution 
to  Equations  (C2)  through  (C S)  are  not  . required.  The  left-hand  sides,  of  Equations  (C2)  are 
the  components  of  acceleration  with  respect  to  the  coordinate  axes.  It  is  assumed  that  the 
components  of  acceleration  in  the  x  and  y  directions  can  be  neglected,  and  that  lateral  de¬ 
flections  are  small  compared  with  the  thickness  of  the  panel.  In  the  present  case,  it  is  assumed 
that  the  stress  in  the  z  direction,  o/z,  may  be  neglected  which  implies  that  all  loads  in  the  z. 

,  direction  arc  applied  to  the  midplane.  Equations  (C2)  then  become,  for  zero  body  forces, 


3x°xx  +  3y°xy  +  °z°xz 


(Cl  8a) 

0  =  dX«Xy  +  dyOyy  +.  azOy2  (C  l^) 

paz  =.  3xoxz  +  dydyi  (Cl  8c) 

From  Equations  (C 1 8)  the  equation  of  plate  bending  under  normal  dynamic  loading  is 


WD^  +  2axxy/Dw)  +  3y>yy(Dw)  *  cfh 


(Cl 9)  ' 


The  derivation  of  this  equation  is  not  presented  here  but  may  be  found  in  several  refer¬ 
ences  27,28  In  addition  to  the  lateral  loading,  a  plate  subjected  to  thermal  loading  develops 
thermal  stresses  in  the  plane  of  the  plate  which  must  be  taken  into  account.  Thermal 
stresses. result  in  traction  forces  at  the  boundary  which  produce  moments  about  the  de¬ 
flected  surface.  The  resulting  modification  to  Equation  (Cl 9),  again  derived  m  the  litera¬ 
ture,  for  a  constant  plate  stiffness,  D,  is: 

D  jaxxxxw  +  23xxyyw  +  °yyyywj  =  fattw  +  Nvxaxxw  f  2Nxyaxy‘w  *  Nyydyyw  (C20) 
Assuming  w  =  AW(x,y)cos(icont+ tf)  Equation  (C20), becomes 


n  )3  W  +  20  W  +  0  \v) 

u  {axxxxw  ^uxxvy  wyyyy  ■ 


!Vsxxw+2Nxy»xy*  + VyvWi 


(€21) 


The  determination  of  Nxx,  Nxy,  and  Nyy  can  be  accomplished  by  noting  that  Equation  (Cl  1) 
and  (Cl  3)  yield  the  following  relations 


rxx  ~  axx^>  €yy  ”  ayy^’ 


exv  ~ 


(C22) 


-M2 


The  forces  per  unit  length  may  be  computed  by  integrating  the  stress  over  the  thickness,  thus 


Nxx  = 


Nyy  = 


Nxy  ” 


c/2 

f  axx 
-c/2 

c/2 

/  V 

-c/2 

c/2 

J  * 

-c/2 


xy 


dz 


dz 


dz 


(C23a) 


(C23b) 


(C23c) 


Substituting  Equations  (C22)  into  the  generalized  Hooke’s  law  (Equation  (CIO))  and  per¬ 
forming  the  integrations  indicated  in  Equation  (C23)  the  forces  per  unit  length  are  expressed 
in  terms  of  the  strain  potential  function  as: 


Nxx  =  a*x^  +  vdyy*  ~  ( i+z;)  aT] 

Nyy  =  ^2  [3yy^  +  -  O+*0  «TI 

N  =Ec_a 
1Nxy  -|+p-°xy^ 


(C24a) 

(C24b) 

(C24c) 


Equation  (C21)  may  now  be  solved  for  W  in  tern  'f  the  frequency  and  potential  function  4/ 
for  given  boundary  conditions.  It  should  be  noteo  nat  W  in  Equation  (C21)  was  normalized 
with  respect  to  the  coefficient  A  which  must  be  included  to  obtain  the  actual  mid-plane  de¬ 
flection.  The  assumption  that  lateral  deflections  are  small  compared  to  the  thickness  allows 
us  to  neglect,  for  this  case,  the  stretching  f  the  mid-plane.  It  is  also  assumed  that  all  in¬ 
plane  loads  are  applied  to  the  mid-plane  oi  the  plate  and  there  are  no  thermal  gradients 
through  the  thickness  of  the  plate  (i.e.;  the  plate  is  initially  flat).  The  peak  stress  at  the  sur¬ 
face  due  to  dynamic  loading  and  inplane  forces  is: 


°xx  =  ’ 


yy 


Fr 

'20^2)[SXXw  +"3yy'vl 


aXy  -  2Gc  3XyW 


(C2Sa) 

(C25b) 

(C25c) 


If  the  iit-plane  forces  and  deflections  are  very  large,  such  that  the  mid-plane  is  stretched, 
second-order  terms  should  be  included  in  Equations  (C25).  These  equations  then  become 29 


Ec 


] 


°xx  =  "^2j  ^xxw  +  +  ^yy  w  +  ^dyw^ 


[dyyW  +  l/4(dy  w)2  +  vdxxw  +  |-(3xw)2] 


Ec 


oxy  -  2Gc  (3xyw  +  dx'*/9yW) 


(C26a) 

(C26b) 

(C26c) 
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APPENDIX  D 

COMPUTER  PROGRAMS  FOR  THERMAL  STRESS  CALCULATION 

Print-outs  of  fou»-  computer  programs  used  to  calculate,  subject  to  various  restrictions, 
thermal  stresses  in  panels  are  presented  in  this  appendix.  These  programs  are  written  in  the 
FORTRAN  language  and  in  a  manner  compatible  with  the  GE430  time  sharing  computer. 

The  four  programs  presented  here  can  be  merged  into  one  program  when  sufficient  storage 
capacity  is  available.  In  the  present  case  four  separate  programs,  each  within  the  capacity  of 
the  GE430,  are  used.  The  outputs  from  one  program  are  read  into  a  storage  file  and  later  re¬ 
called  as  inputs  to  subsequent  programs.  A  block  diagram  showing  the  basic  inputs  and  out¬ 
puts  of  the  four  programs  is  presented  in  Figure  42;  The  potential  program  is  presented  on 
pages  1 22  through  1 25 ,  the  strain  program  on  pages  1 26  through  1 28 ,  the  normal  stress  pro¬ 
gram  on  pages  129  through  130,  and  the  shear  stress  program  on  pages  131  and  132.  A  list 
of  the  symbols  used  in  these  programs  is  presented  on  pages  120  through  121  of  this  appendix. 

'iCi 

The  first  program  in  the  series  is  used  to  calculate  the  strain  potential.  In  all  subse¬ 
quent  statements  it  is  assumed  that  a  state  of  plane  stress  exists.  The  basic  equation  to  be 
solved  is 

9xx^  +  9yy^  =(1+^T  +  C  '  (Dl) 

where  ij/  is  the  strain  potential  function  and  C  is  a  constar*  of  integration.  Figure  43  shows 
the  model  used  in  computing  the  potential  function  by  the  finite  difference  technique.  31 

In  the  potential  program,  the  user  is  requested  to  supply  a  boundary  code  of  “1”  for  a 
free  boundary  and  “O’*  for  a  fixed  boundary.  The  user  is  then  asked  to  supply  five  coefficients 
of  tiie  temperature  distribution  in  the  x  direction.  The  form  of  this  distribution  is 

T(x,y)  =  A ,  +  A2  IX  i  +  A3  IX2 1  +  A4  IX3 1  +  A5  IX4 1  (D2a) 

-7.5  <X  <7.5,  -  IXI<Y<  IX I 

T(x,y)  =  Aj  +A2  IYI  +  A3  IY2I  +A4  IY3  I+A5  i  Y4 1  (D2b) 

-7.5  <  Y  <7.5,  IYI<  X  <  IYI 

In  lines  390  through  520  of  the  potential  program  the  two-dimensional  temperature  distribu¬ 
tion  is  computed.  This  distribution  has  the  desired  centerline  distribution  and  a  constant 
boundary  temperature.  Lines  80  through  100, 300  tlirough  320,  and  400  through  520  may 
be  modified  for  various  other  temperature  distributions. 

The  user  is  next  asked  to  supply  the  value  of  center  temperature,  Poisson's  ratio,  the 
coefficient  of  thermal  expansion,  the  length  of  the  increments  used  for  finite  difference 
calculations,  and  the  number  of  increments  in  each  direction.  Hie  dimension  statements 
may  be  modified  to  NX+1  by  NY+1  arrays.  The  user  is  then  asked  to  supply  a  code  which 
is  “1  ”  if  the  coefficients  have  been  normalized  to  A|  and  “0"  if  they  have  not.  This 
corresponds  to  a  temperature  distribution  which  is  normalized  to  the  center  temperature. 

To  evaluate  the  potential  function  it  is  necessary  to  determine  the  constant  of  integra¬ 
tion  in  Equation  (Dl ).  This  constant  is  evaluated  from  the  boundary  conditions.  Die  value 
of  tho  constants  for  free  and  fixed  panels  with  a  constant  boundary  temperature  arc 

(l  -  v) a (D3a) 
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FIGURE  42 

THERMAL  STRESS  PROGRAMS 
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FIGURE  43 


and 


-(1  +p)aTB 


(D3b) 


respectively.  If  boundary  conditions  other  than  fixed  or  free  are  to  be  used  the  value  of  the 
constant  C  (TPC,  lines  540  and  560)  must  be  computed  from  Equation  (D1 )  and  the  rela¬ 
tionship  between  strain  and  potential.  If  the  strain  at  a  point  is  known  then 


C  =  e 


xx 


a,b  +fHa,b 


(1  + v)  aT  I 


[a,b 


(D4) 


If  stress,  forces,  displacements,  etc.,  are  known  then  the  other  equations  of  elasticity  must 
be  applied  to  obtain  exx  I  .  and  e , 


a,b 


■yy  ia,b 


The  solution  of  the  differential  equation  will  be  restricted  to  a  bounded  plane  region  R 
with  boundary  S,  and  where  ^  is  defined  and  piecewise  continuous  on  S.  The  problem  is  to 
find  a  function  t//  (x,y)  which  is  continuous  on  R  +  S,  is  twice  differentiable  in  R,  and 
satisfies  in  R  the  linear  second  order  partial  differential  equation  (Equation  D1 ). 


Tire  finite  difference  representation  of  Equation  (D1 )  is: 

<MI,J  +  l)  +  iMI,J  -  1)  -  2^(1,  J)  .  *(I+1,J)+*(I  -1,J>~  2*(I,J) 


DY2 

=  (1  +  v)aT(l,  J)  +  C 
Solving  Equation  (D5)  for  \p  (I,  J)  we  obtain 


DX^ 


(D5) 


1 


0  (I* J)  =  Ti DX^TSy2)  |DY  { *(H. • J)  +  *(l  +  J)1  + 

DX2  {\ME  J-l )  +  ^(1,  J  +1)1  ~  DX2DY2((1  +M)aT(I,J)  +  Cl 


(D6) 


This  equation  may  be  solved  by  using  assumed  values  within  the  region  R  and  the  known 
values  at  the  boundary  S  for  the  function  ^(h  J)  on  the  right  side  of  Equation  (D6).  This  pro¬ 
cesses  is  then  repeated  by  replacing  the  old  values  within  the  region  R  by  the  new  computed 
values  of  \KJ.  J).  The  iu.ations  are  continued  until  the  old  and  new  values  of  \Sr(l,  1}  are 
sufficiently  close;  that  is,  until  ij/(l  I)  solves  Equation  (D  l )  and  the  boundary  conditions, 
with  sufficient  accuracy.  The  above  method  is  the  Jacobi  method.2-  An  extension  of  this 
method  is  the  Gauss-Siedel  method.  This  differs  from  the  Jacobi  method  in  that  the  new  values 
of  J)  for  adjacent  points  are  used  as  soon  as  they  are  computed,  litis  technique  increases 
the  convergence  rate  by  a  factor  of  two  over  the  Jacobi  method.  The  equation  utilizing  this 
method  is  given  in  tine:;  900  through  990. 

"it  order  to  compute  the  potential  function,  the  boundary  conditions  for  .nation  (D3) 
must  be  supplied.  In  the  case  of  a  free  panel  with  constant  boundary  temperature  the  traction 
forces  per  unit  area  are  zero  at  the  boundary  and  the  following  condition  for  stress  results: 


pxx  43  °yy  “  tfxy 


0 


(D7) 


From  the  generalized  Hooke's  law  and  tire  strain  displacement  relations  it  can  be  shown  that 
this  condition  reduces  to 


dxU  +  t>dv  =  dyV  +  udxu  =  { 1  ■  +  u) aT 


8yU  +  dxv«  0 


(DS> 

(D9) 


117 


The  simultaneous  solution  of  Equation  (D8)  and  the  application  of  the  shear  restriction 
(Equation  D9)  results  in  an  expression  for  the  displacements  at  the  boundary. 

u(x,  y)  =  aTgx  +  a  +  by  (DIO) 

v(x,  y)  =  aTBy  +  d  -  bx  (Dll) 

The  coefficients  a,  b,  and  d  represent  rigid  body  translations  and  rotations.  These  coefficients 
were  zero  for  the  cases  analyzed  in  this  study. 

When  the  expressions  for  the  strain  potential  are  substituted  into  Equatiorts  (DIO)  and 
(Dll)  and  the  results  are  integrated,  the  expressions  for  the  strain  potential  at  the  boundary 
are  as  follows: 

4>  I  =«TBy2/2  +  C,  (D12) 

lx  1  =  a 

=  aTBx2/2  +  C2  (D 1 3) 

ly  1=  b 

Since  the  potential  can  only  be  known  within  an  arbitrary  constant  and  this  constant  does 
not  enter  into  the  strain  or  stress  calculation,  we  assume  Cj  =  C2  =  0.  Thus  the  boundary 
conditions  for  the  free  panel,  used  in  the  potential  program,  are 

*  =oiTBy2/2,  *  =<*Tbx2/2  (D14) 

lxl=  a  lyl  =  b 

These  conditions  are  supplied  in  lines  720  and  760  through  790  of  the  potential  program. 

The  boundary  conditions  for  the  fixed  panel  are  u  *  v  =  0  at  the  boundary.  This  im¬ 
plies  that 

VBVsVa0  tl>15) 

*1  *ffy>  +  C{  .  (D16) 

lxi=  a 

I  =  f(x)  +  (D 17) 

ly  I «  b 

which  implies  f(x) s  f ( y >  =  0.  As  with  die  free  panel  we  may  choose  C  =  0  tints  #  =  0  at  the 
boundary.  The  boundary  conditions  for  the  fixed  panel  are  given  in  lines  740  and  760  through 
790. 

In  addition  to  the  use  of  the  Gauss-Siedel  method  of  successive  displacements  the  con¬ 
vergence  rate  was  increased  by  using  a  successive  over-relaxation  technique.33  The  coeffi¬ 
cient  for  relaxation  operations  used  in  this  program  is  0.443.  as  suggested  in  reference  33. 
When  different  increment  sizes  are  to  be  used  the  value  of  the  relaxation  factor  sltouid  bo 
changed.  The  over-relaxation  technique  is  employed  in  line  1030  of  *he  potential  program. 

Once  the  potential  Itas  been  calculated  with  sufficient  accuracy  it  is  placed  in  a  storage 
file  for  use  in  the  strain  program. 

The  strain  program  is  used  to  compute  bath  normal  and  shear  strains  front  the  po¬ 
tential  function.  The  strains  are  related  to  the  potential  function  by, the  following  equations: 

«xx  e  H’y  “  d>y**  exy  =  ^  CDW 

The  finite  difference  approximations  of  the  above  equations  are 
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(D 19) 
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”  '  DY2 
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(D20) 


(D21) 


DXDY 


These  equations  are  solved  in  lines  670  through  720  in  the  strain  program. 

The  solution  of  Equations  (D!9)  through  (D21)  can  only  be  found  within  the  region 
R.  In  the  program  for  free  panels  the  strain  at  the  boundaries  is  known  and  thus  may  be 
supplied  directly.  For  the  fixed  panel  the  normal  and  shear  strains  at  the  boundary  are  zero 
corresponding  to  zero  deflection  at  the  boundary.  From  Equations  (D12)  and  (D13)  it  can 
be  shown  that  for  the  free  panel. 


-e. 


xx  !b  yy  la 


=  aT 
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exy  Ib 


-  0 


(D22> 

(D23) 


The  boundary  conditions  for  the  fixed  and  free  panel  are  supplied  in  lines  340  through  480 
and  5 1 0  through  650,  respectively  .  The  values  of  normal  and  shear  strain  are  then  stored  in 
a  file  for  use  in  subsequent  programs. 

Tlie  normal  stress  program  applies  the  generalized  Hooke's  law  to  obtain  the  stress  dis¬ 
tribution  from  the  normal  strains  and  the  temperature.  The  relations  for  calculating  normal 
stress  are 

E  , 

ffx*sCU,2)  1^\+^yy“Cl+^)  (D24) 
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°y yji^) 


%  !«yy  +  Wxx  .**  ^  t 


These,  equations  are  solved  in  lines  4 1 0  through  470  in  the  program.  A  routine  to  compute  the 
maximum  normal  stresses  in  the  x  and  y  directions  is  given  in  lines  630  through  810.  If  the 
value  of  maximum  stress  occurs  at  more  than  one  point  only  the  first  point  is  retained. 

The  shear  stress  program  uses  the  Hooke's  law  relation  between  shear  stress  and  shear 
strain  sliown  below: 


"sy  =  0txy 


026) 


Tlie  equation  is  solved  in  line  220  of  the  program.  A  routine  for  calculating  the  maximum 
shear  stress  is  given  in  lines  390  through  5 10.  This  routine  is  subject  to  the  same  restriction 
as  tlie  maximui  normal  stress  routine. 

Tim  output  of  a  typical  computer  run  is  presented  on  pages  148  through  1 56.  Responses 
which  the  user  types  into  the  computer  are  underlined.  The  example  analysis  is  for  a  15  x  15 
inch  steel  panel  with  free  boundaries  and  a  temperature  distribution  with  TtX)  of  the  form, 

T(X)«  800  -  13.33X  -  8.88X*-  (D27) 

Tin's  corresponds  to  a  880°1;  center  temperature  and  280°F  boundary  temperature,  since  the 
temperatures  used  in  the  program  arc  referenced  to  am  bkmt. 


(The  following  symbols  appear  first  in  the  potential  program) 
SYMBOL  MEANING 


XCI) 

Y(I) 

A(I) 

PfU) 

POT(l,J) 

D1F(I,J) 

T(1,J) 

TP(U) 

NU 

MAX 

BOUND 

TC 

AL 

XL 

YL 

NX 

NY 

DX 

DY 

NXI 

NYI 

Nxe 

NYC 

DX2 

ov: 

RE  fi 
CO 

m: 


CHECK 
TEST  2 
DIRECT 
WEIL 
TfcMPF 


coordinate  in  x  direction 

coordinate  in  y  direction 

coefficients  of  temperature  distribution 

old  value  of  potential  function 

new  value  of  potential  function 

difference  between  new  and  old  potential  function 

temperature  distribution 

nonhomogeneous  term  in  potential  equation 

Poisson’s  ratio 

maximum  percentage  value  of  DiF(I»J) 

code  identifying  boundary  type  » 

center  temperature 

coefficient  of  thermal  expansion 

length  in  x  direction 

length  in  y  direction 

number  of  increments  in  x  direction 

number  of  increments  in  y  direction 

element  length  in  x  direction  for  finite  difference  calculations 

element  length  in  y  direction  for  finite. difference  calculations 

number  of  points  in  the  x  direction 

number  of  points  in  the  y  direction 

x  index  value  at  the  center 

y  index  value  at  the  center 

square  of  DX 

square  of  DY 

code  to  identify  whether  the  Aft)  m  normaliaed  to  A(l)  or  not 
coefficient  tor  relaxation  operations 
constant  of  integration  in  potential  equation 
counter,  determined  number  of  integrations  per  juries 
sequence  code  to  obtain  another  series  ofitCfatiOfts 
printing  code 

sequence  'code  to  obtain  another  series  of  iu».; ,ons 
print  code,  used  to  file  data  or  not 
file  containing  final  values  of  potential  function 
file  containing  values  of  temperature  distribution 


(The  f blowing  symbols  appear  first  in  the  strain  program) 


SYMBOLS 


MEANING 


HU) 

values  of  potential  function  from  P01TIL 

EPSXX(U) 

strain  in  x  direction 

EPSYY(U) 

strain  in  y  direction 

liPSXY(I,J) 

shear  strain 

DELX 

same  as  DX 

DULY 

same  as  DY 

li 

Young’s  modulus 

DXY 

one  half  product  of  DX  and  DY 

120 


SYMBOLS 


MEANING 


\ 

\ 

v 

EPSXXF  file  containing  values  of  strain  in  x  direction 

EPSYYF  file  containing  values  of  strain  in  y  direction 

EPSXYF  file  containing  values  of  shear  strain 

(The  'ollowing  symbols  appear  first  in  the  stress  progum) 

SYMBOLS  MEANING 

NXP  same  as  NX1  in  potential  program 

NYP  same  as  NY  1  in  potential  program 

NX  1  value  of  x  coordinate  of  the  maximum  stress  in  the  x  direction 

NYl  value  of  y  coordinate  of  the  maximum  stress  in  the  x  direction 

NX2CH  value  of  x  coordinate  of  the  maximum  stress  in  the  y  direction 

NY  XK  value  of  y  coordinate  of  the  maximum  stress  in  the  y  direction 

SIGXX(LJ)  stress  in  the  x  direction 

SIGYYUJ)  stress  in  the  y  direction 

ALPHA  same  as  AL 

NX1  dummy  x  coordinate  ?o  compute  location  of  maximum  stress  in  x  direction 

NY  I  dummy  y  coordinate  to  compute  location  of  maximum  stress  in  x  direction 

NX2  dummy  x  coordinate  to  compute  location  of  maximum  stress  m  y  direction 

NY2  dummy  y  coordinate  to  compute  location  of  maximum  stress  in  y  direction 

CK  VAL 1  dummy  variable  to  compute  maximum  stress  in  x  direction 

CKVAL2  dummy  variable  to  compute  maximum  stress  in  y  direction 

TEMPT  1  temperature  at  paint  of  maximum  stress  in  x  direction 

TEMPT 2  temperature  at  point  of  maximum  stress  in  y  directum 

(The  following  symbols  appear  first  in  the  shear  program) 

SYMBOLS  MEANING 

SSGXYCUl  drear  stress 

SX1  dummy  x  coordinate  to  compute  maximuin  drear  stress, 

SYl  dummy  y  coordinate  to  compute  maximum  shear  Stress 

CKV^L  dummy  variable  tired  to  compute  maximum  shear  stress 

SXE  value  of  x  coordinate  of  maximum  shear  stress 

SYE  value  of  y  coordinate  of  maximum  shear  xtress 


P0TENT I  Ai- 


IO  DIMENSION  X(17),Y{17),A<5) 

20  DIMENSION  P(  1  7,  1  7>,P3TC  1 7,  1  7  >,  DIF(  1  7,  1  7  ) 

SC  DIMENSION  TCI  7, 17) 

40  DIMENSION  TPC17, 17) 

50  PEAL  NU,MAX 

60  PRINT,"  FOR  FREE  TYPE  1,  FOR  FIXED  TYPE  0" 

70  READ,  BOUND 

80  PRINT,  "SUPPLY  THE  5  COEFFICIENTS  FOR  TCX>." 

90  PRINT, 

100  READ,  AC1  ),AC2),AC3),AC4),A<5) 

110  PRINT, 

120  PRINT,"  SUPPLY  CENTER  TEMP,,  POISSON'S  RATIO,  COEF .  OF  THERMAL 
130  PR I NT, "EXPANSION,  LENGTH  IN  X  AND  Y  DIRECTION,  NUMBER  OF" 

140  PRINT,"  INCREMENTS  IN  X  AND  Y  DIRECTION." 

150  PRINT, 

160  READ,  TC,  NU,  AL,  XL,  YL,  NX,  NY 

170  PRINT,  "IF  THE  TEMPERATURE  COEFFICIENTS  YOU  SUPPLIED  UERE" 

180  PRINT,  "NORMALISED  TO  CENTER  TEMP,  PRINT  1,  IF  NOT,  0" 

190  DX  SXL/NX 
200  DY*YL/NY 
210  NX  1 *NX* 1 
220  NYl»NY+t 
230  NXC»  NX/2+1 
240  NYCaNY/2+1 
250  DX2*DX*UX 
260  0Y2»DY*DY 
270  READ,  FEB 
280  C0«*443 

290  IF  CREB.EQ.O.)  G0  TO  9 
300  DC  10  K  =  1 , 5 
310  A(X)»A(-K)*TC 
320  10  CCNTINUE 
330  9  PRINT, 

340  DC  l  1  **  1 ,  NX  1 
350  DC  I  ! , NYl 

360  X<  1  )*-<XL/2»  )♦< I~  1  )*DX 
370  YCJ)e-CYL/g.  )♦<,»-! T*DY 
3S0  l  CONTINUE 
390  DC  14  I  «MX  C,  NX  1 
400  DC  14  J •NYC* NY  1 

410  T<  I, J)*A(  1  )+A(2>*X{  I  )+A(3>«<XU  )*>2  >  +  A(4)*CX(  I  >**3> 

4204  +A  ( 5  ) MXf I >**4» 

430  M*<a*NXC)-J 
440  T(1,H)*T(I, J) 

450  K*C2*NYC) - 1 
460  ?(X,J)«TCI,U) 

470  T(K,H>»T(1,J> 

480  T(J,I)«T(I,J) 

490  T<N,I)«T<!,J> 

500  T<J,K>»TU,U> 


•  ws  w*.*4 


l 

8 


Yf 


1: 

i 
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.? 

l! 


POTENTIAL  CONTINUED 


530  T(M,K>=T(i,J> 

520  14  CONTINUE 

530  IF (BOUND. FQ.O)  G0  T0  738 

540  TPC=- (AL*T ( 1, 1 >*(1 . -NU) ) *(DY2*DX2/<2 . * (DX2+DY2 ) ) > 

S50  G0  T0  739 

560  738  TPC=(AL*T(1, J >*  (  1 +l\IU)  >*  (DY2*DX2/(2  .  *(DY2+DX2  >  >  ) 

570  739  CONTINUE 
580  D0  2  1=1, NX! 

590  D0  2  J= 1 , NY 1 

600  TP( I, J>=-(DX2*DY2*( 1 . +NU ) *AL*T ( I , J > )/ (2 .* (DX2+DY2 ) ) 

610  TP < I , J) =TP( I, J >  +TPC 
620  PCI, J)=TP(I, J) 

630  2  CONTINUE 

640  15  CONTINUE 

650  IK=0 

660  7  CONTINUE 

670  MAX=u 

680  675  CONTINUE 

690  P0  88  1=1, NX 1 

700  D0  88  J=  1 , NY  1 

710  IF (BOUND. EQ .0)  G0  T0741 

720  POT ( NX  J , J ) =AL*T (NX  l ,  J) *Y(  J) *Y(«J )  /2  • 

730  GO  TO  742 

740  741  POT  C 1 7>K) =0 « 

750  7<.2  CONTINUE 

760  P0T(1,J)=POTCNX1,J) 

770  POT ( I  ,  NY  1 ) =POT (NX  1 , J) 

780  P0T ( i ,  1 ) =POT ( NX  1 , J  > 

790  88  CONTINUE 
800  POT  (NX i  ,  J )  =>0 . 

810  DD=BX2/ (2 • * ( DX2+DY2 ) ) 

820  DDY5sDY2/  <  2 .  *  (  DX2  +  DY2  )  ) 

830  NX2=NX- 1 
840  NY2  =NY- 1 
850  D0  20  1=2, NX2 
860  D0  20  «J* 2,NY2 

870  P0T<X,J>»(  DDY*(P0T(I- !,U)+P(T+1, J)))+CDD*(P0T(I,J-1 ) 

880&  +  PC  I,  <J+ 1  )))}  +  TP{I,J) 

890  20  CONTINUE 
900  m  24  <J=2, NY2 

910  POT (NX, J)= ( ( DDY#(P0T(NX2, J) +P0T(NX 1 , J) > )+ ( DD* (POT (NX, J  -  1  )  + 
9204P(NX *  J+ 1 ) ) ) )+TP(NX,  J) 

930  24  CONTINUE 
940  D0  26  I=2,NX2 

950  POT  C I, NY>*  <  C DDY*(P0T  < I , NY2 )+P0T (I, NY !)))  +  ( DD* (POT ( X- I , NY)+ 
960&P( 1  + 1 , NY) ) >  >  +  TP( i ,NY) 

970  26  CONTINUE 

980  POT  (NX, NY )  =  ( •' DDY*  (P0T(NX2,  NY>  +  P0T  (NX  1 ,  NY )  > )  + 

9906  (DD*(P0T (NX,  NY2  )  +  ^0T  ( NX, NY1  )  > )  >+TP(NX,NY) 

1000  00  49  I « J , NX  1 
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1010  C0  49  J=1,MY!  - 

1020  DIF <  I,  J>=P0T(  I,<J)-P(  1 ,  J) 

1030  PCI,.<J>*(  1  .+C0)*P0T(Z. J>-C0*P(I, J) 

1040  IF  <ABS(  (DIFU,J>/P0T(I,  J)  > *  I  00 . > .LE.MAX)  G0  T0  49 
1050  IF  (P0T(  I,J>.EQ.Q»)60  T0  49 
1060  MAX«ABS<  CDIFU,  J>/P0TCI,J>>*10O> 

1070  49  CONTINUE 

1080  PRINT  2001, IK, MAX, DIFC NXC, NYC) ,P0? (NXC, NYC >,P0T < NX! , NY I > 
1090  2001  FORMAT  <3X,  16,  4E1  1  3) 

1100  PRINT, 

1110  IKaIK+i 

1120  IF < I K, LE.2G )  G0  T0  7  : 

1130  PRINT  2000,  MAX  ... 

1140  PRINT,”  SUPPLY  STEP  SIZE  F0R  PRINT  STATEMENTS  F0R  I  AND  J” 
1150  READ,  INCRX, INCRY 

1160  2000  FORMAT <”TH£  MAXIMUM  PERCENTAGE  DIFERENCE.  IS  ”,E1C.3) 
1170  PRINT,”  IF  MORE  “ITERATIONS  ARE  DESIRED  TYPE  1  ,IF  NOT  0” ' 
11.80  READ,  TEST 
1190  IF- < TEST.EQ.OJGC  T0  15  . 

1200  38  CONTINUE 

1210  PRINT,”  IF  0NLY  THE  CENTERLINE  VALUES  ARE  DESIRED  TYPE  t>” 
1220  PRINT,”  OTHERWISE  TYPE  0” 

1230  READ,  CHECK 

1240  PRINT,”  X  Y  POT" 

1250  IF (CHECK. EQ. 1 )  G0  T0  95 
1260  D0  78  1  =  1 , NX  1, INCRX 
1270  D0  78  J=!,N Yl,  INCRY 
1280  PRINT  79,  X ( I) , YC J ) ,P(I , J > 

1290  79  FORMAT  <5X, F4 . 1 , 5X, F4 . \ , 5X, E 1 2 . 5 ) 

1300  78  CONTINUE 

1310  95  CONTINUE 

1320  DO  85  I- 1 , NX  1 , INCRX 

1330  PRINT  79,X(I), YCNYC ) , P( I , NYC ) 

1340  85  CONTINUE 

1350  D0  90  J“l, NY1 ,  INCRY 

1360  PRINT  79,XCNXO,  Y(U> • P(NXC, J) 

1370  90  CONTINUE 

1380  PRINT  79,X(NXC),Y(NYC>,P(NXC,NYC) 

1390  PRINT, 

1400  PRINT, 

1410  PRINT, 

1420  PRINT, "IF  ANOTHER  ITERATION  IS  DESIRED  TYPE  01  OTHERWISE  1 
1430  READ, TEST2 

1440  IF(TEST2.EQ.O)  G0  T0  15  .  '■ 

1450  PRINT,  ”IF  IT  IS  DESIRED  TO  WRITE  INTO  A  FILE,  TYPE  0  IF” 
1460  PRINT,  ”N0T,  TYPE  0.  (THIS  IS  F0R  B0TH  P0TFIL  .&  TE'MPF)". 
1470  READ,.  DIRECT  ...  ' 

1480  IF  (DIRECT. ESUO)  G0  T0  71  v-‘“ 

1490  CALL  SPENF  ( 1 ,  "P0TFIL”,  7 )  ..  .  . 

1500  D0  72  1=1, NX1  ' 

1 24  '  '  '  ' 


POTENTIAL  CONTINUED 


1510  D0  72  J=l,NYi 
1520  WRITE  (U70>  P<I,J> 

1530  70  FORMAT  (E14.7) 

1540  72  CONTINUE 

1550  CALL  0PENF  <4, "TEMPF", 7) 

1560  D0  3  K= 1 , NX  1 
1570  D0  3  L* 1  *  NY 1 
1580  WRITE  <4J4>  TCK, L) 

1590  4  FORMAT  <F15.5) 

1600  3  CONTINUE 
1610  71  CONTINUE 
1620  PRINT, 

1630  PRINT,  *'  X  Y  TE 

1640  PRINT, 

1650  00  5  I « 1 ,  NX  1 ,  I NCftX 
1660  D0  5  0*1, NY1, INCRY 
1670  PRINT  6,  XCD,YCJ>,T(I>J) 

1680  6  FORMAT <  3X, F4 » I , 5X,F4» 1 , 5X, E22 . 5) 
1690  5  CONTINUE  '  • /. 

1700  END  ■■ 


STRAIN 


JO  REAL  MO 

20  DIMENSION  P<l7,i7)  .  -  -  ..  ■ 

36  DIMEM5H0N  EPSXX C  17,  17 ) , EPSYYC  !  7,  1 7  ) ' 

40  DIMENSION  EPSXYU 7,  17>,T*T7,  17>  ;> 

50  DIMENSION  X(J7)„Y<J7> 

60  PR1NT,"SUPPLY  Y0UN&I 5  MODULUS.* >P0 1 SS0N * S  RATIO,  C0EF..0P 
70  PR  I NY,  '‘THERMAL  EXPANSION,  i  LENGTH  IN  X  AND  Y  DIRECTION,” 
80  PRINT, "NUMBER  OF  INCREMENTS  IN  X  AND  Y  DIRECTION*** 

90  READ,  E*NU,  AL,  XL,  Yt.,.NX,  NY 
100  D£LX»XL/NX  -  ;  = 

110  DELY»YZi/NY 
120  NXi*NXU 
130  NYlsNY+i 

140  CALL  0PENF  <1 >"P0YFILM>  7> 

.150  CALL  E0FTST  <1,M) 

160  DO  *  I»1,NX1 
170  D0  ?.  *J« 1,MYI 
180  READ  (113)  P(LJ) 

190  3  FORMAT  CE14.7) 

200  2  CONTINUE 

210  CALL  0PENF  <4,"TEMPF",7> 

220  CALL  E0FTST  <4,M> 

230  DO  57  I«l,NXl 
240  W  57  J=i,NYl 
250  READ  <4i59)  T(I,J> 

260  59  FORMAT- (FI 5. 5) 

270  57  CONTINUE 

280  PRINT, "PRINT  1  F0R  FREE  BOUNDARY,  0  FOR  FIXED  BOUNDARY. 
290  DX2  =*DELX*DELX 
300  DY2 = DEL Y* DELY 
310  DXYaDELX*DELY* • 5 
320  READ,  B0UND 
330  IF  (B0UND*EQ* 1 )  GO  TO  53 
340  D0  4  K* 1 , NX  1 
350  DO  4  L* 1 , NY 1 
360  EPSXXCK, 1 )a0* 

370  EPSXXC l,L)®0s 
380  EPSXX (K, NY 1 )»0. 

390  EPSXX<NX1,D«0. 

400  EPSYYCK, 1 >»0* 

410  EPSYY < i,L)“0» 

420  EPSYY (K,NY1)®0. 

430  EPSYY(NXI,L>*0* 

440  EPSXYtK, 1 )*0* 

450  EPSXY< 1 ,L)«0. 

460  EPSXYCK,  NY1 )*0. 

470  £PSXY<NXl,L>*0* 

480  4  CONTINUE 
490  00  TO  54 
500  53  CONTINUE 
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510  DO  61  K= 1 / NX l 

520  00  61  L*1,NY1  . 

530  £PSXX< K,  1>*AL*T<K,  1)  "  - 

540  £P$XX<1,L)*AL*T<1,L> 

550  EPSXXOONX1 )«AL*T(K>NX1 > 

560  EPSXX<NX1,L>=AL*T<NXI,L> 

570  EPSYY<K,NX1)*AL*T(K,NX1> 

580  EPSYY (  1  *  L  ) “AL*T (  1 ,» L  ) 

590  EPSYY  <  NX  1 ,» L  ) »  AL*T (  NX  1  >  L  ) 

600  EPSYYC17,K>=EPSXX(17,K) 

610  EPSXY(  1**.'**0« 

620  EPSXY(NX1,L>»0. 

530  EPSXYCK, 1)=0. 

640  EPSXY(K,NX1 >=0. 

650  61  CONTINUE 

660  54  CONTINUE  ...  - 

670  00  5  I«2*NX 
680  00  5  J*2,NY 

690  EPSXX<I,<n«<P(I+l, J)«2*P{I# J)+P < I- 1,J)>/<DELX**2> 

700  EPSYYC 1/ J>* <P< I>J+1>-2*P(I>J)+P(I/J-1) }/(DELY**2) 

710  EPSXY< I*  J)“<P<  I_»«J>+P(I  +  i*J+l)“P(I*J+l>-P<I  +  l,»J>  >/DXY 

720  5  CONTINUE 

730  NXC«NX/2+l 

740  NYONY/2+1 

750  PRINT  8, 

760  00  13  K»1,NY1 
770  X<K)-0. 

780  YOO»(<K-l>*DELY>-CYL/2> 

790  PRINT  9*X(NXCj^Y(K)#EPSXXCNXCjrK> > EPSYY <NXC^K>*EPSXYCNXC^K> 
800  13  CONTINUE 
810  00  12  K» l *  NX 1 
820  Y(K>»0. 

830  X<K)=<<K~l)*DELX>-<XL/2> 

840  PRINT  9>X(K).»Y<NYC),»  EPSXX<K* NYC )# EPSYY <K# NYC >/EPSXY(K» NYC) 
850  12  CONTINUE 

860  8  FORMAT  <3X,MX".,6X*"Y",6X,"EPSXX",6X,  "EPSYY",  6X,*EPSXY"J 
870  9  F0RMAT<IX,F4.1,3X,F4.1,3X,E1O.3,2X,E1O.3>2X,E10*3> 

880  CALL  0PENF  <2, "EPSXXF", 7) 

890  00  14  1*1, NX1 
900  D0  14  J-l,NYl 
910  WRITE  <2J  3)  EPSXX(I,<J> 

920  14  CONTINUE 

930  CALL  0PENF  <3,"EPSYYF", 7) 

940  00  15  I«l,NXi 
950  D0  15  J*1,NY1 
960  WRITE  <31 3)  EPSYY{i,J) 

970  IS  CONTINUE 

980  CALL  0PENF  <S,"£FSXYF", ?> 

990  DO  710  I“1,NX1 
1000  00  710  J»1,NY1 
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1010  WRITE (5;3) EPSXY  <  I  *  J ) 

1020  710  CONTINUE 

1030  PRINT,»E PSXX,  EPSYY,  AND  EPSXY  HAVE  BEEN  FIEED" 
1040  6  CONTINUE 
1050  END 


STRESS 


iO  REAL  NU, NX 1 CH,  NX2CH,  NY I CH,  NY2CH 
20  DIMENSION  E?SXX<  i  1,  17),  EPSYYC  1 7,  i  7>,  T<  1 7,  1?> 

30  DIMENSION  SXGXXU7,  17>,SIGYY(17,  17) 

40  DIMENSION  XU7>,Y<17> 

50  PRINT, 

60  PRINT,"  SUPPLY  YOUNG'S  MODULUS,  POISSON'S  RATIO,  COEF." 

70  PRINT, "OF  THERMAL  EXPANSION,  LENGTH  IN  X  AND  Y  DIRECTION," 
80  PRINT, "NUMBER  OF  INCREMENTS  IN  X  AND  Y  DIRECTION." 

90  READ, E, NU, ALPHA, XL, YL, NX, NY 

100  DELX=XL/NX 

110  DELY*=YL/NY 

120  NXP=NX+ 1 

130  NYP=NY+ ’ 

140  NXCeNX/2+1 
150  NYC=NY/2+l 

160  CALL  OPENF  (2, "EP SXXF", 7) 
i70  CALL  E0FTST  <2,N> 

180  D0  10  I»1,NXP 
190  DO  10  Jai,NYP 
200  READ(2;23)EPSXX(I,J) 

210  23  F0RMATCE14.7) 

220  10  CONTINUE 

230  CALL  0PENF  (3, "ErSYYF", 7) 

240  CALL  E0FTST  <3,N) 

250  D0  15  I«l,NXP 
260  D0  15  J*i,NYP 
270  READC3J23)EPSYY(I, J> 

280  15  CONTINUE 
290  PRINT, 

300  17  FORMAT  <3X,"X",5X,"Y",6X,"TEMP",6X,"SI6XX",6X,"SIGYY"> 
310  PRINT, 

390  CALL  0PENF  <4, "TEMPF", 7) 

330  CALL  E0FTST  <  4, N) 

340  D0  21  1*1, NXP 
350  DO  21  J»I,NVP 
360  READ  <41 22)  T<I,J> 

370  22  FORMAT  CF15.5) 

380  21  CONTINUE 

390  19F0RMAK  1X,F4, 1,3X,  F4.  1,6X,  F7. 1,3X,E10.3,3X,E10.3) 

400  PRINT, 

410  00  16  I ■  l, NXP 
420  DO  16  J*1,NYF 

430  SlGXX(J„J>«E*fEPSXX<I,J>+NU*(EPSYYU,J)>-Cl+NU>*ALPKA*T 
440*  ( I,U))/<  1-NU**2) 

450  SIGYY(I,  J)BE*(EPSYY<  I,U)4-NU*CEPSXX<  I,  J>  >-<  l+NU)*ALPHA*T 
460*<I, J))/<S~NU**2) 

470  16  CONTINUE 

480  PRINT, "SIGXX  AND  SICYY  HAVE  BEEN  COMPUTED." 

490  PRINT, 

500  PRINT  17, 


129 


STRESS  CONTINUED 


510  PRINT, 

520  D0  57  1=1, NXP 
530  DO  57  J=j,NYP 
540  X<I)*<I-l)*DELX-<XL/2.> 

550  Y<J)*<J-l)*DELY-CYL/2.) 

560  57  CONTINUE 
570  DC  .  K=1,NYP 

580  PR  INI ;  %X<NXC>,Y(K>,T<NXC,K),SXGXX<NXC,10,SIGYY<NXC,X> 

590  18  CONTINUE 
600  DC  20  K«1,NXP 

610  PRINT19,XCK>, YCNYC),T<K, NYC), SI6XXCK,NYC),5IGYY(K,NYC) 

620  20  CONTINUE 
630  NX1-1 
640  NY1 = 1 
650  NX2*1 
660  NY2»1 
670  CKVAL1-0. 

680  CKVAL2=Q.  * 

690  DC  .29  X“l,  NXP 
700  DC  29  J«1,NYP 

710  IF  (ABS( 5 JGXX( J, J) ) *LE*  CKVAL1 }  60  T0  30 
720  GKVAL1 “ABS< SIGXX< I, J) ) 

730  NX1=X 
740  NY1 *J 
750  30  CONTINUE 

760  IF  <ABS<SIGYY<  I, U)>«L£.CKVAL2)  60  T0  31 
770  CKVAL2»ABS<SIGYYU,J>> 

780  NX2«I 
790  NY2=J 
000  31  CONTINUE 
810  29  CONTINUE 
820  PRINT, 

830  PRINT,  -MAXIMUM  VALUES  FOR  STRESS  ARE  I** 

840  PRINT, 

850  PRINT,  "  X  Y  SI6XX  SI6YY  TEMP- 

860  PRINT, 

870  TENPT 1  *>T  <NX1,NY1  ) 

880  TEHPT2=T<NX2,NYS) 

890  NXlCH*<<NXl-I)*DELX)«<XL/2.) 

900  NY1 CH*<CNYl -l >*DELY)« (YL/2» ) 

910  NX2CH=<  <NX2-  1  )*DELX) -  (XL/2.  ) 

920  NY2CM«< (NY2-1 )*DELY)-(YL/2«  > 

930  PxINT  32,  NX1CH, NY1CH, SI6XX<NX1,NY1 ),SIGYY<NX1,NY1 ), TEMPT 1 
940  PRINT  32,  NX2CH,NY2CH,SIGXXCNX2,NY2),SXGYY<NX£,NY2>,TEMPT2 
950  32  FORMAT  <  1X,F4*  1, 3X,F4. 1,2X,  310.3,2X,E10,3,2X,F6 a  > 

96 0  28  CONTINUE 
970  ENT' 


SHEAR 


10  DIMENSION  EPSXYU7,  i  7 >, SIGXY< 1 7,  17),X<m,Y<  17) 

20  INTEGER  SX1,$Y! 

30  REAL  NU 

40  PR  1 NT ^ "SUPPLY  E,  NU,  THE  NUMBER  0F  INCREMENTS  IN  THE  X  AND  Y" 
50  PR INT, "DIRECTIONS,  AND  PANEL  LENGTH  ALONG  X  AND  Y" 

60  READ,  E,NU,NX,NY,XL,YL 
70  NX  1 =NX+ 1 
80  NY  1  =»  NY+i 

90  NXC*NX/2+ 1  .  V/ 

100  NYC=NY/2+  t 

110  DX=XL/NX 

120  DY°YL/NY 

130  PRINT  201, 

140  201  FORMAT < 6X,"X",9X,"Y",8X. "STRAIN", 8X, "SHEAR”) 

150  CALL  0P£NF(5, "EPSXYF", 7) 

160  CALL  E0FTSTC5, N) 

170  G=E/(2.*<  l  .+NUXT 
180  D0  19  1=1, NX! 

190  D0  19 

200  READ  (5,*23)  EPSXY(I,J) 

210  23  FORMAT <£14, 7> 

820  SIGXYU,  J)*fi*EPSXYU,U> 

230  19  CONTINUE 
240  DO  SS  1*1, NX! 

250  X{1-}*<<X-I>*DX)*<XL>8.} 

260  21  CONTINUE 

270  DO  22  U«1,NY1 

280  Y(J)»<  <U“  1  >.*DY>«<Yl,/2,  > 

290  22  CONTINUE 

300  PRINT, "SUPPLY  PRINT  INCREMENT  FOR  X  AND  Y" 

310  READ,  INCRX,  INCRY 

320  D0  27  K«1 » NY 1, INCRY  — 

330  PRINT  200,Xf  NXC>,  Y  (K),  EPSXY (NXC,K).,SIGXY CNXC,X> 

340  27  CONTINUE 

350  DO  29  K«1 ,NXl, INCRX 

360  PRINT  2Q0,X(K),Y< NYC > , EPSXY 4 X, NY£  V,  S I GXY< K, NYC) 

370  29  CONTINUE 

380  £00  FORMAT <3X,fr'7.2, 3X, F7, 2, 3X,EiO» 3, 3X, El  0,3  > 

390  3X1*1  . 

400  SY 1  *  1  * 

410  CKVAL*0. 

420  D0  290  1*1, MX! 

430  D0  290  «J  « 1  ,  NY  1 

440  lF(A8S<St3XYU,U)),LE<CKVAL>  G0  TS  30 
450  CKVAL=A8S<SIGXY(I,U) > 

460  SXl»! 

470  SY1-U 

480  30  CONTINUE 

490  290  CONTINUE 

500  SXF»( (SX1 - 1 . ) *DX)  - (XL/2 • ) 
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SHEAR  CONTINUED 


510  SYF“  ( < SYl - l • )*DY) - <YL/2 . > 

520  PRINT  2  70, 

530  270  FORMAT <3X,‘*X",  7X,*,Y",2X,"  MAX  SHEAR  .  MAX  STRESS**) 

40  PRINT271,SXF,SYF,EPSXY<SXI,SYi >, SIGXY(SX1,SY1 ) 

550  271  FORMAT < 3X,F4* 1 , 3X>  F4« 1,  3X,E10»  3, 3X, EiO»  3) 

560  END 
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